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FREQUENCY ANALYSIS OF BEAM 
GIRDER FLOORS 


By 1M. ASCE 


JA procedure for —_— the frequencies and | modes of such | three-dimensional - 
— 


‘caample, to building f frames or to airplanes. can 
be used with \ advantage f for the frequency ana analysis of of any ‘complicated structure, 


and are easier to determine. If n necessary, t ‘this 
process ca can be used supantetiy. a: he principle i is theoretically exact but leads 
: to the use » of infinite : series; ; for practical p purposes excellent approximations are 


. obtained by using only the first few terms of these series. _ - A reader who wishes 7 
| to > apply the method should have some knowledge of the theory of oscillations 7 
and, in. particular, of the fundamental notion on of ‘normal modes” and of their 


_ ‘The paper is is ‘presented in in four sections including 2 a numerical example 7 


‘(Gection 4). 7 Section 1 deals with the theory and its ts application tc to statically 7 


‘determinate systems. ‘Section 2 an extensi sion of the theory to 


statically indeterminate systems. Section Bisa collection of methods 
formulas for reference in determining frequencies and modes of stant systems. 


“Although some of these formulas ar are known, many | of them a are new. The 
Purpose of Section 4 is to serve as a guide for the arrangement of practical ¢ com- 
-Putations. For this reason the simplest, sy system on n which all the steps of the 
method can be demonstrated was chosen, thus enabling the: writer to 


3 every ‘detail fully. ‘The method was s devised (and has been applied successfully) 
ease; NOTE —Written comments are invited for immediate publication; to insure Publication the last dis- 
 Cussion should be submitted by March 1, 1950. 


1 Associate Engr., and Hanover, Cons. Engrs., New N. Y.; Lecturer in Civ. Eng., 
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in Section 4. 


discussions, ‘conform with the American Standard Letter Symbols 
for Structural Analysis, prepared by a Co Committee of the American Standards =e 
_ Association with ASCE participation. . An important exception is that the 7 
symbols m and M i in this paper are used for masses and not for bending mo- 
ments. The latter are designated by mand M. defined where in 


in 
Pay first appear, in the text or by il illustration. The most significant sym-_ a 

eal 


c= = an on factor; a coefficient; c; = a specific value corresponding 


t= = a number; 
j=anumber; 
k = a number; a panel, or section n point number on a beam; — _ 
_M = moment corresponding to an externally applied load; Eq 
m= value of M for a unit load; we 


M= concentrated m reactions; 
m = distributed masses; 
ch 
deflection: 
= = deflection in a 
ux = value of u due to redundant moments ae * 
) = unit weight, in pounds p per foot; also elastic weights; ; 
th 
X = redundant moments; 
Ke 
y= amplitude of a mode: |e 
= successive approximations of y;_ 
= = displacement, afunctionoft; 


er - a substitution constant in Eqs. 26 and 41; 
A= determinant; 


= = amplitude o of the normal mode, ¢; corresponding to wi; er 
‘$s = amplitudes for systems A and B, respectively; 


” ‘deflection due to redundant moment X = 1; and ae 


circular frequency: 
7 wa and wg = frequencies for systems A and B, respectively; 
and = specific v values of w. 


 Pundemental Facts Concerning Oscillations.—If any mechanical system is 


a in a state of free oscillation the only forces acting on the e system are the internal — 
elastic forces and the inertia forces. If such a system consists of the masses 
m and the eR are yu), the inertia forces are 


| 


= y(t) 
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‘The displacements y(t) are > functions ns of the ti time and can b be e expressed as” 


_A mechanical eystem can ‘oncillate freely only with certain distinct fre- 
quencies Wie ; The amplitudes 4 Ys corresponding to each of these frequencies are 
: called the principal oscillations. — It can be shown that the amplitudes y; and 
corresponding to t different frequencies w and « and wi 


Eq. 4 is integrated o over all the masses the 
be: multiplied by an arbitrary factor c, because the product C Yi still — 
amplitudes of a possible state of ‘Principal oscillation this factor so 


symbol = is called the corresponding to and 
Eq. 4b is an expression for the ‘ ‘normalizing condition. - Different normal 
modes? and are orthogonal; that is, 


fe dm = 


_ Systems with an axis of symmetry | have two types of normal ‘modes—(1)_ 
1etrical modes for which ¢; has the same value in symmetrically ; situated 7 
“Points 9 and (2) asymmetrical modes for which the value of ¢; in ‘symmetrical — 


Points is equal but of opposite sign. It is usually found advantageous to 
_ consider and determine symmetrical and asymmetrical modes separately ye 7 
| Twoin important properties of the normal modes ¢; are: tO 
(a) If any mechanical system is subjected to to the inertia forces ws m gi as 


static load, the resulting deflections are just the normal mode os; and 7 
(b) The deflections u of any mechanical system loaded with aanekel arbi- 


os external forces P can be expressed as an expansion by all the normal 
vw 
Modes of the system— 


od: *“*Mechanical Vibrations,”’ by J. P. Den Hartog, McGraw-Hill Book Co., Inc., New York, N. Y., 7 
,1940,pp.192-196. 


7 
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own | jin which w is the circular frequency and ¢ is the time, in seconds. a 7 -_. 
quency, in cycles per minute, is 
ards 
here 
a 

| (4c) 
all 
a 

§g§g§$$@O 


| —provided that the forces 4 act only at points where masses m of the es 
chanical system. occur. The values of the coefficients c; are 


Fora approximate computations it is possible to consider o only a a few t 


Eq. 5 corresponding to the lowest frequencies of the system. :_ 
All the integrals | oid dm, ete., have been written in this form to include 


distributed masses m and concentrated masses M. For numerical -computa- 


tions such integrals 1 must be converted into an integral over the distributed | 


“masses plus a sum over all the concentrated ma masses; ‘that i is, Dinners 


Complementary Si Systems. —The .e following theoretical treatment involves two _ 


concepts: The first is the concept of systems, hich is 


Floor systems usually a a of beams ‘supported 
by main floor girders w hich h rest either on external supports like columns and 
walls o1 or, in complicated cases, on other - beams. Restricting the subject to 
statically determinate arrangements, it is required to show how the complicated — 
problem of oscillations of a floor system can be solved if the fre- 
quencies and the normal modes of a simpler system, having one main ye 
known. ‘If the simpler system is still too complex. for a direct solution, 
‘it can be reduced further r until, in the end, a system is derived which consists: 
only of ordinary beams, the frequencies and normal modes of which can be 


found by a number of well known methods. 


a One main girder i is selected w hich, while carrying ¢ other g girders of the sy ystem, 4 
rests in turn on external supports. In Fig. 1, for example, girder A should be ; 
selected. _ The girders normal to it would not be suitable because one a of 


each of girders rests on girder A. 


“a girder A to be I4 = oo. Thus, this not deflect and all girders” 


supported by i it can be antian resting on fixed supports. . This new system — 
(simplified ; system B), therefore, has one ne girder less than the original sys stem. 


it has normal modes ¢;— pand frequencies w, WipW hich, , temporarily, : are assumed to 


_ Apart from system ~ a system A is { formed which is identical with the — 


er z original system as far as dimensions and mass distribution are concerned. The — 
moment of inertia of girder A , Fig. 1, is Za. as in the original system. In all 


an 
4 however, the moment of i inertia is I[=o, A has cert: ain 
na 
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“or An important property of th - the two two artificial systems ‘Aand E Be ean n now be 


Theorem I.—Thé deflection u of a any , point in the original system, re- 
sulting from an arbitrary external load, is equal to the sum UA + up of the © 
deflections of the two caused d by the same load. 


if the 1 moments oments corresponding to the external load are called Mm, and | ~ 
em 


Zi 
Ds 
Fia. wall 


plated, are called m, { m, then the « deflection of the original | wi n consisting of of the 
gi , can be expressed as 


girders a, a, b, ¢, 
ich the three integrals are taken over girders a, b, and « eC, respectively. - 
uations can be written down for systems A and B. _ For system A, 


= 0, whereas, for or system B, Ta = Consequently, 


and, by comparison of E of Eqs. 8a and 9, 


Because of this systems A and B are named “ complementary systems.” 
e 


ae In the case of symmetrical systems it is necessary to modify the procedure 
slightly; otherwise the advantage of the symmetry would be lost.- Ina case 
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A will of inertia and I’, at their original the ] 
in all the other beams the value of J will be (©, . The important relation woul 
_ expressed i in Eq. 10 remains true. If the complementary systems are chosen 


7 in this way, the two ‘systems A A and B 
‘remain: ‘symmetrical and their normal 

nodes can be ‘separated in into symmetri- fF 

Por: 
cal and asymmetrical ones, which is of 


aymmetzical structures, system A con- rep 
sists. of only one flexible girder A which § use « 


1 supports the masses of the original sy om the 1 
by further girders” that are as- | met! 


J sumed 1 rigid. Fig. 3 shows the simplest I 
case possible, in which the flexible girder the : 
supports one rigid girder B. The am- shov 


-plitudes bs of then normal mode a are indicated i in this ‘diagram, . Ifthe amplitude | 
of the normal mode at point 1 is ga(l), all the am amplitudes along girder B can 


be expressed in form: 


— the system to be loaded with its inertia forces w 2b, dm, the reaction 


If sy stem A contained a concentrated mass M, at point - 1, the inertia force ie 
* 


F rom a comparison | of the two values of R, in Eqs. 12, is to spe 
cause the same inertia force as a concentrated mass: 


the heading, “ Facts Oscillations” (Section of g 


1), it has been stated that, when any mechanical system is loaded with ie = 


4 inertia forces w; m $i corresponding toa normal mode $i, the deflection of th 


+53 system is is equal i the normal mode gi Therefore, if the inertia force exerted a 

n by the fictional mass" M, is equal to the inertia force exerted by t the masses: a 

on girder B, it is ; obvious that the mechanical system A must have the same Qe 


normal modes and frequencies as a system consisting only of girder A if this E: 


girder carries the additional mass My, assumed concentrated at point | 
symbolic purposes the mass M; is termed the dynamic reaction” of girder B. 


dynamic reaction” of girder » 
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The value of the dynamic reaction is different from the value obtained by 
the law of static reactions; the latter, expressed in a form similar to Eq. 13, 


would have been. 


For a wana distributed ma mass m the dynamic reaction becomes m m=: 


In other more complex cases the procedure is similar. _ The system is loaded 
with its inertia forces and the reactions of the rigid girders on girder A are 


replaced by the inertia forces of fictional mannan @ or “dynamic reactions.” The 
use of this concept enables the designer to system to single 


| In ‘certain ee the > dynamic reactions become functions of 


irder § the amplitudes ¢ at other | points s of the girder. An. example of this condition is - 
shown in Fig. 4. In system A for this structure the members w ill have 
tude 

Can 


ments of inertia Ip = Ic = =Ip= Girder De carries a mass 
op in the center, carry distributed masses mg and me, re- 


~The amplitudes at points 1 and 2Qare $(1) and $(2); the dynamic 


| ‘Teactions due to the masses on girders B and C are found, as before, to obe — 


(13) and The amplitude at point 3 due to the 
tion of girder A is , which in an inertia force w? Mp 
the 


This force has reactions w? Mp ——__—— at both point 1 and point 2 To 
obtain the dynamic reactions, these values must be divided by w? ¢(1) od 


| 
| w* (2), respectively. ~The result gives total dynamic reactions at 1 and 


Mp o(1) + 6(2)_ -™ lp, Mo 


ypers October, 1049 1099 — 
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— 


‘This: result is not as simple as before, because M, and M, now now depend « on the 
ratio of the amplitudes at points 1 and 2. _ These ‘ ‘variable’ dynamic reactions 
can be u used without difficulty, » howe ever, to find the normal ‘modes, applying any 
the methods described ‘subsequently i in Section 3. 
Bor symmetrical structures system A may consist of two flexible girders, 
but dynamic re reactions can be e found in a similar ma manner if symmetrical and 

os asymmetrical modes are treated separately. = For example, to obtain the nor- 
- mal modes of the system shown i in Fig. 5 consider first the e symmetrical modes 
only. Thea amplitudes ¢ of both girders a a and a’ , Fig. 5, ‘must be alike—that i is 
$(1) = ¢(2). amplitudes for all points of girder b therefore, are 
(x) = = $(1). ‘The inertia loading of girder b is w? m = = and 


- Th 8 reaction is identical | with the inertia force of a mass s M, assumed acting 


point 1 and that of : a mass acting at point x 


Eqs. 16 are the expressions for the dy namic reactions for symmetrical modes 
which, i in this } case, could be derived from the law of static reactions. For 


asymmetrical “modes, (1) and the amplitudes of girder. b are 
= o(1) —— ioe The reactions of the inertia forces are 


ware 


and the dy namic on for asymmetrical modes become 


ly 


=, Eq. y ‘yields 


By using the concept of dynamic reactions any system A can be reduced — 
to a one-beam system, the normal modes of which can be found. a sore 


_ ‘The N Normal Modes s of t the Original System. —Thus far, the paper has shown 
how a system of 


f girders can be separated into two ‘idealised systems A A and B 


— 
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; called “complementary systems.” It remains to demonstrate that the fre- 
7 quencies ai and normal modes of the , original system can be built up fr from those 
15d) § of systems sAandB. The deflections of the original system due to a any external > 
an load are given by Eq. 10. _ Furthermore, the deflections | of any mechanical = 
the system can bes expressed i in terms 0 of it its no normal modes by Eq. 5. this 


Lions “equation to each of the two systems 1s A and B, 
and ug Sen as) 


and Eg. 18, substituted i in Eq. 10, yields 

~~ Eq. 19 defines the d deflections of ‘the original system | due to an arbitrary 


external load. The coefficients c can be from Eq. 6; thus, 
and 


- Any normal mode ¢ of the original system can be considered as the deflection 


the system due to the inertia forces w? m Applying 19 and 20 to 


this loading, Eq. 19 yields 

Pig dia + > cis Dis.. wee (21) 


Substitution of dz of aP = « = wd dm in m in Eq. 2 20 | gives, a 


The amplitude ¢ is now introduced from Eq. 21 except th that, to avoid confusion, 
a subscript k is substituted for in the summations: 


dm + | Deep dm )....(23a) 


pa dm + Les | bia den dm...... (24¢ 


dis dm... 


| 

4 

and 

16a) 

| Cia = dm and = cig = —— | din ddm.....(22) 

a a 
17a) 

iw 

17) 

4qs. 23 can be written in 

a 

— 


if 


simplicity, substitute 


Ba = | dia des 


_ and the equations for r the coefficients ¢ (Eqs. 24) become “ 


After changing the “dummy” indexes in Eq. 276 from k to j and from 7 tok, 


coefficients car can be expressed: 


4. k 
e coefficients CiA A ritten 


= = 


° 


¢ Such” equations a are capable of ‘solutions other than « Cia 
only if the determinant A of their coefficients \ vanishes; the condition A. a=0) 


is therefore the frequency equation. 


1102 FREQUENCY ANALYSIS— Papers | Oe 
i 
CB, 
(97 
If 
by 
“be 
4 and, by substitution in Eq. 27), 
ow 
— There are as many such equations 
fully, these equations become for be 
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ons: obtain 28, it was ‘necessary to divide Eq. 27b by 


‘This procedure i is only if that factor does not in other words, 


shows that w be if all the values Bis vanish, but this bia 


is trivial and of no o importance. The only | other possibility that « Wee May 


25b) become a frequency of the original system is when w;z is a multiple frequency ._ 
‘system B. Ifwxs = By the the possibility of this frequency beit being a root w 
@ = = W(k+1), By Kgs. 27 will be satisfied if all the coefficients a and 
(26) CB, except and C(k-+1), B vanish and provided that 


’ The number of equations detent to Eq. 30 i is ; equal to the number of functions 
gia used. If there is more than one mode Eq. 30 has solutions 


(27b) § which is the necessary condition for w = wep being a frequency o of the system. 
- if only ¢ one mode ¢za is used, Kq. 31 is automatically satisfied, and then every — 


to k, ~ double: root enn creates a root of the original system. : This root is not — 


by 286 and can only be deduced b by referring to Eqs. 27. 
The « question of multiple frequencies requires attention in practical “cases” 


(28a) because similar girders, or similar groups of girders, are sometimes used re- = 
~ ‘Peatedly i ina structure; and this fact leads to multiple frequencies ORB of system 7 
A case of frequencies Wep iS treated subsequently in Section 4. 
The Frequency Equation for Special Cases.—In general, | there will be as 
a ‘many ‘Eqs. 29 for the coefficients cia as there are e normal modes. pa. As ex- | 


plained in the next two paragraphs, the number of these modes to be considered © 

| Sy stem A is equivalent to a girder ca carrying a certain distribution of mass. — 
Ift there are n concentrated masses such a system has n different normal modes. a 
- In addition, if part of the mass | is distributed, the number of normal modes _ 


infinite. All these normal m nodes can be arranged i in order of their 
frequencies W1iA< < < ++ @nA. Fora system consisting of one girder 


will be. found that wea is considerably higher than w14; and much higher 
a than we, wea; ete. |For: an equally distributed load alone the | sequence \ of fre- 
quencies is = 1:4:9:-- 


If there were many coefficients, Cia, the computation would become > cumber- 
some and impracticable. Luckily, however, the expansions Eq. a, show 
marked convergence and the first few terms give excellent approximations. 


a he values of the coefficients decrease quickly due to the factors —— in 

=0 Eg. 22. For an equally distributed mass these factors decrease in the ratio 


T , and the. third one is negligible. 


pers 
| 
| | 
(288) 
ay 
| 
ritten 
a 
| 
| | 


rule presented subsequently under the heading, “Practical 
Bate of Modes ¢a and $3” (Section 1) can be used to decide how many 
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10des dia and das should ‘ed to obtain reliable  resulis. In normal 


= this rule requires es the use of f only on one or two modes $i . If, in addition — 
to a number of modes dp, only « one mode ‘gia is ; considered, the frequency 


equation, A = 0, becomes 


; Eq. 28a is used to compute the e amplitudes of the e normal mode corresponding 
to any one of the roots w w of this e 


or 


é By using Eqs. 25 and 26, ther normalizing condition can be transformed; thus, 


3 - Making further use of Eq. 27a the final form of the —_— owen is 


Eqs. 33 and 35 are sufficient to determine all the values C14. and cp. 


In addition to the roots of Kq. 32, any multiple frequency kB = W(k41))B 
a system Bisa a frequency of the original system, w= ‘ORB. 7 To obtain the — 
amplitudes of the mode corresponding to such a _ frequency, , Eqs. 27 must be 
used, as demonstrated subsequently i inSection4, 
32 to 35 enable the designer to com mpute the modes and frequencies 
of a system 1 if the modes and 1 frequencies of the simplified systems A and B are 7 
known. | System . A can always be reduced to a simple beam, the modes of which — 
3 can be: found by various 1 methods. — For ra system B that consists of a number 
of ordinary beams the modes. ip ‘and frequencies wip are the total of the modes 


‘ 
4 | : and frequencies of the separate parts of the system. 7 ~ If one of these beams has" 


a certain frequency Wi and a mode ¢;, system B as a whole has the he same fre- ‘ 

quency WOiB = Wie ~The amplitudes the corresponding for an any y point 

within the beam in question are ¢ig $i, whereas for any other point ¢iz = 0. 
"tt system B is not of such a je nature but contains beams supporting 
one another, a preliminary computation of the modes must be made by applying j 
the method L of this paper to system B. _ Through repeated application of this — 
“ method, one will reach a stage at which system B becomes a number of simple Z 


beams the modes of which can be found. oe 
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Discussion of the Frequency Equation. —The > frequency equation mn (Eq. 32) can 
be represented graphically in a v very simple and instructive manner. The 
abscissas in Fig. 6 are the frequencies w and the ordinates are the values of the 
left-hand ‘side of Eq. 32. For w = = 0, all terms of the sum vanish; if w (@ ap- 
proaches any one of the values wz, one term becomes + . The full curves 7 
in Fig. 6 represent the left ction Eq. 32, and the vertical dotted ed lines are. 


‘asymptotes to these curves. 
_ The dotted curves marked J, IJ, and 1 III in Fig ig. 6 indicate the value of the 7 


= side of Eq. 32 for three different values of 14 


"begins: with an an infinite value for w = 0, vanishes forw = w1a, and and then steadily 


i. - The natural frequencies of the system are represented by the intersection 
points of the full and of the dotted curves. Fig. 6 shows that there is always 


one frequency that is lower than any frequency of systems A and B; ; and ae 


is one frequency between each of the consecutive frequencies wiz, woz, 


: If there are n frequencies, wip, the total number of frequencies will be n + 1 bie 


Fig. 7 shows an example of this condition. Bis B 3 loaded 


a the special case | in which the complete system ccataine only masses" which | 


This total number can include a trivial root w =o, however, which cm " 


already form part of system Me on 


"Therefore, original system should have we one, or 
f 


three essen _ The system contains two masses only, however; and it has" 


_ two degrees of freedom because vertical displacements only are to be counted. 


— For mechanical reasons such a system has two natural frequencies, and no more. 


Th frequency formula (Eq. 32) for r this case, therefore, must give two proper 


The Flooring —Thus far only systems of beams supported by one another 
have been considered. _ The last element of a floor, however, i is the flooring, | 


Toots. ‘The: third root must be a trivial root w = «© ‘without physical meaning. 


which may be tl the open s steel type, checker plating, or reinforced concrete slabs, — 
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: and it is necessary to show | how to introduce the flooring into tl the computations. 
Because of the ‘small span between the secondary beams, the flooring v will not 
‘deflect appreciably from beam to beam, but it must move with the supporting | 


beams and its mass may be a considerable e part of the: total n mass sof the st aoer pel 


tion is justified, it = to pipe a more precise e computation. . If the 
flooring is not continuous over its supporting beams, it can be considered as a 
number of narrow beams spanning the secondary g girders. This assumption 
reduces the system to a | type of structure that has already been treated. eed 


relatively rigid and have high natural frequencies. is justifiable, therefore, 
regard these beams as entirely rigid—that i is, I 


_ Without presenting details, the result of computations, using finite difference 

equations, for a structure consisting iz of n equal secondary beams 

ort, Ta, wit! with equal ‘spacing and loading, is as follows: 7 The values of the f fre- 


inci mp is the mass of the beams per unit length; ‘m, is the mass of the 


0 flooring per unit area; bi is the span; J is the moment of inertia; and a is is the 

Spacing of the beams. The quantity cos cos can vary only 


and it and, therefore, all the frequencies must lie b betw een t the mae 


To compare Eq. 36 with the frequency obtained if the m mass of the flooring 
te concentrated on the beams, mg = mg + am, and m, = 0 are introduced in 
mp + ams 
‘The a approximate value w derived from Eq. 38 is equal to the. lower limit for 
- the exact frequencies w;. Ifn > 2 or if Ms is small compared to mg, there will 
: be little difference between the fundamental frequency @1 as duthend from Eq. 


36 and the approximate value of w derived from Eq. 38. . a 


‘Therefore, the concentration of the mass of the flooring on the beams is a 
good approximation . for the fundamental frequency if there are more than two — 


a ™, is small compared with ms; but, if higher frequencies are required, a more 
exact computation must be used. * This conclusion remains s true if the secondary 


- beams beams carry unequal masses or are otherwise different. — 


secondary beams. — If there are two secondary beams, it is still justifiable 5 | 
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If the secondary beams are supported by o other beams forming part o of a 
more complicated | structure, the approximation | can still be used as long as 3 the 
frequencies which are ultimately determined remain below the lowest frequency 
@ of any of the secondary beams alone. — On the other hand, if reliable fre- 


quences higher than the lowent are required, a more exact computa- 


eases, however, it imply that the structure has a considerable 
emer of resonant frequencies within the range of the running frequency which 
_ would make the: structure unsuitable for its purpose. The difference between 
i the approximate e and the more exact solution is sh is shown own subsequently i in Section 4 
Computations— 


Mass Distribution. —To obtain correct results it is of the utmost importance 
that the mass distribution used for the ‘computation ‘shall represent the actual 
distribution under working conditions. It would be quite erroneous to include _ 
the mass of the accidental live load, as as specified | for t the design, in in the customary 
manner. In most cases only the ‘permanent load needs to ee considered; but, 

‘if considerable live loads, such as those in water tanks | toa coal bunkers, can 


a Another consideration is the matter of accounting for the mass of machines. _ 


_A machine can usually be regarded as a rigid body supported by the structure 
at three or more points. _ Structural engineers are likely to assume that the 
total mass of such a machine distributes itself over the points « of support in 
the same ratio as the static reactions. This assumption is rarely justified. 
A correct procedure for replacing a machine ne by its ‘ ‘dynamic re: reactions” at the 


points of support is treated in Section 3. 7 aa apa 
_ Number o of Modes GA and ¢z.—As stated, it is not necessary to use all the 


frequencies and normal modes of | the complementary systems A and B, but 
the higher modes can be neglected. — _ The decision as to which modes und be 


retained depends on the value o of the frequencies to be investigated. To be 

‘sure that a certain: frequency w, Wi, found by these computations, is reasonably — 

exact, all the frequencies" up to 2 should be used. It is quite obvious, 
= furthermore, that at least one mode of each system m must be considered. If 


system B consists of a number of separate girders, one mode of each of these 


oe To check the suitability of a a structure for a certain frequency w, ‘w, only the 
‘Tesonant up to about 15, are Tequired. for the 
fre 

7 this rule requires the use of only one mode of system A, and of wile one ene 
of each of the girders forming system B. 


Factor of Safety Against Resonance. —After finding the natural frequencies 
ofa system they must be compared with the frequency of of the machines carried. 

it if the machines run at variable speed, a an entire range of frequencies n must be 
garded as dangerous. : To obtain the smallest possible vibratory | deflections — 


and: to avoid resonance , during the starting o of the e machines, the running | fre- 
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quer one hei sufficiently lower (40% being a safe margin) than the funda- 
mental frequency of the the structure. 
If this condition cannot ; reasonably be fulfilled, either because the running = 
speed i is high or because the span of the structure is large, the system will still e 
prove satisfactory if the running speed is kept clear “of the nearest resonant _ 
frequencies. It i is essential, _ howev er, that as few frequencies as possible be 
below t the running speed. 
i: _ Concrete and Reinforced Concrete.—If the structure is concrete or if steel 
beams incased in concrete ; are used, an uncertainty is introduced because the 
value of E for concrete varies considerably. The only safe procedure is to 
determine the frequencies for limiting values o f E and to regard the entire 
range between these frequencies as dangerous. — When hen fixing these limits, it 
should be kept i in mind that vibrations of structures must be restricted to onal 
amplitudes and , therefore, to low vibratory stresses. — ~The value of E for 
such | low stresses is likely to be ; to be higher than the velo of E used for stress 


computations, 


Subsequent Alteration of Sections. the frequencies of a system are 
found to be unsatisfactory, the moments of inertia of some of the members are 
eae and the determination of the frequencies is repeated. : It is not nec — 


& to repeat at the entire = computation if, the e following property of the modes 
isused: 


‘Theorem II.—If all the moments of inertia, m4 of the beams of a system 
are changed to a I and if the ratio @ is the same for all beams and the mass" 
distribution remains the same, the frequencies Wi, We, ***, Of the original 
system are changed to ¥a@w1, ¥ @ws, The numerical values of the 
amplitudes of the normal modes $1, $2, «++, corresponding to these fre- 
quencies remain entirely unaffected by sucha change. | 


_ Theorem II ma may be applied to any system. — It can be useful if changes are 
- unde 1 in the moments ts of inertia which affect only the beams of one of the 
complementary ‘systems. By the terms of this theorem the modes of this 
4 complementary system cones used ag again, and only the frequencies re remain to 
beadjusted, 
Systems That Include Statically Indeterminate Beams. —At the beginning 

of Section 1 the discussion was restricted to 0 statically determinate arrangements — ; 
of floors; but the method developed can ke. applied without modification to 
floors that include le statically indeterminate members as long as the arrangement 


of. the beams i: is such that the re ‘reactions from each beam on the o others 


__ Consider the floor shown i in plan i in Fig. s. 1 and assume beam A to be part 
ofa portal frame consisting of beam A and its supporting columns. | ‘The struc- 

ane as a whole is statically indeterminate, but the reactions from the secondary 


beams on beam A remain statically determinate. 
This type e of structure can be dealt with by the method of ‘complementary 
a _ systems” because the fundamental Eq. 10 applies. By’ this method it is pos- 
; sible to reduce the floor system | gradually to one of the basic ¢ types of statically 


indeterminate structures, like continuous beams or portal frames. Such struc- 
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tures can be dealt with by k known methods—as, for example, those explained 


inSection 3, 


-Statically indeterminate of other types require different treatment and 


are considered in 2. 


2. THEory APPLIED TO (STRUCTURES 


The procedures explained in Section 1 apply only to structures that are 
statically determinate or that: belong to a special type of indeterminate e systems 
mentioned at the very end of Section LL - Frequently beams are continuous 
because their end connections are rigid and able to transfer bending moments. 
To estimate the influence of rigid connections and to treat floors that ory 
continuous beams, or or are otherwise indeterminate, the method developed in 


Section 1 is extended, as follows: 


First , imagine | that all the continuous beams are cut into simple beams, 


leaving a floor system of statically determinate elements. Frequencies w; and _ 
normal modes i of such a system can be found by the procedure i in Section | 1. 
—_ Next (as i in the theory of statically indeterminate structures), assume that 
moments X; act on each of f the points k where the continuous beams were cut. 
- Each of the moments Xi = 1 acting alone on the statically determinate basic 


struc es certain deflections, yz. . If a system is symmetrical, it is 
to treat t symmetrical and asymmetrical modes 


t selected as statically indeterminate values Ww hereas for 
X,= —X,=1isa suitable indeterminate value. = 


Since the deflections u of any statically indeterminate structure are equal — 


the deflections Ua of the statically determinate basic structure due to the 


‘same load, plus the deflections ux due to the statically indeterminate moments 
the former can be expressed by the normal modes of the basic structure 


= Le di (see Eq. 5) and the latter a are ux = UX. Vee | The amplitude o of 
te normal mode, ¢, of t f the statically i indeterminate system can be found as the 


& “deflections u ¢ of the system due to the inertia a loading « wo dm, from which _ 


= Leg + 


The coefficients c; in this expansion are defined by Eq. 6 § which (after introduc- _ 


ing dP = w? #m) yields 


¢o;dm = $j dm + | $; dm... 


‘Using the orthogonality relations (Eq. 46 and Eq. 4c) the first sum on the right- 
hand side can be reduced to C; Then, replacing the dummy index ; i by 2 iand- 


introducing: the he abbreviation: 
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the form i is 


Additional sions be obtained from the theory of statically indeter- 
-minate structures. If the moments due to any loading Pona , statically inde- 
terminate structure are M p and if the moments created by the action of any 
statically indeterminate | moment X; on 


structure are m;, , then the 


0 


L.4 $ ie, the moment M, can be e expressed by ee second derivative of the — 


i= AL "Therefore, 

| 


If these values are substituted for. a, and mj; in a Ba, seal 


Jee y";de = 


fies 


he 


_ The integral in the first term of Eq. 46 is the well known expression for the 
virtual work done by tl the loading w; i Pit dm, which creates the deflection pis 
if the system undergoes. a deflection Yj. . According to Betti’s law 
a work must be > equal to the virtual work done by the moment x; 
system, undergoes a deflection 
The virtual work done by the loads w; dm on the deflection is obviously 
fe Vi dm; and, with reference to Eq. = Bij The 


in 1 which Vii is the : a between the two tangents to the curve te at point j. 


> 


far 
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‘integrals i in the second term of 46 can be interpreted i ina ‘similar 7 

They represent the angle between the two tangents to at 


* 
“There are as many such equations as there are static: ally indeterminate 
values of X;. Eqs. 42 and 49 form. set of linear and homogeneous equations 
for the values of and Xi. Such equations have : solutions ~Oand 


s. The condition 


‘obtain normalized solutions the Ci onl X; must be deter- 
mined in such a way that Eq. 


4b is satisfied. The normalizing condition can 
be brought into the form: 


o apply the formulas derived in Section 2 it is first necessary to find the 


values of the coefficients defined by Eq. There would be “no 


mental ‘difficulty in computing the value of the integral in Eq. but 
Eq. 47 gives an alternative relation 


which avoids the determination of 
any integral. From the -right- 


Vii 


in which ¥;; is the angle between the two tangents of he at t point j j. 


is Positive if if xX; does positive wi work k during its formation. _ 


: 3. FORMULAS FOR FOR REFERENCE 


fe ‘ormulas s for Frequencies ond M shee of Simple Beams.— 


—A 
uted mass m a i, as shown in Fig. 8, has the napenien 7 


42 2 


t, I, in inches‘, w, i in » pounds per foot, and 


The value a,; can easily be computed by differentiation of ‘Substitution 
of Eqs. 47 and 48 into Eq. 46 yields 
| 
4 This angle = 
— | 
a | in which is measured in fee 
= 30,000,000 Ib per in.2.. Th 
= — sin — (i = 1, 2,3, (58) © 
| 
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jue of the following integral i is frequently required: 


—The masses are at point 1, Ms at point 2, ete., a as ‘shown i in Fig. 9. 
such a § a sy ystem n vibrates with the circular frequency ¢ w and the amplitudes are 


Yay the inertia forces acting at points 1, 2, --- are Mi yi, M2 etc. 
These inertia forces cause the deflections of the otis: - and, if the deflections 


due to these forces are determined, the result must be ana to the original 


an amplitudes Yi, ete. This relationship can can be used used to find t the natural 
quencies s and modes of such a system. J a 


to a av unit load at point are ‘the deflection mys due to the inertia 


ean be expressed i in the form: 
— 


= = w? Yi bu w? Y2 bo + M3 ys + (5 5) 


a If there are n “masses, there are nm such equations for the n nm values of 4 Yr 
are and have only trivial solutions v= 0 unless the 


rat Therefore, A = 0 is the  fre- 
quency equation: 


7) 


Eqs. 57 ar are of the 1 nth order i in te terms of 1/w? and have n roots which need 
q “not all be different. 7 It can be shown that all these roots 1/w? are positive a and 


q real. Each a be introduced in Eqs. 56 56 yielding a set of amplitudes yi 
for each frequency w. 


%**Mechanical Vibrations,”’ 


+ J. P. Den Hartog, McGraw-Hill Book Co., Inc., New York, 
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- 
Every set of "values S Yi can . be multiplied by an arbitrary factor. If _ 
~ factor i is suitably chosen, the normalizing condition i is satisfied and the ne ampli-— 


tudes y; become the normal mode The condition is 


This method is v very — e for sy ais a small number ¢ of masses. - 


if the number is high, the computations become rather involved, but approxi- 
mate solutions can be derived by the matrix method.‘ ' The latter is not 


- explained i in this | paper because such cases can be solved by the Stodola-Vianello 
method, which is based on the same principle. The Stodola-Vianello a 


is the more advantageous because it is more akin to engineers than the rather 


The Method 1 of Stod Stodola- Vianello—1 —TI The principle of this method i is as follows: 


To obtain. a series of approximations of the _ fundamental mode assume & 
first approximation y™ by pure guessing. ‘Then, compute the of 


the system due to the loading m m y®, if y® were were the exact mode ¢, o a mul- 


tiple c ¢ of it, the resulting deflection y® would be equal to —y, because the 


- deflections of any system due to the inertia loading w? m ¢ corresponding to a 
normal mode are equal to this same Since the amplitude is 
‘somewhat like $, the amplitude y 2 y® will still be roughly similar in shape to that 

y® but it will be on a different scale because of the factor 1/w*. 

| Tt can be proved that y? is nearer to the correct shape of the > principal 7 

~ oscillation than is y™. if the | process is repeated, a series of approximations 4 


yy, are obtained. . After a few steps t these curves do 1 not differ in 
shape any more—they differ only in scale and the ratios y®/y®, y y™/y, 
are approximate © wae of w. The most exact value i is the average value of 


Any me y method, graphical or numerical, ca n be us used to find t the consecutive 
deflections y™, y®, yO, --- . The foregoing procedure can be applied to 


determinate as well as indeterminate structures. 


a One suitable method of finding the deflections of beams with 1 either constant. 7 
or variable moment of inertia J is to divide the beams in a number of short 

segments a (see Fig. 10) for which the moments of inertia can be “wes to be Po 
constant, ‘and to concentrate the masses in the division px points 
Then, the moments a ‘ at the division points due to the inertia loading -— 
oe According to the lei of “elastic weights” the deflections of a girder — 
are equal to the bending moments in a beam of equal span created by the 


_ 


e __ 4**Mathematical Methods in Engineering,” by T. von K4rm4n and M. A. Biot, McGraw-Hill Book Co. oe 
Ine., New York, N. Y., 1940, pp. 196-204. 
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‘elastic weights” W1, ws W3, t res 


are given by the formula: 


SFT, 


2m) +t SET 2 + M 


4 If th the intervals ls a; and the moments of inertia J; are wmatonh, ‘Eq. 60a 
3 


Wi = 


Since the masses M; are now concentrated at the fivision p points, t the in inte- 


in Eq. 59 into 


M; 


~The amplitudes c of the principal mode y y are multiplied by a factor cand 
7 become the normal mode = =cy. _ This factor c is found from the normalizing 


condition > 


Ifa system is symmetrical, the e method « can be used to! find the 1¢ lowest sym-_ 
metrical and the lowest ‘asymmetrical modes. These two modes are obtained 
by —* once with a th a symmetrical c curve y®, and once with an asymmetrical 


A “Generalization of the Stodola-Vianello Method. —The Stodola-Vianello 
method can bet used to find the fundamental frequency of simple beams but it 
‘is possible to extend it further and to o determine higher frequencies. Making 
use of the fact that “een mode go and any higher mode ¢ are 
-orthogona 


the designer can can find the second mode by the same me procedure, if he begins 3 with 
an assumed first approximation y™ which is orthogonal to 1; thus, | 


dm = 


— 
Functions y® which 


formula: 


5 “‘Statically Indeterminate Stresses,” by J. I. Parcel and G. Maney, John Wiley & Sons, Inc., New 
‘York, N. Y., 2d Ed., 1986, p.47, 


— 


, these weights in w 
— 
orth 
. _ If this method of sed, the entire computation can be nF I 
— | 
—m 
orth 
but 
.! 
mas 
cept 
4 is st 
— 
4 
mac 
‘poir 
and 
— 


This relation can easily be « checked by : 
= substituting it in Eq. 63d. 
=. Beginning with a y® a according to Eq. | 64 further : approximations 
can be obtained | by following the same procedure as for the first mode. If the | 


first mode i is known exactly, it \ would be sufficient ti to start with a a function y® 


which is s orthogonal to 1 and every ‘subsequent approximation would remain _ 
’ orthogonal to ¢:; but ¢; is known approximately only, and, therefore, Kq. 640 
' gives a function y® which i is approximately orthogonal i in relation to the exact — _ 
; amplitudes: ¢:. Therefore, it is necessary to insure that 1 yr remains s orthogonal - 

to by applying Eq. 64 to the subsequent approximations y®, y®, 
— ‘If further modes are e required the same process can be used, but the normal- 


} 
izing condition— 
| 


red. This procedure insures that the approximations are 
_ orthogonal to all the low er modes which a are already known. The convergence | 
but the first four or five frequencies can. be computed with reasonable ac neounney. 
. Replacing Rigid Bodies by Their Dynamic Reactions. -—The masses supported 
by a | beam frequently comprise machines that can be: considered as rigid | bodies. 
| “Fore computations it is desirable to replace such bodies by ‘equivalent i imaginary 


masses at the of This can be using the con- 


own in 110). The is four 


‘is indeterminate because a1 rigid an only y three sur supports. 
_ There is an important case in which this difficulty disappears. If axis y . 
is an axis of symmetry of t the machine and if the amplitudes | of girders T; ot 


7: are are equal, or - approximately ‘equal, each girder ca carries one half of the load. 


machine (see Fig. 11(a)); Ji is ‘the polar moment of the : mass | M in relation rag 
point O, which is the projection of the center of gravity on the girder axis. ' If 
_ the beam oscillates with amplitudes y, and y2 as shown in (Fig. 1 1(a)), the move- 
ment o of the machine consists of a vertical movement and of a rotation 1 around 
point O. Therefore, the inertia forces are a vertical force M 


and a couple w? yuo 8 . Ther reactions at point 1 and point 2 om 


Imaginary masses at points 1 and 2 would create inertia forces: 
4 1 = and 2 = w? Me Y2... 
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of T, 


- beam; 10.1n by 43-in @ 251d. 


a I - beam; 20-in by 


6}.in @ 65 Ib 


| 
—— 

— 

| 
— 


ao 
can be e replaced by the i masses: 
mu, >M+J/) , — 
+ + ly)? yo 


These fictional m masses M, and Mz [, still depend or on the ratio ; yilys but they can, 
"nevertheless, be used for frequency computations without any difficulty (see 7 
- Section 4). Eqs. 68 apply for the case of equal amplitudes y: = ys, in the form: _ 


Eq. 69 indicates a mass di stribution according to the law of statical reactions. 


, It is is s correct to use veel 69 ) only if he amplitudes of the points of ‘support: are 


| 
| 


a If the machine is not symmetrical or if the girders T, and Ts do not deflect. - 
‘equally, it is necessary w remove one of the redundant supports, analyze the 


‘this s support 4 by the method aaa’ in Section 2 for statically ‘indeterminate 


— To demonstrate the principle, consider again the machine shown in Fig. 


1 (b). ye ‘The machine still has the axis of. symmetry y—y y but t the deflections 
of T, and are assumed to be different. 


~The total mass of the machine i is 2 M and there w will ll be ty two polar moments 
J, and Jy which refer to axes parallel to and y, respectively. In > 
order not to lose the advantage of the symmetry but still to obtain a statically | 7 
determinate case, remove the supports at points 2 and 4 and support the 
"machine instead at point 5 halfway between points 2 and 4 by an imaginary — 
‘rigid girder that rests on girders T, and Te at points | 2 and 4 (see Fig. 1l(c)). 
Expressions for the fictional masses My, M2, Ms, and M, as functions of 
] the amplitudes y:, y2, ys, and ys can be derived in a manner similar to that used 
for Eqs. 68. The frequencies and modes of the e system carrying th these imaginary ; 


masses can then be found by the method given in Section 1. - The modes’ _ 
this statically determinate basic system can then be explained in 

Section 2, to find the e frequencies and modes of the actual statically indetermi- 
nate system. _ The statically indeterminate values X defined in Section 2, in _ 


this case, are a . couple of vertical forces acting at points 2 and 4 on the machine, 


and an opposite couple of equal size acting at points 2 and 4 on the girders T; 
and Ts. ‘The value of these couples: will be such that the imaginary gap be- 


tween the machine : and the s supporting girders i is ‘closed. ee ight 7 


ot should be emphasized that, to obtain fairly exact results, it is really 


‘Necessary to use the procedure i indicated i in| this subsection. . Toc distribute the 2 
mass M over the supports simply following the laws of static reactions and to 


disregard the restraint due to the indeterminate four-point support of machines 
would lead to considerable « errors in ” rrors in the frequencies (particularly in the higher 
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- frequencies) except when the machines represent only a small part of the total — 


The dynamic reactions thus found are used to replace the actual masses * 


rigid bodies all through the computations. To determine ‘the value of an 


expression such as | ¢? dm, the amplitudes y; in Eqs. 68 must be interpreted 


to mean the values of ¢ at points i. In the case om an ‘integral “ba os dm, 


“the dynamic r reactions can be computed by introducing ‘the amplitudes of either 
oa or into these > equations; but da, or én, must be used consistently for all 
pain reactions. — The final value of the integral does not depend on which | 


of the two is is used. 


ILLUSTRATIVE EXAMPLE 


- To demonstrate the application of the theory | developed in slate’ 1, con- 


sider’ a rather which, nevertheless, will serve as a guide to 
procedure in more cases. The structure is shown in Fig. 12. 


to Ts, in turn carry ‘a steel | flooring 13 17 in. deep weighing 14 Ib per 


‘7 The secondary beams Ts, 7. Ts, a and T; ¢: carry two identical machines, 

beams shown? in Fic. 12 are British standard sections, and would be suitable 
to carry the dead load, the weight of the machines plus 50% impact, and a dis- 


tributed live load of 100 lb per sq { ft. Iti is desired to investigate whether this 

structure is safe against r resonance if the machines run at 320 rpm. 
— Only the masses that are on the structure under working conditions s should 


be taken into account. _ These are the masses s of the flooring, of all the beams, 


a appreciable live load occurs eal no allowance will be ‘made for it. For the 
P ‘computation of the lowest frequency of a system 1 the distributed mass of the 
_ flooring can be assumed concentrated on the secondary floor beams (see » Sec 
1). Making use of this consists of only the 

_— The system has two axes of symmetry a —a and. b - as indicated in Fig. 
In general, computations can be simplified by considering symmetrical 
“and asymmetrical modes separately. ? In this particular case it is advantageous 
~~ use the ; axis of symmetry a—a and the modes which are symmetrical in 

* relation to a— —a are considered separately from those which are asymmetrical 
_ inr relation 1 to < a—a. For structures that do not contain continuous beams the 

lowest frequency corresponds, as a rule, toa symmetrical mode a and, therefore, 

1 the designer needs to compute symmetrical modes only. In computations of 
this nature it is essential to keep the dimensions uniform. _ Throughout this 

- example : all distances are in inches, weights ar are in 1 pounds, and time is in 1 seconds 
— definitely stated otherwise. — The acceleration of the earth in these 

imensions is g = 386 in. per sec.2. 
nee Choice of Complementary Systems.—As explained in Section 1, the system is 
- divided into two complementary | systems A and B. © Girder Ty (Fig. 12) i is the 
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only one that rests entirely on external supports. — To form system B, therefore, 
assume girder Ty to have Is whe 00 , , all other girders retaining their actual 


4 moments of inertia. . System A is obtained by leaving the moment ; of | inertia 
n | of girder on at its original value, the value of I to be © for all other girders. 
Referring to system B, the assumption = transforms all the girders 
d | to Ts into simple beams on fixed supports. s. The n normal modes of system BL 
‘therefore, comprise all the modes of pat 
ny the girders T, to Ts. These modes — 
_ § can be found by the methods i in Sec- 7 t a 
tion 3. ‘The choice of complemen-| “Fre. 
ll tary systems, in this case, 
tary sys ems, in this case, has not 
: been influenced bya any - consideration of symmetry; these systems can be u used 
- for both types of modes, symmetrical and asymmetrical. a 
_ The Normal Modes of System B System B consists of the eight secondary 
n- beams 3 which (as a a result of the simplification ¢ concentrating the mass of the 7 
to ‘flooring on the beams) “oscillate independently from each other. The four 
2. — T,, Ts, Ts, and T's (see Fig. 13) are mechanically alike and carry an equal — 
ms ‘mass. Their frequencies” and modes must The equally dis- 
er tributed load wis5 X + 21 = 91 Ib per ft. The total m mass per girder is 
ce. ml =— = = 3.7 The moment of inertia is I = 81. 
2d inf Fort = 1, Eq. 52 gives the fundamental frequency of these beams as: wip 
his = = wep 94.0. The amplitudes of of the modes 
wld] tesiven by Eq. 53: = dun = = 
From Kq. 54, dip dx bap dt = = 
Jec- 0.874. The dynamic reactions of each of these 
beams ¢ on Ty are: M, = M; = Mg 1.26. 
Fig. 
decide if if higher modes of these beams | are required the rule given in 
vo ‘T} Section 1—that only frequencies below 3 @ need to be considered— —is applied. 
in The machines run at 320 rpm which yields w 33.5; and, conse- 
quently, _modes corresponding to frequencies ‘greater than wo = 
a neglected. i _ As the fundamental frequency of these beams is already 94, no 
- of higher modes are required, and the designer can proceed to the more petal nae 
this. task of finding the modes 3 of th the g > girders Ts, T: 3) Ts, and Tr 7 ‘These latter girders — 
onds | £ are again equal in dimensions and mass and must have identical frequencies 
em 
s the 
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by concentrated masses; ar and the present case is one to which 6 68 


As shown i in Fig. 14 each of f these girders carries distributed load and 


one half of a | machine. - Then moment of inertia of these girders is I= 122. 3. 
=5xX 14425 = 95 Ib per ft, the 


m= = 0.246 ———. The we eight of each machine is 6,720 lb, of 


4 


which one half is arried by di ‘The mass, therefore, is M = 


2x 2X 386 


8.70 ——._ To replace the machine by its dynamic reactions at the } points 


of support A and B, the polar moment J of the machine is required. If this _ 


— 
= -_ . 


per 


is equal to 336 330. 6.96 these in 


Eqs. 68 the dynamic. reactions are 
= BX870 46.96 2 x 2X 8.70. 6. YB _ 


70 + 6.96 


_ Application of the Stodola-Vianello Method to the Determination of the Mode 
of Te -—As shown in ‘Fig. 14(6) the span is divided | into eight equal parts and 


fi the masses are concentrated at the division points. The choice of the division 
points has been made in such a way that importan t and large concentrated 
“masses, such as the dynamic Ma Mz coincide with ‘division 


are: 


? points. The masses at these points { 


pound- 


inches 
ails = Ms= = Me = M; = 2 X 0.246 = 


| and 
| 

— 
with 
ag The 
48 28 28 — | Ge 

| 

+ 1.7442. (70a) 
= 2.61 + 1.74 ==. . (706) Cols. 
— 
an 
‘that 
a 


Me 0.246 + 2.61 + 1.74% 


M,=2x0. 246 +261 - + 1.744 = 


+ 1740" 3.10 4 1.7444 


Points A and B in Fig. 14(a) are points 2 and 4 in Fig. 14(b), therefore, ya = ye 
yp = = As explained in Section 3 the Stodola-Vianello method makes 
it possible to obtain the normal modes by successive approximation. — The 7 


first approximation ny is found by guessing x only; then, the beam is loaded | 
with the inertia forces corresponding to saad and the deflections are determined. 


These deflections are the next approximation mn y® of the mode; the process is 


then repeated as often as required. - 


TABLE 1.— SSTopoLa-VIANELLO Meruop, First Step 


mx! [moe] | 
ty}? en EI |. 


ries (6) 


5.908 
010 | 0.22 
4.15 11. . 23, 012 


13.404 | 13.3: 83.236 | 3,396 


542] 1382] 88. 34,145 
3,100 


0.49 0.12 


0. 03 0. 12 


~ 


Col, 2, 3, and 4, Table 1, pees the values of the assumed sctetiieaion: > 
‘the masses M, and M TP. . To obtain the numerical values of the masses M 


at points 2 and 4, it must be carefully remembered that the ratios Ye on and — cE in n 


"the expressions for these masses refer to the amplitudes just considered and 
‘that, therefore, y y2 = 0.80 and y, Y4 1 = 0.95. (see Col. 2 , Table 1) Ww hich yields 


M= 3.10 + 1.74 a 5.19; and M, = 3.10 + 1. 74005 = 4.57, 

Cols. 5, 6, and 7 7 are headed “Tnertia Loading” and the load 


My®, the shear forces | 
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7 
Also_ 
b) 
| 
Point | M (RT 
a | @ | | | | ao | ay | aa 
| 045 | 0.49 0.40 
2 0.80 | 5.19 0.731 
3 | 1.00 | 0.49 093 ea 
4 | 0.95 | 0.994 
er 
In 0.50 | 0.25 | 7.768 | 7.73 | 59.472 | 1,733 | 0.67 
7 | 0.25 3.946 | 3.92 31.694 | 492 | 0.36 
on 
ed 
> 
on 
we a 
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“the reactions at point 0 are found, the shear forces are dete determined, one by one, 


j 
_ by sah subtracting the point load from the shear in the preceding panel. Since 


‘the panel length is is constant, a = 24 in ., the v values ¢ of the moment — can be 


found, step by step, by a wee the shear force to the the moment in in . the preceding 


panel point, with M M = for 0. 


Ws the shear, due to the loads w,— 
all w with suitable ‘multipliers. — Col. 10 is the approximation aad of the funda- 


Col. 11 gives the values of of” M Cy} which are ‘required to find the f fre- 


quency from Eq. 61: wt = ;e= X —> = 0.01126 


66,493 X a® 493 xX a®’ 66,493 493 


30 X 108 x 122. 


2, 990; w = 54.7. obtain the values in Col. Ul, 
the masses at points 2 and 4 “must be recomputed, using the values of ye 


= 64.864 and ys = 88.288; thus, Mz = 3.1041 288 5.47 and 


810-4 174 = 4.38. 


Col. 12, ‘Table 1, finally, contains the values: 


| 


If al had k been the correct mode, w? y 1 would have been exactly equal to toy, 
and beans actual differences en enable : the designer t to judge the degree of approxi- 
mation. As this is s only the fir first st step, there are still noticeable differences. 
a The com computation i is now - repeated, beginning with the values y just found. 
A piimaiged mode can always be multiplied by a factor; it is therefore permissible 


_ to start the next step with the values w? y. This procedure i is conv enient 

-heomuse otherwise the numbers i in the tables would i increase after every step. ‘ 
It is not n necessary to reproduce the co computation for the further steps s since the Doms 
tables are quite similar to ) Table 1. The results of the first three steps | ays 

shown in Cols. 2, 3, and Table 2 

‘The result w? ty? of the second s step (Col. 1. 3) differs very little from wy 

ea 


al 2) ‘and the third step (Col. 4) produces no change. ‘The result of the 
_ third step, therefore, is selected as final. — _ The values in Table 2 represent t the 


_ 


COND 


hud 
a Coho 


‘emplitades of a principal mode and therefore contain an arbitrary factor. To | 7 
find the normal mode this f factor must be chosen to — ened = T 


the: 
4 

| 

a 
sequ 
— and 

— f 
1 
346 EFI... 1 

= 
— 

| 


¥ ma quantity = 83 39 9 has been computed in Ce 
condition will be satisfied by 


= 0.345 yo 
ail 


Col: 9, Table 2, contains the v value of 28 which will be required 


2B = W3B = WeB = W7B => 54. 6; 


mT = 
The exact value of the next frequency ¢ of these | beams i is not t know n without 
further computation, but it will be at least four times the fundamental fre- ‘ 
- queney—that i is, more than 200. As this value is greater than the previously 
F explained limit of # = 100, further modes of these beams need not be considered. 


TABLE TODOLA-V IANELLO METHOD, RESULT 0 OF 


Dynamic Reaction Ra 


Second 


M | My oon] M 
alo 4) | (s) | ©) | s) | (9) | qo) | on 


0.25 
0. 0.375 5 


Non 


43 
0.07 
0.03 


’ 0.01 


4.75 
54.6 | 


_ The dynamic reactions of these girders on the main girder T, are given by 


‘the integral formula Eq. 13. For concentrated masses M this equation yields 
— 


| “with Fig. 3 the dimension z must be measured | Pa the far vad: of the girder— — 


that i is, in this case from point 8, as shown in Pe. “). FW hen determining i 
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l. 6, Table 2, and 
| 
bec 
6 
7 
d 
2) 
— 
1) 
‘1 040 0.396 | 0.396 | 0.49 | 0.08 | 0.137 0.059 
cls 2 | 0.731 | 0.727 | 0.727 | 5.48 | 2.90 | 0.251 1.032 4.26 a roe 
3 | 0.93 | 0.934 0.934 0.49 | 0.43 | 0.322 0.099 0.49 7 cae 
; 40 t 0.994 | 0.993 | 0.993 | 4.37 | 4.31 | 0.343 0 | 0.750 5.71 re 
d 5 0.89 | 0.896 | 0.896 | 0.49 | 0.39 | 0.309 | 75 | 0.057 0.49 . - ee 
ai 6 | 067 | 0.672 | 0.672 | 0.49 | 0.22 | 0.239 90 | 0.029 046 a aa 
0.358 | 0.358 | 0.49 0.06 0.124 | 25 | 0.008 | 0 


the masses Mi in n Col. 11, Table 2, for points 2 and 4, the amplitudes ; y at these of ‘this 


points are required. (Amplitudes 4 y always refer to the Numb 

aspen which are, in this case, Ye = = 0.75 Yo ar and y= 0. -50 Yo.) symme 

Finally, the e dynamic reactions of these girders are: R= = = metry, 

| the tw 
= 4.75 

7 ‘The Normal 1 Modes « of System A— —System A consists of girder Ts, which has | x 10- 


a moment of inertia Is =1, 226 in. 4 and | of eight imaginary beams T, to Ts | i Th 


which are assumed to be rigid, I= 0, “the sal 


_ As explained in Section 1, such a system can be replaced by girder Ts alone, | for the 


carrying (in addition to the nates masses) further i imaginary masses which ar are” ‘numeri 
termed the dynamic reactions of the omitted girders. 


There 
Fig. 15(a) shows girder Ty and ‘its ‘mass distribution. _ The weight of the higher 
girder, 65 Ib per ft ft corresponds to a mass m The addi- | introdi 

_ tional me mass at point 1 is the sum of the: dynamic reactions | of the beams T,and J = 


which been found equal to R, B= 1.26 ~sec? 4 The ad- 


ditional masses at point 2 are R2 = Re = 4.15 ——. . B By symmetry, | ny, the ; 


Substi 
_ masses at points | 3 and 4 are equal to those at at] points 2 and 9 and fix 


the computation, | concentrate the comparatively small distributed mass | that 
t 


at points 1 to 4 and thus obtain the mass distribution shown in Fig. 15(6): write « 

= M; =5 X 0.168 + 2 x x4. 75 = 10.34 by gir 


io of the small number of masses the e use of influence numbers (ex- . - 


> 4 plained in Section 3) gives the quickest result. For the present only the: modes 
that are symmetrical to to axis a—a are of interest; and to take full advantage 
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of this the influence numbers are defined in a special 


Numbers 6: and 612 are the deflections at points 1 and 2 as a result of a 
symmetrical arrangement of unit loads shown in Fig. 15(c). . Because of sym- 
metry, bis = O12 and = Similarly 621 and are the e deflections due to 
the two unit loads shown in Fig. 15(d). The numerical values of these influence » 
numbers are (in inches per pound): bu = = 0.108 X 107%; = $9 = = 0. aed 
10%; and = 0.274 X 107 


The amplitudes of the principal modes are y; = = using 
the same reasoning that led to Eqs. 56 i in Section 3 = nine are obtained — 


‘numerical values for M and 6 in Eqs. 56, oe 
(0.363 x 10-* — 1) 


0.558 X 10-4 + ( 2.885 5 10-4 — 


Eqs. 76 have only the trivial solution 4 yi = yo = O unless their detern minant 


vanishes. Therefore, A =0i is the * equency equation: 


-3. 198 — + 0.071 = 0 


There are two roots: —> = 3.176 and 0.022 or w w = 56.1 and w = 725. ‘The | 
_ a of w is ‘beyond the limit w = 100 and only one root. remains: 

= 56.1 and wy, = 3 149. To find the amplitudes of the normal mode, 
olen: this frequency in Eqs. 76 both | of which give the identical result 
yi = 0. 610 Yo A further relation is is the normalizing condition (Eq. , 58) W a > 


for the present case, reads 
Me + Msy*s + Mays 
ituting y: = = y, = 0.610 yrand ys = Yo, 2 (3. 36 X 0.6102 P+ 10. 34) y ys = 1; 
and finally y2 = 0.207 and y; = 0.610 X 0.207 = 0.127. To express the tact 
that these are ‘the amplitudes of the normal modes of system A one may 
write $4(1) = = 0.127; and ¢a(2) = 4(3) = 0.207, 
_ These are the values of ga at the four points where concentrated masses 
“occur. | Ther next step is to find the ‘amplitudes. of the mode ba for points situ- : 
“ated along the secondary beams. Since in system A these beams are rigid, 
they must also remain straight. At the ends where these bea ms are supported 
by girder T, the amplitudes « are eda = = 0.127 for girders 7. T. Ts, and Ts and 
¢4 = 0.207 for girders To, Ts, Ts, and T7; at the other ends oa must vanish. 


Therefore, the linear expressions | are: For beams Ti, , and 


= 
= 
| 

7 
7 
| 
| 
es 


1, for beams ams Ts, Ts, T 


In Egs. 79, z is the coordinate of the ‘point considered measured — the 
external | support of each girder. © 
Computation of the Coefficients | Bix. —The frequency eq equation, Eq. 32, con- 
tains the coefficients Bi which are yet to be found. Eq. 26 defines these 
coefficients. Intl this example one mode 2 pi and eight modes dis (k = é to 8) 
occur. There are, therefore, eight values of Biz (k = =1 to 8). 
‘The integrals in the equations for Bi, must be taken over all ‘the masses of 
‘the Because system consists of a a number of separate any 


Substituting « oa = 0. 1272 — (which is valid for points on beam 


and remembering that dm = m dz, 
=0. 127 fF giz dm = 0.127 | m (81) 
‘The value of the last integral in in Eq. 81 has already be been found, together with the 
“mode $iz, and By, = 0.127 X 0.874 = 0. 111. The: values of Biz for girders Ts, 


T,, and Ts must bet the same as those for Ti; thus, Bu | = =Bu= = allt = = Bs = = 0. ati. 


In the s: same manner, Bie = = = Bis = Bis = Biz fo = 0.207 | 


207 F don = 0.207 X 2.034 = 0.423. 
we, 
“The Equation.—If the numerical values and the values of the 
frequencies = 3,149; = 4B = = 8, 836; wep = 


= wes = wr = 2, 980—are e substituted into Eq. 32, the frequency equation 


"assumes the form: 
0.049 0.716 _ 3, 149 _ 
«8,836 2,980 


—1 
ieee 16, a w as variable, gives the values of the left-hand side of + Bg. 82 
as solid lines and the values of the right- hand one dotted | lines. _ The inter 
= points : are the frequencies of the system. . The roots of the frequency 
Raton found by trial error, are: = 1 1,654, 7, 950, and 26, 000 or 
7, 89. 2, and 161. _ 


| Papers O 
: 7 
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ant 
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| 
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«dit has been shown in Section 1 that Eq. 32 does not give all the iememnaion 
of a system if some of the values wiz are alike. This occurs in this example: 
W1B = 4B = WB = = 94.0 and = = WB = W7B = 54.6. 
multiple frequency of system | B remains a frequency of the complete system | 
and there are, therefore, two additional frequencies, w = = 54.6 and o= 94. 7 
Al the frequencies found correspond to modes that are symmetrical to 


axis a- —4, , and these frequencies, arranged in order of increasing magnitudes, - a 


4 vit 


are: wi = 40.7, ws = 54.6, ws = 161, 


To enable the designer to appreciate the meaning of these results the ampli- 
tudes of the modes 8 corresponding to these frequencies were computed. © This | 
Procedure 4 will give some enlightening information on the frequencies we re 


ws which would otherwise not be quite clear, 


@ 


7 


Determination of the Mo odes. and ~Te compute | amplitudes of 
x these modes which | belong. to the the Toots of the frequency equation (Eq. 32), im 


| Eqs. 33 and 35 were used. 


‘The amplitudes of the modes o1, ds, and 5 are illustrated i in Fig. 17(a), 


‘17(0), and 17(e). . Referring to the axes a—a and b— -b shown in Fig. 12 these 


modes are symmetrical not only for axis a—a but also for axis b— b. a 


Determination ¢ of the Modes and berm —To find the modes de and refer 


‘to! Eqs. 27 which now become 


& 


+ Buz CkB = = 


Eq. 83b eight equations | 


i ges 
| 
| 
— 
| | 
the 
tion 
(82) 
ency | Be _ 
ps 
e can be any number from 1 to 8. 
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Investigating the root = Ea. gives (for = 1, 4, 5, or 8) simply remai 

Bre ¢ = 0—or ¢ ca = = 0. For k= 2,0 the pr 


and, “because ca = 0, Cop = 0. Similarly, = CoB = cB = 0. 


; these results in Eq. 83a (remembering that Bu = Bu = Bis = Bis), only four 


terms of this e equation remain; thus, + + CoB + 3B = Since the 
‘mode is symmetrical to the a—a axis, ,C1B = and = = = or CB + = 0. 


The numerical values of Cip and csp are ‘computed from the normalizing 


‘condition (Eq. 35) = = 0.500 and sp = con = — 0. 500. ampli- 
_ tudes of the mode are illustrated in Fig. 17(d). A similar computation for 
= gives analogous Tesult—ca = 0; = Cap = CoB = Cop = 
CB = = 0.500; and cep = = = 0.500—as illustrated in Fig. 17(b). 
These two modes (Figs. 17(b) and 1 17(@)). are still symmetrical to axis a—a 
but they are asymmetrical to axis b—b. The significance of these two modes 
now becomes quite clear; each group of secondary beams can oscillate without 
_ interference » by 1 the other beams, because the girders on both s sides s of To move 


a. in opposition; girder Ty ‘remains stationary ‘because there a are no loads « on 
= -_ To, which is is due to the fact that the reactions for each pair of secondary 


+ 
ll 


beams T; and Ts, Tz and Ts, etc., balance each other. 


TABLE 3.—CoMPARISON OF FREQUENCIES 


40.7 | 54.6 158 
41.0 4 55.4 5 166 
= 
— 
7 Reliability of Frequencies.— _ —The me mechanical sys system em shown in Fig. 12 includes 


some distributed masses. _ Therefore, it has an infinite number of degrees of 

freedom. Such a a system has an infinite number of frequencies as well. — The 

"present study has revealed only five frequ equencies corresponding to ‘symmetrical 

modes because the frequencies of the complementary systems greater than 

= 100 have been neglected, thus : artificially reducing the number of degrees 


the: computations; and 2, the frequencies obtained by 
= the number of modes of the com plementary systems. . The first two frequencies 
unaltered, which i is’ s better than could be 


— 
@ 
— 
| 
| 
Th 
(a) 
(c) 
a chec! 
to the rule given in Section 1, this is a good approximation as long as the result | “ tain 
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remains below the lowest frequency of the secondary girders thems 
the present example this lowest frequency i is wee = 54.6. « w1 aye 
are good approximations, whereas the higher frequenci ies are not reliable. As_ 


(a) Normal Mode w,=40 7 


(1) ‘Normal Mode 54.6 


Fie. 17.—(Verticat oF d, AND ¢, REDUCED) 


a check on on this rule, . line 3, Table 3, shows the result of a more exact computa- — 
tion taking into “account the effect of the flooring. © Line 4, gives the result 


obtained if the same number of modes is employed as for line 2 2. The more = 


ng 
ult 


exact Tesults i in 1 line 4 differ from line 1 only by 0. 75% for | the e fundamental. | remail 

fn and by 1.75% for the next frequency, which is quite ‘sufficient for | furthe 

all practical purposes, and justifies the rules given in Section that n 
In the foregoing computations all the connections of the » 
have been assumed entirely hinged . In actual practice these beams will be at 
least partly restrained and this restraint will tend to increase the frequencies. ‘ 
To obtain an estimate of the e effect of partial fixation, the fundamental fre- 
“quency of of the system with fully continuous ‘secondary beams was determined 
by the method explained i in Section 2 2. | The fundamental frequency was found 


to be 43. 8—8% higher than that continuity. ‘The effect of partial 


restraint, therefore, will be less ‘than mathe 


‘machines run run at 320 r rpm, the circular frequency being w = 2m X 5.320 33.5. Ori 


~The factor of safety against resonance therefore is 55> = 1 This result is 

not quite satisfactory. Considering t that a certain margin n of ¢ error in computa- | 
tions of this kind must be expected, on one should attempt to obtain a factor of 
"safety of 1.40 whenever possible. ‘The structure shown in Fig. 12 therefore i is 


a“ . not quite suitable for the contemplated | purpose and should be strengthened. 


Th 
nant f 
of the 
‘Tesona 


1.40 X 33.5 5 = 46.9. 


__ Alteration of the Moments of Inertia.— ed moments of inertia re- 


een the change in J, which i is approximately correct, it will be seen that aie 
all the original computations can be used. 


To The most economical method of obtaining the small required increase in | 

the frequency is to | . alter only the beam, or * beams, with the lowest frequency ) 
wa OF ‘Thus, in the present case only beams T, Ts, Te, a and T; need to be 


= 


altered. This « alteration does not affect the modes and frequencies of system A. . 
the modes ¢13, $48, $58, and and their frequencies which originate 
from girders T,, Ts ‘Ts, and Ts, remained unaltered. 

obtain the new modes bop, $02, », and and their frequencies, 

‘theorem II, Section 1, is used d. It: states that the normal modes of a system 

‘Temain unaltered if all the m moments of inertia are changed, proportionately, 


from I into al; such a change, however, alters the frequencies w;, which become - 


If theorem II is is applied to each of the separate girders Ti, To Ts, and Ti, 


heir moments of inertia are | increased from to a@ Is. This increase alters 


their frequencies Va W3B, wep, and Vaw a w7B, but leaves. the 
bop, and unchanged. As the numerical values of 
none of the modes have varied, the numerical values of all the coefficients Bie 


) equati 
ety 
stituti 

ap Tro 
Th 
of the 
quence’ 
— 

4 4 

= 

— 
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remain the same. Thus, t the frequency equati 

further computations because it is only the , 3B, Wes, and 

that need to be. adjusted: — 


i 


a 
| 
| ‘The « e effect of the « change of sections on the frequency equation is only the 
Pires of the coefficient a. . Since the desired frequency w is known, this 
equation can be reverse to determine the required value of a. Sub- 


stituting 46.92 = 2,200 i in 1 Eq. 85 gives 3,836 +5 980 2,200 


and @ = 2.01. 
Originally, he 122.3 and now ‘it is ‘required hel 122.3 
=246in3, 


Conctusion 
The method presented in this paper permits the determination o of the reso- 
nant frequencies of complicated structures. It makes possible a ‘comparison 


of the values of the ‘operating frequencies of Machinery with the» values of 
resonant frequencies of a supporting floor. Ifan any y danger of resonance is found, a 
appropriate changes can be made before 1 the beginning of actual construction. 
‘The author has applied the method successfully to the frequency analys sis 
of: the complicated fan floors of boilerhouses and, on one occasion, to the fre-— 
quency "analysis of a high radio tower on the roof of a multiple- “story | building. 
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“ROADS AND PAVEMENTS, SAMPSON NAVAL 


TRAINING STATION | 
JACOB FELD," M. ASCE 


During World War II, the design and of temporary access anc 


permanent roads, walks, and paved training areas was one phase of the emer- — 


gency provision | for the United States. Navy personnel training facilities at > 


absence of a any local suitable base or aggregate materials, limited transportation - 
facilities into the site, and unusually severe weather conditions made necessary _ 
“flexible program. The construction methods for almost 65 lane miles of roads 


within 6 months, with little time for - preparation, a and with keen ‘competition - 


for labor and equipment, and yet at a reasonable cost, are recorded because 


| Sampson, N. Y. ¥., 01 on the east shore of Seneca a Lake. _ Unfavorable soil conditions, 
“of the extremely low maintenance required in 4 years of concentrated use as 


1 a training | station followed by 3 years of use as a temporary college. ‘Some of —_ 
the procedures ¢ deviated considerably fi from pre eviously established practice and 


are suggested as as proper methods of design and construction for roads in fine- 
“grained areas. Design predicated on the availability of 


; 


‘ Seneca Lake in New York, w ith about 5 miles of lake frontage, was to be 
converted into a training station for Navy personnel. The basic training areas — 


7 “Eat Extent of of Project.- —An area area of 2, 500 acres, adjacent | to the east shore of | 


(six i in number) consisted of a rectangular 15-acre sodded drill field, with a | 

paved drill area 50 ft. wide on three ‘sides. All six fields with barracks and 

service buildings were connected by paved roads and paths. — Each area had a 

an average of 1.5 miles of main road and a tr aining strip about 0.5 mile long. a 

A number of special service areas (for advanced training, administration, ware-_ a 

housing, civilian housing, and a hospital) a accounted for almost as much paving 


as in the s six training a areas. faced. ‘The 
aced. The 


a As a total, the equivalent of 51 miles of roads 22 ft wide was surf 
Toads ounpelsed (in million square yards): Primary road paving, 394; secondary 


_. Nore.—Written comments are invited for inenetione publication; to insure publication the soll 
discussion should be submitted by March 1, 1950. 


Engr., New York, N.Y. 
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ROADS AND 


PAVEME NTS 


paving, 158; concrete p paving, 48; ‘gravel surfacing, 17; and resurfacing 
existing roads, 49. About 95% of the paved | area required an 8-in. stabilized 
gravel s subbase, which with ‘shoulders involved : almost 300,000 tons of of gravel. : 
Cr ushed stone base was placed on all main roads. Hot lay tar conerete 
to surface all p primary, , secondary, and ‘existing road areas. Penetration 
seal coats required about 800, 000° of tar, and the hot mix required 


Land was acquired in early | June, 1942, before any surveys or plans had 


wees 


been executed. — Construction was well under way by the end of July, and i 7 


spite of weather conditions, all road excavation and subbase Ww ere completed 
by December 15, with about 60% of the tar | concrete i in place. The ‘station — 
was “established” on October 17 and all areas were occupied by the end of 
-= _ Some 440,000 sailors received their basic training at this station. 

_ Local Geology and Topography.—The entire e work was within the strip of 
‘eines shale exposure (Devonian period) (see. Fig. 1 1). | . The surface soil is a 
Dunkirk silty clay | loam (New York State “Agricultural Classification), 

~ designated as a Honeoye. silty clay loam. It seemed to have only two qualities 
-—tough gummy clay when wet, and gaseous dust w when dry. Soft shale was’ 
found e everywhere within a depth | of 3 ft. This soil i is very difficult to drain, 
‘quite unstable when wet, and readily subject to frost heave. | It gives a slightly 
alkaline reaction. _ Known also as “Seneca County clay,” it is - considered the 
worst soil for highway work in New York State. - nial is: 


The developed land strip was about 6,000 ft wide, generally sloping uni- 
formly toward the lake shore. Fig. 2 is @ part reproduction of the 


containing the area occupied by the Sampson Training Station. The edge 
— 00 ft has a drop of 60 ft, the next 1,500 ft has | a drop of 40 ft, and there i is a 
drop of 100 ft in the a telly 4 4,000 ft. The land was divided by a series of 
sharply sloping gulleys, running normal to the lake shore, approximately 2,000 


ft apart. "Natural runoff had eroded these channels from 20 ft to 25 ft in depth; 


less than 30 n min dating a thunderstorm. "Runoff was almost 100%. V egeta-_ 
tion consisted of some tile drained grape and apple orchards, and small isolated 


fields of flax and buckwheat. Most of t of the ar area, however, was: abandoned su! sub- 


marginal farmland, 


Any rainfall in excess of ( 0.1 1 in. stopped. rubber-tired vehicles on unpaved | 


"areas, A heavy rainfall stopped all operations for 3 days. — At all other times, 
continuous sprinkler wagon coverage on all traveled ways was absolutely 


essential During dry periods, the dust cloud | was visible 10 miles away, | and 
was not. ‘unusual for automobiles to drive with full headlights in 


daylight. 
Weather Conditions.— —Local experience indicated that: most of the work must 
be performed before the usual, heavy, early fall rains. The year 1942 did - 
a agree with local experience, however. — In July : and again j in August, there were 
10 days in which precipitation exceeded 0.1 in. . In September and October, 
there were only such da days each m month. Normal construction season in this 


quadrangle Geological Survey (United States Department ‘of the Interior) andl 


and, as far as possible, were left undisturbed. With the area 
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area is less than 6 months; the last frost, was on May 11 hi, and winter arrived 


ConsTRUCTION MATERIALS» > 


Available Materials.—In the immediate vicinity, the. only stone found was a 
soft shale, which slaked easily and dong oy upon exposure. | _ The nearest hard 
shale quarry, near Hayt Corners, N Ay + # ¥. (a 10-mile haul each way), | had been 
abandoned and required a full ‘complement of tools and equipment. Lack: of 
time prevented its development. _ Good limestone and some dolomite. was 

4 available from Oaks Corners, N. Y., 6 miles north of Geneva, N. YY. (a 22-mile 
haul to the station) . This rock was selected for railroad ballast, road base, 


7 nel concrete aggregate, and tar concrete : aggregate. way 
it 


production. few days of work depleted the entire area. Good sand 
_ gravel deposits were being worked along the contact between the limestone ~oe 
shale outcrops, chiefly at Oaks Corners and Phelps, N. Y., about 24 miles 


The Oaks Corners pit produced a , clean, fine, n —double washed. | 
- The glacial deposits at Phelps gave a satisfactory run-of-bank sandy gravel | 
suitable for road subbases, with no processing | or blending. A fine sand was 
ee from a deposit at Junius, N. Y., containing considerable lime in- 


-trusion from outwash seepage. Adjacent to this deposit w: was & ‘sandy y gravel 
- deposit: (54% retained on the No . 4 mesh, and 4% passing the No. 100 mesh) — 
which was developed to produce 2,500 tons per day. ‘Transportation bottle- 
rf necks and depletion of exposures in this area north of Geneva forced | 
toward the south, and in 1 October, a run-of- -bank | gravel pit was opened at : 
North Hector, N. Y.—about a 16-mile haul to the area. Several feet of over- 
ow was blended into the gravel and much large stone had to be “sealped” 
- = the e roads when this n material \ was spread. The finer graded materials - 
were more suitable and ¢ gave better results with less effort. uns a 
_) Transportation Facilities. —The only road feeding this area is New York Toute — 
-96A, which was then a concrete pavement 16 ft wide. The lake road (see Fig. 
was a black-top light. traffic -pass branch of route 96A. Railroad connec- 
ao were built as an early part of the program, , but there on little space on 


cars for anything beyond lumber and metal materials. 
“a By careful scheduling of truck 1 movements over the one one road : access (using 


about 110 trucks on the day shift and 50 trucks at night) n more than 100, 000 
tons of sand, stone, and gravel were | brought into the site per month. | To 


permit higher speed and to increase safety, the road was widened by paving» 
with one course of hot lay tar concrete a 3-ft strip on each side of the 16-ft 


concrete str: strip. No trucks were permitted on the lake road, and trucks -— 


required to. keep off the main road for 30 : min, morning and evening, during ~ 

the peaks of worker travel. _ Almos ost 10,000 people, about 90% in private — 

(fully loaded), came into the station within a 30-m min period a’ at the the beginning — 
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— 
of the day shift. About 8, 500 sven ti from th the n north, with tra affic congestion of 


Existing Roads and Temporary Facilities — —Three east and west roads were 
found crossing the station area (see Fig. 2). | Road A was & built-up surfaced 
farm road, consisting ¢ of a cobblestone bed 8 in. . deep 0 over part of the width, the 7 
remainder being 6 in. of loose crushed soft shale. On this bed had been placed 
several layers of bituminous penetrated crushed limestone (probably a tar 
binder and 1 in. of unscreened stone), making u up 2 a 6-in. thickness. — The total 


surfaced width was 16 ft, bounded by steep ditches. 


~~ Road B consisted of a 6-i in. bed of soft shale with a 2-in. ete 


trated screenings. was. also 16 ft wide. C was similar but the 
screenings had “not penetrated. Locally, these roads are known as 
“shuck” ‘roads and ; are made by spreading crushed shale directly o1 on the » natural 
clay oll and letting traffic iron out the e aggregate into a pavement. _ Later, 
screenings are added to fill out the ruts; and, when traffic warrants it, especially 
to reduce the dust in dry periods, | a ‘penetrated surface is added. The roads 
shed water ' very r readily, are not slippery, and permit a 35-mile speed fi for all 
types : of passenger and light truck tr affic. Roads of this type were used for 


: truck loads of sod (from 3- ton loads to 5-ton loads) 


pote yre these Toads required dig digging new ditches and filling i in the existing 


drainage channels. There was no time for change. Ditches were 

dena and the road surfaces were scarified 1 in. deep, penetrated with grade 
RT-4 tar, and covered with } ‘in. of stone screenings. — No added surface layer aa 
could be applied because the resulting road width decrease was not con- _ 


sidered safe, 


a In the middle of each road ad length, a a traffic circle was constructed to to permit 
direction. 7 The radius of this circle was 30 _ These traffic 


_ Eventually, after construction was almost completed, these roads \ 


“the gravel subbase for use as temporary access continuous main- 
7 tenance ce gang of from 10 to 25 men, with : a grader, a a a bulldozer, a and 2 trucks kept : 
all all roads available for construction traffic, while also operating 3 water sprinkler 
"wagons to keep the the dust within some control. neon 
 All-weather par parking and bus terminal areas were built by pr preparing sand- 
clay compacted surfaces. The procedure wi was to searify the soil about 3 i in. 
— deep, , cover with a 2-in. layer of fine dusty s sand from the limestone outwash 
area, harrow and disk the — mixture to a uniform color, and roll the surface. 
‘Wate was added from a st sprinkler truck. = Control of water xr and of rolling was 
in accordance with a simple instruction: ‘When a light private car left tire marks, 
- More e rolling w was needed; when the t tires raised dust, more water was needed. 
These p parking areas wwe ‘maintained by adding stone e chips as low s spots de- 


There was no subsurface drainage and the surfaces were sloped at 
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~~ terminals still left fairly good condition. 


= Grading ) Excavation.—In general all excavation and shaping of the subgrade 
was by tractor scrapers; fifteen such units of 12-cu-yd ‘size, aided by 60- hp 
- bulldozers in tough ground areas, each moved { from 600 eu yd to 700 }eu yd 


daily to build up up grade. and cut ditches. es. The wo work was then finished by ps patrol | 


_ graders so that practically no ‘euboulk: was needed. . The road 1 cuts were first 
made 5 5} in. below final grade between « ditch slopes : and nd then * ‘boxed”’ out 8 in. 
for the gravel subbase. . For primar ry roads, approximately 500,000 sq yd 


f grading was was required; and, for secondary road areas, , 200, 000 sq yd. :. Peak 


3.—Gaavme Wonx AFTER A — 


production in the ait September when 143,000 sq was 
graded and prepared for gravel subbase. Most of the excavation ‘volume 
pony from the drainage ditches. This. excavation was expedited to provide 
drainage for all road excavation. To drain the box trench excay ated 
a the road pavement bleeder trenches were cut 75 ft apart ¢ and 2 ft wide, on 


each side into the open ditches, and filled with broken stone. Working condi- 
tions after a rainy spell are shown by Fig. 3. 3. 


The same -labor-equipment group prepared subgrades, finished 
“ditches, and nena gravel subbase sige 2, 000 tons per day to 4, 000, tons 


150 men to 185 1 men foremen. 


-.) about 1 on 200 towar as and ; b 
ig 
fr 
4s 
jo 
st 
22 
(2 
at 
esi 
og 
| TO 
d = 
194 
| 


October, 1 1949 AND PAVEMENTS 1139 


; To expedite the road preparation work, some wheel- ‘ieee scrapers were 
_ brought to the site in October because no additional tractors were obtainable. d 
The soil conditions ‘made these : scrapers useless s since the ground was too hard to 
cut. without the help of double bulldozers and the slightest rain eliminated all 
traction. i. ‘The additional equipment was then used chiefly | for ‘spreading gravel © 


from stock pile and other comparatively easy sy excavation. —_— sm. 


48 in. average ‘diameter being 203 Pig Small sizes were laid with open 
joints and trenches were backfilled with broken stone and large gravel. _ The - 
total length of all drain pipes used was 71,000 ft, with ‘about 10, 000 tons of 
stone fill. _ Maximum production in pipe laying was 25, 600 ft in October, 
about half of f which was laid i in trenches dug by back hoes; | the remainder was 
laid in sections of o open sn ditches prepared by the road ga gang . In both Septem- | 


ber and November more than 15,000 ft of pipe waslaid. 


~ Pavement Design —The general schedule for - road construction planned all 


- permanent road work for September, October, and November to permit the 
buildings, railroads, and subsurface utilities to have a good start in June, July, 
and August. Although the size of the project was almost doubled, no extra 
time for 1 road construction was possible because of weather conditions. Some 
Work ¥ was } performed 1 in December, but about a third of the final road surface 


Primary o1 or main roads, consisting of 8 in. of gravel, 24 ft wide at ‘the 
a 3-in. stone base, penetrated? with RT-12, 24 ft wide at the bottom; a 1}-in. 7 


binder « course and a 1-in. wearing course of plant mixed tar concrete, 22 ft wide; 
A and 6-in. gravel shoulders, 9 ft wide on each side. _ oii 


b. roads, paved 11 ft for affic with the six 


cconerete. On narrow roads, gravel were ft wide o1 on each side. 
pavement on paths and light traffic areas, consisting ofa6-in. 

eons subbase; a prime ¢ coat of RT-3 (0.4 gal per sq yd); and a 1}-in. binder 


course (rolled with pup’ ’ roller) of plant mixed tar concrete. 
All subgrade excavations were crowned, allowing for for a final 3-in. crown on the 


_-22- ft-wide roads. . To discourage high speeds, no ‘superelevation \ was provided 
at curves. _ The foregoing designs were agreed upon late in July, after alternate 7 
estimates of cost and labor requirements eliminated the possibility of concrete | 
roads, , and material controls p prohibited the use of asphalt in any form. 

7 Subbase Course. —Daily subbase gravel deliveries averaged more than 2, 500 


cu yd with maximum deliveries at about 4,500 cu yd per day in October, so that: 


2? “Military Roads in Forward Areas,’’ ASCE Manual of Engineering Practice No. 23, New York, N. Y. 
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stock pile « of more than 32,000 tons was ready for use when weather slowec ed both 


§ gravel production | and trucking, 
Bids had been asked for the following items: 


1, Run-of-bank gravel, ma: maximum s size 23 in., with more e than 30% p passing 


the j-in. screen; 


2, Run- n-of-c rusher limestone- -dolon ite, passing a 2-in. screen to 
3. Run-of-c usher hard shale, passing a 23-i 

48 — that will not weather or ‘soften when ‘exposed to air; and he 
4, 


‘Slag from industrial furnaces, all passing a 23-in. screen. 


a: Each vender was to agree to furnish. 50,000 cu yd and 1 to deliver ata rate of of 


3 ,000. cu yd per day of 20 hours. A theoretical preferred. specification was 


"Sieve size by (% 


10 to 30. 
7 The —_— passing the No. 200 mesh was specified to be between 5% a: and 15%, 
less than half of the fraction passing the No. . 40 mesh. — The t binder includes 7 
“any blended binder passing | No. 40 mesh with a liquid limit less than 25 and a 


ion at each pit 

and the acceptability of the gravel shipped w: was determined by the men placing: 
it into the subbase. Considerable gravel from different ‘sources was blended — 
- at the site to reduce the necessary blading and rolling for proper compaction. 
The cost varied from $1.45 per cu yd to $1.75 per cu yd delivered to the site, — 
“measured ‘in trucks upon delivery. As the previously developed pits were 

depleted, prices went u up icaneneeiie for the e purchase of land. and for stripping | 
and exploratory expenses. Hauling costs remained stable during the entire | 
period at approximately 5¢ per cu-yd mile hauled, with no allowance for the 
—— trip. 7 Practically every available truck within 50 miles of the station 
employed int this haulage, as subcontractors to the gravel venders. As 
an average, a ton of gravel yielded 2.67 sq yd of subbase, compacted toa 

of 8in. This is equivalent | toa density of 125 1b per cuft. 

~~“ On main and secondary road areas, a layer of fine sand was first spread, 
a from 1 in. to 2 in. thick, to seal the clay, before placing: the gravel unless the 
gravel as delivered d contained a surplus of fines. In = training ar area built, 
the subgrade ‘was covered with straw cut from the area. This straw was laid 
about 2 in. to 3 in. high and the gravel was dumped on it. _ Objections to this 
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procedure rewalted in the | change to fine run-of-ban silty as a 
“on the clay subgrade. __ Inspection of one of the earliest roads b built, after 1 year of f 
use and after 3 years of use, indicated that thes straw had served its purpose. No 
clay had worked up through any of these roads. — The other roads, with the 
_ sand screen, also stood up very well; curves, grades, and crown remained in 
excellent. shape. e. The volume of the sand used is iadkeded ta in the quantities. 
listed as gravel. = ‘The satisfactory | behavior of the  subbase i is partly explained 
by the « screen layer er placed on the subgrade. Shoulders of gravel remain in 
good shape (Fig. 4) although weed growth wa was not controlled. ow 
Stone Base.—The contractor began work i in the stone base and penetration 
7 on September iL. 1942 (Fig. 5). . The stone used was ¢ crushed dolomite, delivered 
by rail and distzibuted from the siding by trucks. The grading specification 
100% passing a 24-in. screen, , 90% to 100% passing a 1j-in. . screen, 
and 0% to. 157% passing a 1- -in. screen. After | spreading ng and rolling» to the 
r RT-12 depending upon air temperatures) | 
was applied by pressure 25 gal per sq yd. The surface was then 


covered with 50 lb p per sq yd from tin. stone to 3-in. 1. stone, ne, rolled, and broom 
bed. This was followed by a second 


After a week of trial, it was s found expedient to apply | the stone base i in two 
7 layers, which reduced the total rolling time. On | November 20, the tar was 
Seas changed to RT-10 because of low temperatures, and work continued through 
1 7 most of December. — The original plan was to allow an interval of at least a 
a - week between the completion of the gravel subbase and the a application of the 
stone. After November 1 (with: freezing weather starting | on October 27), 
with the danger of rain freezing within the gravel, ‘the stone was: applied soon 
after the gravel had been shaped. required some extra stone, since 
traffic had usually pounded out some depressions in the gravel subbase which 
Total: amounts used approximated | 60,000 tons of No. 3- -As stone; 10, 000 tons 
of No. 1 stone, No. 1-A stone, and No. 1-B stone; 1,000 tons of No. 3 stone for 
7 drains i in shoulders; and almost 700, 000 gal of tar (RT- 12, ‘RT-11, and RT- -10). 
- typical gang consisted of 30 men with 2 spreader b boxes, 2 rollers, and 2 
shapers, and they used about 400 tons of stone to « cover 900 ft of ‘out 24 ft 
= At full speed, five gangs were ‘operating. - Stone was ‘delivered to the 
work by | 15 trucks, and 2 cranes unloaded t the cars. _ The stone cost 27¢, the eo 
placing: cost 70¢, and the tar cost 21¢ per t per sq yd of base. . The: stone delivered 7 


to the site on cars cost $1.68 per ton, 


Black- -Top Pavement. —A premixed t tar co concrete test: road section, | 12 ft by 
250 ft long, was laid « on | a . gravel subbase to test the feasibility « of this material — 


for the intended traffic. Even though no base was laid down, heavy trucking 
did no not seriously damage the test strip. Upon deciding to use tar mix (asphalt 
was not available), it was 1s imperative to move a st suitable plant to the site. _ Ape ~_ 
— 75,000 tons were required and the — of traffic on ‘the 
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of hot 


y access roads would seriously handicap | any attempted delive ery. schedule 


located, and it was moved bodily to the contractor was 
engaged to » peodeses all binder and top course material at $7.65 per ton yn complete 
in place. This price covered all materials, equipment, and labor, e except that 
a railroad siding was ‘provided to the contractor without charge. _ About it 20 
men including truck drivers (the | average haul distance was 3 miles) were 
The plant | was an integrated mixing unit with feed bin and reciprocating 
feeder. 7 Two aggregate driers, equipped with oil burners and a cyclone dust. 7 
collector were set up in tandem, with hot and dust elevators. Screens i in the 
mixing ng tower unit were 2i in. » in., » and $ in. The mixer was a a4 ,000- lb- 
capacity twin pug mill, equipped oth -5,000-Ib aggregate batch scale and 
~ 500- ‘Ib tar bucket scale. Storage of aggregate w as in a 50-ton bin with four 
compartments. Tar storage was in tank cars with a 10, 000-gal_ reserve storage 
tank. Steam generated by a 150-hp locomotive steam boiler was used in the 
ter heating kettle and i in the jackets surrounding the tar lines and mixer box. 
The mixing cycle, electrically controlled, was 1s 45 sec, including dry and wet mix. " 
Temperatures of the mix leaving the p plant were from 180° F to 200° F for the 
binder and from 200° F to 230° F ' for the top course. — Plant production aver- 


100 tons of binder and 80 tons of top course per hour, 


__ The tar concrete w was mixed, laid, and rolled according to New York State 
tions? — | top course were in 
with the following limitations: 


13 in 


The aggregate sizes used in the mix were of New York State highway desig- : 


nation and were of the description: 
y The No. 2 , No. i, L, and No. 1-As sizes, which were limestone, as well as the 


sand, or limestone screenings, were shipped by rail. ‘The sand was 
trucked in in from a pit 22 miles distant. proportions of. each aggregate 
Were set up to provide a impervious, surface c course, retaining nonskid 


Hot Lay Tar Concrete Pavement—Two Construction Mixing Method—Type 6,” Y. 
' State Highway Specification, Albany, N. Y., 1939. 
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he mixing after some trial proportions were laid and 


odified, were as follows (Table1): = | 


TABLE: 1 —MIxine ForMULAS For BLAck- Top. 


Nove: EMBER 19 9 NoveMBER 19° 


@ 


4) 6) 
a 1,560 
ie 


28 120 [ | 
6 


89 09 


a. substituted for RT-12. Check analyses were made daily of — 
from the bins and of the mixed product. — Where it was not possible, due to 
construction scheduling, to. cover the binder cour: course within 2 tha der 
mix (Col. 5, Table 1) was changed. This denser binder * mix withstood con- 
struction traffic better than the more open mix. <. The better practice is to cover © 
the binder course the same day with top mix, butl in order to keep traffic flowing, | 
and to avoid changing from binder m mix to ) top mi mix at the plant during 1 the middle | 
of the day, the binder course was left uncovered for periods : as long as 3 days 


or more in some few instances. . Wh hen cold weather arrived all the binder was 


of the denser mix. 
Tar n mix production began on September 15 and closed down on December | 
9 for the winter, after 40,000 os had been placed. An illustrated progress ’ 


‘report ¢ of this work to this period 1 was published? i in January, 1 1943. | The plant — 
was opened again in ‘May and completed the remaining 35,000 tons by the | 


middle of July, 1943. At first the work progressed on a 10-hour day, 7- day- — 
-per-week basis, and changed to a 12-hr day, partly after nightfall. — ¥ Only the 
binder mix was laid ur under artificial light. Production of 1,300 tons in a 12- 


br day was not unusual. 


unt _Navy’s order to stop all binder placing when the temperature fell 
¢ below 32° F and top course W pron it w was below “al F was only generally followed. 


yr October (with more than 2 in. 7 

of snow) and 41.7° F for November (with almost 5 in. of snow and 7 days Ww ith : 

precipitation greater than 0.1 in.). W ith below zero weather | early it in De- 

_cember, it was impossible to continue the surfacing. | 


__The contractor placed tart mix using 4 4 paving machines to spread and 


paving machines. which pushed | the trucks ahead. / Additional equipment 
used included 8 rollers, 21 trucks, 1 grader, 2 bulldozers, power broom, and 


Construction Methods, January, 1943, p. 62. 
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O h th 10- ‘ton to ‘12+ -ton, 3- wheel ‘rollers, or 8-ton to 10-— 
to. produc a compression of 250 Ib per in. of tread width. 


Surface tolerances were } in. from crown template and } in. from a 10-ft straight- 


_ edge parallel to the axis of the road. — It was pointed out by | "the vender in 1 in a. - 
middle of October that tar mix was being placed with air temperatures far 
below the 50° limit specified by the N. Y. State Highway Department and 


that the base was often i in rather soft condition. cases, a 2- ‘ton roller 


work was continued immediately At times, the 
"spreader had to be pulled forward by the trucks delivering the mix (instead of ; 
the spreader pushing the trucks) to pass through soft spots. In spite of all 


“these conditions, the road surfaces held up very well, the gravel subbase drained - 


‘Fic. 6.—Texture or -Roap ‘SuRFACE (Rien?) AND AssEMBLY AREA (RULE Is 8 In. Lona) 


out th the excess moisture, and the roads w were used for 4 years with practically 
no maintenance. _ An illustration. of the pavement surfaces photographed in 


1946 i is shown in Fig. 6. 


- Theoretical yield was 7. 0 sq yd per ton as compared with an actual yield | of 
about 6.5 sq yd. ‘The binder course of 13 in. compacted ¢ gave 10.93 sq yd per 

ton and the topping course of 1 in. compacted gave 15.61 sq yd per ton. ‘The 
‘density of the top mix, as tested from samples i in place, was 141.1 lb per cu ft. : 
The he yields include the losses in material at soft spots and the filling of rutsin 
‘the : stone base from traffic. The total cost of primary pavement surface was” } 
$1.08 per sq yd. 


Pavement. 


trucking and vehicle servicing was 
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4 with a concrete surface. Concrete design was specified in accordance with the 
State Highway (1939) -Tequirements_ for plain concrete, permitting 
. "crushed gravel as coarse aggregate. 7 Cement content was 6 bags of Portland 
and 1 bag (80 Ib) of natural cement per cu yd. . With a water content of 5.5 
8 al per bag of cement, the slump was held to 2 in. a When the central c concrete 
plant was not being run to full capacity by the other o operations, transit mix 
concrete used for road construction. To prevent stoppage of work 
— small batch bin was set up and a paver mi mixer, with cement in in bags, was used 
to maintain continuous | concreting. Nor reinforcement was used in the 6-in.- 
fy slab, with transverse joints filled with cane asphalt premolded strips, 
in. by 53 in., and transverse dummy joints: 1} in. deep | spaced not more than 


_ About 50,000 sq yd of concrete pavement w was nits cold weather 
"forced a change to black-top pavement for the remainder of the area. . Water 
for the mix came from Seneca Lake and during the period of water main steril- 
ization with chlorine, very rapid setting of the co concrete mix was noted. With 

a residual of 5 ppm of chlorine in the water, co concrete had to be ‘spread and 
finished within 20 min of mixing, e even with air temperatures ‘between 45° 


_ _ and 5: 55° F. Little ‘surface spalling | was found in this concrete e paving after 4 
years” of use and exposure. Paved surfaces were cured with straw covering 


> 


the — was done k was done by an an experienced concrete road gang. 


PERSONNEL ACKNOWLEDGMENT 


_ All work Was | ‘was under the direction of ‘Capt. J.C. Gebhard, M. ASCE, (CEC), 
wv. 8. Navy, with: the firm of Shreve, Lamb a and Harmon associated with Gilmore 
Clarke, M. ASCE, Michael Rapuano, | and Malcolm Pirnie, ASCE, as 

Te the roads and pavements was assigned 
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PAPERS 


AERODYNAMIC THEORY OF 


D. . B. S STEINMAN? 1 M. ASCE 


Oo o 


SYNOPSIS 


General equations are are derived, expressing the aerodynamic “forces of lift. 
and torque on any oscillating | section of two independent degrees of freedom. — 


The parameters. representing th the distinctive aerodynamic characteristics of the 


‘ection. and the coefficients ‘representing the variation. with | wind velocity ore 


derived from static wind- tunnel tests on straight and curved section models. | 
The effects: of vertical velocity, angular displacement, and angular velocity ae 


“ included, as as well as the effects of the angle of incidence.’ The complete e: expres- 
sions for the aerodynamic forces 2 are equated to the expressions 


“restoring force, inertia, and d damping. The resulting ; general equations yield 
solutions for all practical problems of vertical, torsional, and. coupled oscilla- 


ies, instability» response se (rate of ‘amplification), an on amplitade response. 


ples are given, to illustrate the practical spetention of the — 


‘dined sitions and to confirm the validity of the underlying theory. Graphic 
“methods: are presented to facilitate application and visualization. 


is of the general “aerodynamic” 
for the forces (L and M) acting © on any section, together with the general 
“dynamic” equations for the s same e forces. Once these ; general. equations ai are 


v velocities, , amplification ¢ or ; damping, limiting (maximum or t minimum) ampli- 
tudes, amplitude re response curves, and other relations for vertical, angular, and 


coupled oscillations. 


4 Nore. —Written comments are invited for immediate publication; to insure publication the last dis last dis- 


_ cussion should be submitted by March 1, 1950. — 


1Qons. Engr., New York, N.Y. * 
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For the complete solution, the aerodynamic equations must include the 


effects of vertical velocity, angular position, and angular velocity « of the sec section. 


dynamic equations must include the elastic (and gravitational) re- 


sistance to to ) displacement, the inertia resistance to acceleration, and the damping 7 


force or resistance Ge with the velocity). 
Theoretical equations for the limiting case of a flat plate or streamline 


ee have been derived by Theodore Theodorsen,? R. Kassner and H. Fingado,* 


Glauert, 4.66 H, Kiissner,” Herbert Theodor von 
M. ASCE, and J. ‘M. Burgers, ® and others, using tl the methods of classical 


aerodynamic theory (based a flow and the 2 application of the Kutta 


for which the Kutta ‘condition 1 is not applicable, it appears impossible 

to write the aerodynamic equations by purely theoretical derivation. - Some of 

the necessary coefficients must be supplied by experiment. “elegans 
~ Consequently, for a practically feasible solution to cover all cases, ‘it is 

necessary to deduce e rationally the correct | t general form of the aerodynamic 

equations, with the u unknown coefficients or parameters (representing the 
distinctive characteristics of different sections) to be supplied experimentally. 
: a”. Since only the aerodynamic characteristics of the section need to be deter- 


mined experimentally, full-scale models incorporating the dynamic character- 
istics of the complete | structure are unnecessary ; simple section models suffice. 


_ Moreover, to the extent. ‘that static wind-tunnel tests reveal the distinctive — 
- aerodynamic characteristics of a section, oscillating model tests are not required. 


thesis that. the distinctive aerodynamic ic characteristics of of oscillating sec- 


tions can be. mathematically predicted from | the known or measurable aero- 
ade characteristics of stationary section models has now received cumula- 


ive and consistent confirmation. 
Notation. —The letter symbols used in this paper customary notation 
as far as practicable. 
Vertical v vectors are measured positive e upward; angular vec vectors are measured 


positive clockwise; wind velocity i is to the right, and wind inclination is measured - 


As used in this paper, b is the chord or or width of the section (in feet), and - q _ 
the « dynamic parameters, m, and. K, are per full width. flow 
Theory of Aerodynamic Instability and the Mechanism of Flutter,’’ by Theodore Theo- Spac 
dorsen, Technical Report No. 496, National Advisory Committee for Aeronautics, Washington, D. C., 1935. 4 are 
3‘*Das ebene Problem der Fligelschwingung,”’ by R. Kassner and H. Fingado, Luftfahrtforschund, 
7 . November 20, 1936, pp. 374-387. (‘‘The Two-Dimensional Problem of Wing Vibration,” by R. Kassner the | 
a a Fingado, Journal, Royal Aeronautical Soc., Vol. 41, 1937, pp.921-944.) ya 
4“The Accelerated Motion of a Cylindrical Body Through a Fluid,” by H. Glauert, Reports and Memo- _— 
randa No. 1215, British Aeronautical Research Committee, 1928-1929, Vol. 1, January, 1929, pp. 118-127. posi 
: «$$ “The Lift and Pitching Moment of an Aerofoil Due to a Uniform Angular Velocity of Pitch,” by lin 
H. Glauert, Reports and Memoranda No. 1216, British Aeronautical Research Committee, 1928-1929, Vol. 2, 11nes 
January, 1929, pp. 636-644. 
6**The Force and Moment on an Oscillating Aerofoil,”” by H. Glauert, Reports and Memoranda No. 1242, 
- British Aeronautical Research Committee, 1929-1930, Vol. 2 2, March, 1929, pp. 742-758. anibia t 
%*Zusammenfassender Bericht tiber den instationiren Auftrieb von Fligeln,” by G. Kassner, 
8 ‘Ueber die Entstehung des dynamischen Auftriebes von Tragfligeln,” by Wagner, Zeitschrift 
angewandte Mathematik und Mechanik, February, 1925, pp.17-36. 
_ *“General Aerodynamic Theory—Perfect Fluids,’’ by Theodor von K&rm4n and J. M. B M. Burgers in 
Theory” (W. F. Durand, editor), Springer, Berlin, 1935, Vel. 2. ~ 
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Numeral subscripts distinguish corresponding parameters for vertical. and 
torsional oscillations: w:, Ki, 51, 81 (and 83) for vertical oscillations; Ky, 
92, 8 (and s4) for r torsional oscillations. ons. 


Use oF STATIONARY SECTION 


— of Curved Models. s—By a a principle of relativity, the arm of the secti ~— on 


~ (regarded as an immersed boundary surface) may be translated into a relative 
tilt or distortion of the flow lines in the fluid, or into a contrary relative tilt or 
"distortion of the boundary surface, so that in either case the boundary ool 


then be treated as as stationary. 


graphic representation (Fig. 1) ‘consists of the: simple composition of 


‘Fic. 1.—TiLT1InG oF Fretp or FLOW FOR 
VerticaL Motion or Section 


“Fie. 2.—Curvine or or For ANGULAR 7 


Morton or Section 
son oF Suction 


flow lines with spacing 1/V and relative (upward) vertical flow lines with 
Spacing 1/v. _ By this graphic composition, the horizontal flow lines in the fluid 


_are simply tilted up toa positive angle of incidence v/V. Instead of tilting 
- the horizontal flow lines upward, an equivalent procedure i is to tilt the section — 
“downward the same angle v/V; this is also given (Fig. 1) by the graphic com-_ 


‘position of the s same superimposed flow lines with the directian of the vertical 


» eal Analogy 1 —The tilting o of the stationary model in a wind d tunnel t takes 
a the place of a of the — 


When the same procedure i is 3 applied | to angular velocity of the section, the - 


graphic r representation (F ig. 2) consists of the simple co: composition of the hori- 
Zontal flow lines with spacing 1 L/ V and a series of concentric circular flow lines 


| 
a 
q 
— 

35. 
mo- 
a 
— 
| 
rift, 


with graduated spacing 2 3 b/(c ‘v)—in which vis is the linear ¢ downward velocity — 7 


at the leading edge and c is the distance of any flow filament from the center. : 
This graphic composition yields a series of concentric circular flow lines , convex 
upward. The angular velocity o of the section is thus translated into a curved 
field of flow, as in an arching wind tunnel, with the immersed section model kept _ 
stationary and unchanged . Instead « of curving the field of flow, » an equivalent 


procedure i is to curve the immersed section 1 model s an equal at ando opposite amount. 
This curvature is s also { given . (Fig. 2) by the graphic composition of the same 
The equivalent warping 1g of the immersed section model is / symmetrical about 

_ the center of angular rotation, with horizontal tangent at this center and with 
slope change v/V at each end; the radius of curvature is }b V/v. Accordingly: 


__ Analogy 2.—The curvature of the stationary model in a wind tunnel ; 


takes the place of an angular velocity of the section. 
With a curved section ‘model, : as in the more familiar tilted section model, the 


_ added 1 angle of of incidence at the e leading edge remains of) Ve 


#&F or wind- -tunnel | ‘duplication o of ‘this condition, t the section model m must it 


Ww arped to the circular (or parabolic) ¢ curvature indicated above, with horizontal 
_ tangent at the center of rotation and the curvature changed for different de- 
_ sired values of the » slope B at the leading edge (angle of ‘incidence at 1 the leading — 
edge equivalent, to angular velocity of section). > This curvature is called the | 
_ ‘The first w wind- tunnel | tests on curved models of bridge » sections were made 
in 1947 at the Vi irginia | Polytechnic 1 Institute at Blacksburg by FJ. Maher, er, 
in cooperation w with the | writer. The results yield the needed additional data 
for the analysis of torsional instability. 
In the absence of such tests on curved models, the analysis’ for the effect 
of angular velocity has previously been based on on the static torque graphs 
ob tained from straight section models (simply ti tilted “instead of curved). 
Qualitatively, ‘the results and conclusions should substantially unchanged, 
3 since the major disturbance of the incident stream flow is determined by | the 
7 > effective angle of incidence at the leading girder. Papin the indicated 
—— from straight section models to curved section models is in the 
direction of reducing stability torque or incensifying instability torque. ‘al 
od it should be noted that the straight, tilted section model has a dual signifi 7 


cance. re. It represents the effect of angular position, 1, a, also: the identical effect 
verten velocity, v. Accordingly, in aerodynamic formulas, the two 
interchangeable; the conversion relation i is simply a =— The two 


tributions may be bracketed, as (a — »/ with the same factors and ‘coefli- 


cients applied to both. 


The curved section model represents solely the effect of angular velocity, @. 

The conversion relation i isB = — b a/(2 V) 

a For. ordinary values of B, it is ; sufficiently accurate to use a circular or para- 

bolic are for the curvature of the model. _ For higher values of 8, the curve is 


“more accur ately : a catenary, to to keep the elements of the e section width unchanged : 
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in 1 length. For the limiting case of V = 0, yielding B = — 90°, the eatenary. 


becomes a folded straight line (zero span, and = b/2). 
a ‘Curvature of the section n model i in static wind-tunnel tests, as the equivalent 7 


of angular | velocity ¢ of the : section, is (for t this application) a new concept intro- 
' - duced by the writer. 7 A search of aerodynamic literature reveals that a similar 
q idea has been applied i in the study of curvilinear flight of airships. 10 The same 
idea, but intended to be applied analytically instead of experimentally, i is also | 
briefly suggested by Messrs. Kassner and Fingado,’ w ho state that one of the 
terms of the lift force due to angular velocity of an airfoil “could be deduced 
replacing the rotating airfoil by a correspondingly cambered a airfoil.” This 
relationship is confirmed by comparing ‘the respective aerodynamic 
-formulas.® A ‘simple an: analy tic © derivation 1 yields a parabola for ‘the dynamic 


camber. 


a 


RAPHS FOR AND Curvep OF AN 


H-Secrion 0.2 -20) 


_H-section of section ratio d/b = 0.20 (approximately the section of the original — 
(Wash.) bridge). The respective slopes are s; = + 4.00 for the 
‘straight model (corrected for drag) and s3 = +2 .00 (approximately) for the 
curved — model. Accordingly, for this section, the vertical lift resultant pro- 

duced by a a given \ B (at the leading edge) due to angular velocity i is approxi- 
‘mately one half the ver ertical lift resultant produced by an equal a due to 
vertical or to angular position (qs = s3/s: = 0. 50). 

Fig. 4 shows the static moment graphs obtained from the iniie straight and | 
curved Section mo models. -The respective slopes are = — 0.50 for the straight 
model and s, = - 00 for the curved model (corrected for drag). _Accord- 
ingly, for this. section, the aerodynamic instability torque produced by 
: leading edge) is twice e the aerodynamic instability torque produced by loud 


equal due to vertical velocity -or to ar angular position (qs = “84/82 = = 2.0), 


Ones of Curved Models and Its Application to the Study of Curvilinear Flight of Airships,’ a 
. A. Gourjienko, Technical Memorandum No. 829 (Pt. 1) and Technical mente 3 No. 830 (Pt. 2), 

‘ational Advisory Committee for Aeronautics, Washington, D. C., 1937. 
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— Since the lift is halved ‘and the torque is doubled, the eccentricity, of the 
resultant (or center of pressure) on the curved model is four times as great 
= 84/83 


“—- on the straight m model, the respective values for this section being 4 
=- — 0. (placing the resultant at the leeward edge of the. section) 


Torque Curve, 


curve, 


= 


- @. ).125. Accordingly, the ratio M/L for angular velocity is not 


re. 4—Sranc Torave | FoR AND CuRVED MODELS OF AN 
Secrion (d/b = 0.20) 


the same as for vertical velocity or angular | position. The resultant center of — 


pressure does not remain stationary (along the section) during a cycle of angular 
oscillation, nor during a cycle of coupled vertical and angular ¢ oscillations. — — 


As shown by the positive value of s. s, and the negative values of s2 and 8, 
0.20 is vertically stable and torsionally unstable. such 


See 


the section d/b = bal 
= the curved model yields increased negative sl slope « of torque graph 


> th: representing intensified ene instability. 

TABLE 1—Gaarn SLOPES FROM AND CurvED SEcTION 


— 
0 (theoretical) 6.28 1.57 | —3.14 0.25 
‘1.10 | -1.75 0.55 | 0.21 | —0.31 

(1.00 | —0.125 


—0.5 
0 0.50 


Section (d/b) 
—0.33 


0.02 (test) | 
—0.50 | 2. 
4.00 | —0.50 | 2.00 | 


0.20 (test) 
torsionally stable s sections, as by a flat plate, 


a (Fig. 6). For comparison, the thon and their peli ratios, ‘for d/b = = 
d/b = = 0.02 (flat plate tested), a and d/b =0. 20, are 


‘in Table 
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Sratic Torque GRAPHS FOR STRAIGHT AND CurvVED MODELS OF A 
Fuat Puate (d/b = 0.02) | 


cad 


6. 


TQyosanjeA 


Fiat PLATE (d/b = 0.02) 


Fie. 5.—Sraric Lirt Grapus For SrraicHt AND CurRVED MODELS OF 


Fia. 5. 
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For Oscillating Sections. —For the derivation of the aer | 


= initial angle 0 of incidence (for the rest); 
; = initial static lift at as; 
= in initial static torque at ai; 
= vertical. displacement; 
dn/dt = vertical velocity; 
= = vertical al acceleration; 
angular displacement; 


da/dt = angular velocity; 
= @a/d? = angular acceleration. 


Then’ the general aerodynamic force e equations moving 


and 
M = Mi + M(a) + + M(@) + MG) + 


.., denote the contributions to L due to @, ua snap 
MG), the e corresponding | contributions to M. All moments 
_ are taken about the midpoint of the width b. _ Rees 
wae The angular displacement a is measured from the initial angle. of incidence 
The slopes and are measured at a; (for small displacements). 


For » verticall ‘Ss mmetrical sections, as a flat. late or an H-section, in 


“horizontal wind, a; =0;L;=0;andM;=0. 


‘The ordinates of the static lift and torque. graphs are Cz, and Cy, respec: 
tively, defined by: (3 p V?b) C1; and J M = (3p b*) Cu. 
The slopes of the static lift and torque graphs are : as ie: 


slope of static | lift graph, straight model = 

= = = slope of static torque graph, straight model = 0Cy/da 

slope 0 of static lift graph, ¢ curved model = 


= slope of static torque graph, curved model = dC u/Op 


The contributions due to to angular displacement a will then be 


Cr=sa and 82 a 


rly, the contributions due to ver 


since a and i i/V) are equivalent i in aerodynamic effect. a 
- Similarly, from the slopes g3 and s, of the static lift: and torque graphs 
obtained from the curved model, the contributions due to angular velocity & 


will be 
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which B= =- 


Ce 


to a model held stationary in an curving 
or to a straight model 1 moving in an opposite pig path in a straight 
airflow. In each case the radius of curvature is R = V/d. 2 his ee 
joa radial 9 acceleration or centrifugal force involving Mass Ma of the cylinder 
of air enveloping the section. on. This radial acceler ation i is V a, andthe 
resulting inertia force is — m. Va. Since this is an extra force created - 
the curvilinear flow or motion, represented by the curved model, “it must be 
subtracted (algebraically) from the lift L given by the curved model. = It does 7 
not affect the moment M (da), since this lift correction is centered at the m idpoint ; 


(b/2). Consequently the of this correction for force dead 


Ma and M(j) =0 


The contributions due to angular acceleration & are 

= = and M(&) = 

in w n which J I, is the siiieacaiain polar moment of inertia of the enveloping cylinder _ 


7 Collecting the e component values of L and M and ‘substituting in Eqs. 


la and 1b, 


[o 


a 


in which Cy > Cs, and are are phase-correction factors the effect of time 


with the that om: pening the terms ‘ma and 
4 contain no phase-lag factors. _ Similarly — or constant contributions, 
“like L; and contain no 10 C-factors. 
Eqs. 2 are of fundamental importance. | They, are the pivotal formulas of = 
” analysis presented. Theories fo for evaluating the coefficients C1, C2, Cs, 
and C, vary; but, whatever evaluation of those coefficients is adopted, the basic 7 7 
validity of the general equations (Eqs. 2) is an independent proposition. ae 
The applicability of Eqs. 2 1 is without restriction as to the form of the sec- 
tion and without restriction as to the nature of the motion or oscillation of the 7 
Section (uniform, irregular, or harmonic). Moreover, the section may be hori- 
stationary in a steady ‘wind (as in the case of a bridge) or moving 
horizontally relative to the wind (as in the case of an airplane) ; V is the relative 


ba&/V. For this case, however, one correction is required. _ a 
a 
O augme ) Q) Mla V 
The contributions due to vertical accelerauion 7 are 
b) 
4 
mer 
M= Mitte Vd | 
| 
7 
ie 
| 


In hace case of harmonic oscillations, the coefficients Cy Ca Ca and Cy are 
functions of the velocity ratio - NO (or of k = 5 
case, the C-coefficients are of the vector form, G, representing an 
angle of lag equal to tan“! ‘G/F. As derived the values of F 
pe G ma may | be calculated and } slotted from the pressure distribution graphs for 
the section (Fi, F., Gi, and ¢ G2 f the straight model; Fs, G3, and 
‘the curved model), 


In Eqs. 2 the terms and within the brackets are dimension- 


i ynles fac tors. 


_ (a) Straight model 
-} 


Fia. 8.—PressurE DistrisuTIon GRAPHS “FOR: 
CurRvED MopELs OF A 
Fuat Puate (d/b = — 


STRAIGHT AND CuRVED MODELS OF 
H-Secrion (d/b =0.20) 
toe a thin flat plate (d/b = 0), the mass of the enveloping. cylinder of air i 


= «= For deeper sections | (H, deck, or through sections), the | 


' 


- made by simply ine including m, and I, of the circumscribin inder of air in 


sponding value is = —pr,b*, in which = 1 (d/b)?, \ w vhich 
reduces tor; = 1 for the thin plate. q 
D By aerodynamic theory, for a thin flat plate, T. = = —— pb. This be id 
: value of I, is o1 one fourth of the polar moment of inertia of the air cylinder, or fF E 
| one half of the moment of inertia about the vertical axis. For deeper sections, Js ' 
sufficient accuracy, also, = 35 Me = b4. (This last 
considers that the effective radius of is based | on and not on on 
Because the m mass Ma of thee enveloping cylinder of air is s small in comparison 
._ | with the mass m of the bridge section, the terms _ senting the effects of a 
(7 and may be omitted. ~The rrection ired) is = 
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are ‘the natural ‘frequencies i in still air, not in a vacuum. 
tional | additive corrections to m and I a 


Pressure Distribution Graghe: —Fig. 7 shows typical pressure distribu- 
tion graphs" for a straight and a curved model, respectively, for ‘the section | 
d/b = 0.20 , as obtained by manometer readings in static w ind- tunnel tests. oe 
_ Fig. 8 shows the corresponding pressure distribution graphs for a flat plate 
(d/b 0.02). Fig. is similar to the curve of pressure distribution by 


‘aerodynamic theory,? except for the « correction to zero pressure at the » leading -_ 
edge; aerodynamic theory fails to provide this correction. For the curved 


model, Fig. 8(6) differs from the (symmetrical) curve of pressure distribution 


yielded by aerodynamic theory, particularly i in the region near the leading edge, 


where aerodynamic theory i is known to be defective and inconsistent. aye 


‘The of the are designated by y: for the straight 


el. By Ap 
ACTOS the width b (or from x z = 0to 1), 
5p VB 


-2) dx 


Kutt 


These values of 81, 82, 83, and s,, refer to the stationary section, and should 


be identical with the slopes of the static lift and torque graphs (as i in Fig igs. 3 


and 4,or Figs.5and6). 


ad Effect of Phase Lag Across the Width of the Section.— —Let oi - total phase a 

~ across the width b; d= phase lag at any element dx of the width b, and ~ 

& = abscissa along the width } (taken as s unity). a Then, in the case of harmonic _ 

oscillations of offrequency N, 
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Each pr Pressure ordinate, 2 or Ya, lags an amount behind the respective 


produces it. Hence, for the oscillating 
section, _ the static ordinates 4 yi and ys oes replaced by the lagging values— 


= y (cos — isin); ; and = ys (cos $ — tsin 


: static values $1, 82, 83, and sq are replaced by th the modified lagging values: _ 


1 


| 


Cs y3e~** dz.. 


For each force or moment vector C. 8, Fs is the cosine ¢ component and 


“a Gs is the sine component. . The angle of lag of the resultant aaa or 


moment vector is tan! for the force vectors, 


ys dx (= 85) fous (= 83) 


yi (} — x) cos 1 (3 — 2) sin de 


(= 82) J n@- 


vs (2 — 2) cos = — x) sin dz 


5 


G 
re 
ond 
— 
| 
3 
— 
(5) 
and, for the = 
| 
3 
Ps 
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In Eas. 5 —s denominators (by Eqs. 3) represent 81, a 83, and s, and the 
numerators: represent Fis, Fe 82, Fs 83, | and P48, and G, 82, Gs 83, and 
Ga 8 For all cases s of FP, G, and 
For infinite | V or zero 


F =1; = =0; 1: Cc = 8. All phase difference and 
the forces and moments “acting on the section coincide in phase with their 
respective causative displacement or velocity vectors, and coincide in magni- 
tude with the values given directly by the static lift and torque graphs (for 7 


the e straight and cu curved models) without — correction for ie, lag. 


Fie. 9.—F-GRAPHS FOR FOR AN = = 0. 1.20), 


The F-graphs and G-graphs for the H-section of d/b = 0. 20, computed by 

‘E Eqs. 5 from the pressure distribution graphs of Fig. 7, are plotted i in Figs. 9° 

a nd 10; and the F-graphs ; and G-graphs for the flat plate. of d/b = 0.02, aie 

e puted by Kgs. 5 from the pressure distribution graphs of Fig. 8, are plotted in 

Figs lland120 | 
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For Sections.— —The —The general dynamic equations for a symmetrical 


M = (Mi) Kea. (6b) 


in which and are measured from the initial static and ai, Ky, 
are the spring constants or restoring-force c coefficients, and ¢; and ¢2 are the 


Values of 


No and a, are the vertical and : angular amplitudes, respectively, 


an 


S 
q 
> Fie. 10.- —G-Grarus 8 FOR an H- Secrion (d/b = 0.20) 
damping coefficients. ‘The terms Ly and M; are canceled by the identical 
"terms in the aerodynamic equations. For the uncommon case of unsy sym- 
"metrical bridge sections, , Eqs. 6 require modification by including the late lateral 7 
‘eccentricity of f the center of gravity and of tl the : stiffness axis. 
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‘in a which | vi is the phase « differer ence (angle. of lead of the ver on oscillations) i in 


_ The s structural damping (within a a cycle of oscillation) can be represented 
as a , force i in phase Ww ‘ith the velocity but of a , magnitude proportional to the 
restoring force (corresponding to Mr. -Theodorsen and I. E. Garrick" 
suggest the form: (1 + 7 and I (1 +7492) a, but this would not 
_ yield gi = 6:/m and g2 = 52/m for th the general case. _ The form used in Eqs. 8 


Combined Force Equations for Oscillating Sections. —Equating the aero- 
dynamic forces given | by Eqs. 2 with the dynamic resisting forces given by 
Eas. 8, the general force equations for an n oscillating section are obtained; thus, 


= ) +1 J2 - 


ca 


2 


‘The j i- -terms and i-terms by including m Ma in m and I. 


T1 — Css. 83) a 
11 ‘*Mechanism of Flutter: A Theoretical and Experimental Investigation of the Flutter Problem,” 

- Theodore Theodorsen and I. E. Garrick, Technical Report No. 685, National Advisory Committee = 

_ Aeronautics, ‘Washington, D. C., 1940. 
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All the ene are in dimensionless form. —_ 


_ Eqs. 10 are ee e key equati ati ions for solving the v 


namic oscillations, 


Initial Static Lift and Torque.—The initial static values of lift and torque 
(Li and M;) at zero amplitude are given by the ordinates (Cz); and (Caz): of 

: the straight model static lift and torque graphs, respectively, at the angle of 
incidence a;. This angle of incidence 1 is given by 


in “which Qe is ‘tn angle of ateiiitio of the wind and az is the static pite ho 
static angular "displacement produced by M; By 


Eqs. | 11 a positive initial torque : augments the apparent angle of 
incidence aw by the amount az to yield the corrected angle of incidence a;. 
This ag-correction, commonly overlooked, should be considered in experimental 
determinations of critical angles of incidence. ~The K relations and Kr 
relations in Eqs. 11 are written for a section model (uniform ‘, M, ”; and a). 
‘For at bridge ora scale model of a bridge, aa has sinusoidal variation and t] the 
‘mean value of ¢ Qa is reduced. 7 Consequently, section models and complete scale ; 
‘models yield different apparent values of critical angles of incidence. | Incon- - 
 Sistent results, inadequately explained by the induced angle correction for 
ratio, may be explained by the a,4-correction. 
The initial torsion ag and its “accompany ing vertical deflection “Na are func-— 
‘tien: of V and Oye ‘This explains how a change i in velocity or in inclination of 
the wind may initiate the oscillations of a Vstructure. 
‘The values of L; and M; are aerostatic forces and are balanced by the 
equilibrating deflections and restoring forces (Ki ‘na and i. a) evoked in the 


structure. Their sole influence in aerodynamic relations is in determining the 


true : angle of incidence ay from which the angular displacements @ are measured 
7 
and at which the slopes s; and are taken. 


il 

163” 
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— 

dy ol 4 
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of Incidence: ‘Numerical Example.—Given: a» = + 5°; (Cm)w - 
+ 0.070; s2 = + 1.00; Ks = 60,000 ft-lb; b = 50; and V = 100 (= 68 miles 
per hr). T hen 3 p = 30, ft-lb. By Eq. + 0.07 = + 
a= =5°+ 4° = 4+ The initial static torque, in this case, increases 
the angle of incidence from te 5° to aj + 9; half of the increase 
o due to the positive slope, ge. The static torque is thereby increased from 
= 2,100 ft- lb to M; = 4,200 00 ft-lb. (All moments are “mo are “moments per foot 


Vertical Oscillations. pure vertical ing degree of free- 


dom, or ° torsional oscillations either not initiated or damped ou wants with a@ =0 


F 81; whence 
> 


(F, i) 


in 1 which 


real terms of Eq. 12 furnish the equation: oh 
—— 


6 Gy 81, 1. This modification n of vertical frequency i increases ith V/(Ne b), 
is only a small percentage change—usually negligible. 

Vertical Instability: Numerical Examples.—For an H-section of d/b = 0. 20, 
with: b/w = = 0.268 (hence p= 0.0206), a and 81 = +4. 00, and using t the ordi- 
nates of the F,- graph presented i in Fig. 9, Eq. 13 yields. the vertical instability — 
graph plotted in Fig The highest instability range (at zero structural 
damping) extends from V/(N b) = a 08 to 2.57, with lower (“subharmonie’ _ 


— instability ranges at V/(N b) = 0.50 to 0. 68, 0.33 to 0. 40, ete. . (In the cr 


Tacoma Narrows Bridge, with d/b = 0.205, the corresponding ranges of re- 


“Rigidity and Aerodynamic of Suspension Bridges,” by D. B. B. Steinman, — 


Vol. 110, 1945, p. 450, Eqs. 29,0 
p. 558, Eq. 1580. 


Values of 6; 


inman, 


orl 

(0) 

g 

7 

t 

12) 

ee... 

0.0766 b?/w. Eg. 13 is identical with the equations 

itten for | 
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4 165 

representing substantial agreement subject to admitted inaccuracy the 

vertically “stable” 

(or 


original observations. ) ) As shown | by the graph in Fig. 13, this section is a 


(Fig. 3). 


cuc a 
The maximum value of 6 (the logarithmic increment) occurs at 
V/(N b) = 1.7 (= 6 miles per hr for the original Tacoma bridge) and amounts” 
higher of b/w. 


_" oscillations v were V/(N 6b) = 0.93 to 2. 35, 0. 53 to 0. 89, 0. 37 to 0.4 43, ¢ ete. - 


section— -that. is, one that is ‘stable | in the high velocity 
“catastrophic”) range, as predicted from the positive slope of the lift g 


raph 
to 0. 0152; all ordinates of of the 6- graph wo would be proportionately increased for 
At V/(N b) = = = 17 representing max 


tion, the graph recor ded in F ig. 10 gives G; = 
Eq. 14 ‘yields N = 1.0026 Ni, an almost negligible i increase above 
New 


m vertical instability for this sec- 
Unstable Ho 


0. ).46. ith these values, 


e the natural 


 -0.015 


Fig. 13.—Verticat Instapitiry GRAPH FOR AN 
“4 


(d/b = 0.20) 


For a flat plate of d/b = 


Fie. 14.—VerticaL GRAPH FOR 
Puate (d/b = 0.02) 


0.02, w ith s, — = + 5.25 and oe 0.268, and 
using the ordinates of the Fi-graph represented in Fig. 11, Eq. 13 yields the 
_ Vertical instability graph plotted in F ig. 14. W ith all ordinates 1 negative, this” 
shows vertical stability at all wind velocities. 
Torsional Oscillations. —For pure a angular “oscillations (single degree 
| md vertical oscillations either not initiated or damped out), with 


= 
7 

| > : 
| 
4°, | 
: 
ses 
ASe 
— 
| 
| 
: 
0 1 — 2 — 
(14) “Values of 
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= | > F 


in n which u = pb?/m = 0.0766 b?/w, and a, = 472/b% 


Eq. 16 is identical with the equations previously published for 6 in torsional | 


oscillations, except that ( B so) is now represented by the amplified 


F484 
sane! lag factor. _ This last expression is significant. The first term represents 


_ the integrated cient of angular position, a contribution previously overlooked 


or dismissed as apparently canceling out over the eyele. T The second term 
represents the integrated effect of angular velocity, but is now based on ‘static 


Ww ind-tunnel tests on curved section models. els. At At high wind velocities, although 
approaches zero as a limit, 2;and limiti—o F, Hence, for: 


sections yielding So, limiti—0 So 


it Was PI previously know that : a a $= = =4 high velocities was 


needed in order to check the known test results on torsional oscillations. The: 
— foregoing an aplication of the theory supplies this missing factor of 4, and yiel lds 


a perfect check on the recorded results of oscillating model tests. ‘- _ 
real terms of Eq. 15 furn nish the equation: (22) Frys So 


; in which 6 was the prior symbol for the 


Values of 


— 


As a rule, for negative Fas so (t tor 2 wind 
velocities, this increase in torsional frequency n may ¥ 
trophic vertical oscillations, the corresponding effect would bea small reduction | 


im frequency.) 


Torsional I I nstabilit y: Numerical Examples.— —For an H-seetion of d/ b = 0.20 
with s: = — 0.50, ss = — 1.00, b?/w = 0.238, and a; = 0.476, and using 
the F ,-values and G-values represented in n Figs. 9 and 10, respectively—Eq. | 16 


_ yields the torsional instability graph plotted in Fig. 15. This graph shows 


4 catastrophic instability ‘starting ata 1 critical velocity ra ratio of V/(N b) kde 2.50. 
1 For 


ASCE, Vol. 110, 1945, p. 460, Eqs. 51. 


Ibid., p. 559, Eq. 1600. 


18 Tbid., p. 563, Eq. 164b. 


19 “Simple Model Tests Predict. Aerodynamic Characteristics of Bridges,” i a D. B. Steinman, Civil 
ineering, January, 1947, p. 21, Eqs. 4 and 5. ieidees wis. 


’ by D. B 
20 ‘*Rigidity and Aerodynamic Stability of ‘Suspension Bridges,” by D. 
ASCE, Vol. 110, 1945, pp. 463-465. 
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Louis G. Dunn? j in 1941 from oscillating tests on a section 1 model of the : same — 


d/b and b?/w. Ww ith an indicated vertical shift of the ¢ experimental graph to 
correct for a residual positive damping that was overlooked by the investiga- 
tors,” the two graphs compared in Fig. 15 are found to be almost ee 
Without the indicated correction, Mr. Dunn’s graph showed a critical V/(N_ /b) 
of 3.38, which is now known to have been too hi high. More recent tests (1945- 7 
1946) on n oscillating models | of the Tacoma section yield values of critical 
VAN b) of 2.50 (at the California, ‘Institute of Technology in Pasadena) and 
from 2.43 to 2.56 (at the University of W ashington i in Seattle), confirming the a 
correctness of the theoretical graph i in Fig. 15 and the validity of the ~~ 


0.025 


Values of 


Oscillating tests on | 
Tacoma model 


=0.175 


_ [Corrected axis for model test on 
of residual damping(=0.052) 


15. GRAPH FOR FOR | —Tonstowat Srasiuiry GRAPH FOR A Fiat 
an H-Secrion = 0.20) PLATE rE (d/b = 


Tepresented by Eq. 16. _ The striking agreement of the slopes of the two — 


; in Fig. 15 is further confirmation of the correctness of the > theory here presented. 
W ithout the angular position term and the use of the curved model for the 
- angular velocity term in Eq. 16, the slopes and the ordinates of the co computed 
6- -graph would be only one fourth a as large. If instability graphs | can be pre- 
dicted scientifically, as in Fig. 15, without expensive and time-consuming 
oscillating model tests, the theory is justified, 
The indicated critical V/(N b) = 2.50, applied to the original 
bridge, represents | a critical wind velocity of only 13.3 miles per hr (at zero 
; damping) for ‘catastrophic torsional instability. Sue Such earlier torsional on 


tions of the : span in the. first asymmetric mode were prevented by the ‘midspan 
4 **The Failure of the Tacoma Narrows Bridge,” a report to the Hon. John M. Carmody, ars 


¢ > Federal Works Agency, Washington, D. C., March 28, 1941, by a Board of Engineers consisting of - 
1. Ammann, Theodor von K4rm4n, and Glenn B. Woodruff, Appendix VIII, Fig. 


, _ 2“ Rigidity and Aerodynamic Stability of Suspension Bridges,” by D. B. Steinman, Transactions 7 
_ ASCE, Vol. 110, 1945, p. 
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7.5 (corresponding to the wind velocity 
that destroyed the Tacoma span), the graphs i in Figs. 9 and = 0. 48 
and = =+0. 59. . W ith these values and b/w = 0.268 (or uw = 0.0206) and 
a= = = 0. 476 , Eq. 17 y ields N ers 9995 Ns, which is an almost imperceptible 
decrease below the natural frequency N». = For the same section, at V/(N b) 
= = 18. 75 (corresponding | to 100 miles per hr for the same assumed value of N), 
— F.=+ 0.91 and G, - +0. 26, and Eq. 17 yields N = 1.21 Nz, which is. a 
increase above the natural frequency. 
For a flat plate of d/b = = 0. 02, using the -graphs and -G-graphs repre- 
‘sented in Fi ‘igs. 11 and 12, with 82 = + 1.10, ss = +0. 55, b8/w = 0. .268, and 
as = 0.333, Eq. . 16 yields ‘the torsional instability graph plotted in Fig. 16. 
W ith all ordinates negative, this ne shows torsional stability at all wind 


17.—Vercror D1AGRAM FOR TORSIONAL OscILLATIONS (DRAWN FOR 8: = 4 


_ Vector Diagram for Torsional Oscillations —The terms in Eqs. 2b, 9b, and 


100, representing the aerodynamic moments in angular he al are ‘shown 


in vector relationship in Fig. 17 ¢ s-vectors = — yon =) 
moment, due to an angular position a, lags an amount: go behind a. ‘The 
(C4 84)- moment, due to an angular ve ocity « &, lags an amount ¢4 behind — a. 


k 


The | 


The componen nts the vertical (imaginary) axis yield G2 82+ Fa: 84 

opposing a. Th they represent positive damping if Gz s2 and Fs, are” 

4 positive, and negative damping (amplification) if G2 2 $2 and Fy $4 are negative— 

all in with Eq. 16. The components along the horizontal (na 
axis yield F282 — Gas opposing Therefore they represent a re- 


duction in frequency if Fes 82 is positive, and an increase in frequency if F282 is: 


_ negative, unless outweighed by the term G4 s;—all in accordance with Eq. WV. 
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Limiting AMPLITUDES AND AMPLITUDE Response CURVES 
Effect: of Curvature of the Static Lift Graph- —The equation o of the static ‘lift, 

- graph, straight ‘model, ¥ with or rigin assumed at ai, may be w ritten in the general 


form: Cr = A, +Aia a+ + Asa +: in which A, is the ordinate 
) | (Cz) and A, the slope (s:) at “iia No is the an amplitude | of vertical oscillations, 


n= No sin (w t) 
=e No COS (w t) 


The work done in a is AW = = L ndt = —F,V 
cos (w t)] dt. the expressions and a, and per forming 
Ai + + 


the integration, AW =F, b 


- For constant slope $1, AW =F (81 ao) Hence, for the curv- 


Equivalent = Ai+—A;3a%+ 34 lsat, 
Equating this expression on to ze zero ‘0 determines the limiting amplitude. 


‘equation, ¢ or if. terms in n that equation of higher - order than ¢ at may be dropped, - 
then (for zero structural damping the limiting amplitude, defined by « 


= 0, will be given by Eq. 18as 


— 


4 


wh 
"" If the slope | (F; 81) at a; is negative, representing instability, Kgs. . 19 define a 


maximum limiting amplitude, terminating the instability range. 7 If the slope 7 


~ at a; is positive, representing stability, Eqs. 19 define a minimum limiting 
Oe amplitude, beyond which potential in stability commences. For real values, 


= ine either case, by ‘Eqs. 1 19, A; :(= $1) and 4 A; must be be opposed ir in sign. . The 
de usual identification is a section of the graph curving away from the vertical 
i axis, 7 A section of the graph curving toward the vertical axis will not yield a 
e— limiting (maximum « or minimum) amplitude en 
eal) By Eg. 198, the limiting ve vertical amplitude 9m increases with V; but by 


‘Eq. 19a, the range am along the graph, defining the limiting vertical amplitude 
, is independent of V but depends only upon the geometry of the graph. 
Extending the AW-integration over a segment ¢ of the e span, using the sine- 
carve wave form, and equating to zero, the | limiting vertical amplitude at 
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Limiting Amplitudes of Torsional Oscillations. 
‘torque 1e graph, straight model, with o origin assumed at a, Qi, may , be written in ~ 
general form, Cy = + As ++, in which ‘is the 
ordinate (Cm) ‘A, the slope (82) and the equation of the. static 


torque graph, curved model, with origin at B=0 may be written similarly, 
Cu = =Bo+B 1B + Be +B B; B® + in which Bo is the ordinate (Cm) 
and B, the slope (ss) at B= 0. By the same integration as used above for 


+ - (20a) 


19 are in agreement with corresponding equations published. 


+ Asa’. + 


\ Gos. + = 0. Substituting equivalent and equivalent 

pst, the terms of fourth and higher order, and noting that 6, = 


= limiting Qo = - 4 


Eq. Bla B, = G o/k = 2, F,= 1, and with the negligible term 


Analogous to Eq. 19¢, the limiting angular amplitude at midsegment ofa a span 


wile 


21 correspond to equations previously published, with s, now added 
to 82. For s, ~ 2 $2, the calculated value of the limiting amplitude in in torsion 
is thereby multiplied by V2. This supplies a fa factor previously missing 1n in 
checking the theory against recorded oscillations new confirmation 
of the basic validity of the theory, 


ASCE, Vol. 110, 1945, p. 452, Eqs. S4aand 34.00 


= % ‘Rigidity and Aerodynamic Stability of Suspension Bridges,’ rr D. B. Steinman, Transactions, — 


Ibid, 


p. 461, Eqs. 54a and 54d. 
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For d/b = static torque curves (Fig. yield = — 0.50, 


+3.7. Hence, by Eq. the limiting segment 
40°, which checks the reported value for the Tacoma 


ey 


— Maximum limiting amplitude 
Minimum limiting amplitude 
-Minimum amplitude 
INIMUM AMPLITUDES 


Values o a 


Fic. 19.—Grapatc Mersop ror Limitina AMPLITUDES 
dded 


Psion 


This is illustrated in 


e range PQ as 


— 
| 
| 
by 
efor zero positive 
g in damping. Positive damping 4, (structural and friction) will reduce the maxi ss 
amplitudes and raise the minimum amplitudes. 
—— Graphic Methods for Limiting Amplitudes.—In any lift or torque graph, let 
— a 8q be the respective slopes at the center O and the ends, P and Q, of _ es! 
ff any range of oscillation, and let s, be the equivalent slope over thi a a 
Poo 


Papers 


affected by the of the “Then, b i by 18 and 
80 
yields as a simple graphic ¢ method (Fig. 19) for 
‘ene defined by s.=0. The two points P and Q, marking the limits — 
; _ of the range and spaced at equal horizontal ee am from 0, are given by 
= limiting amplitudes nm ¥ will then be 2 (V/w) Qm. Fig. . 19 yields a maximum 
“4 imiting amplitude, and the inverted reflection of Fig. 19 would yie a minimum 
limiting amplitude, and the inverted reflection of Fig. 19 would yield a _— 
The approximate value of the coefficient As for substitution in Eqs. 


= 


may | be similarly obtained from the ‘geometry of the lift graph: Oo 


in which a, is any convenient assumed amplitude. ‘The: same geometric rela- 


_ tion applied 1 to the two a — yields the cc corr responding coefficients As 


of limiting amplitudes, 18 and 20 should be used. 
— Use of Slope Graphs.- —Since the significant feature of. static lift and torque 
5 
graphs : appearing in the foregoing equations (Eqs. 18 to 24) is the magnitude © 
and variation of the slope of the graphs, a directly indicated geometric device — 
is the construction and application of slope graphs. — By plotting the slopes 
of any lift or torque graph as the ordinates of a new curve, the corresponding 
slope | graph is obtained (Fig. 20). _ Any range below the zero axis represents 
instability. — ~The four respective slope graphs for any section will be referred 
to as the graph, sxgraph, ss-graph, and ‘se-graph. the lift or torque graph 


graph is is a cubic equation, the slope | graph will be a par abola. 

a ie The center of oscillation O is marked on the c curve, at a; for the s;-graph 
and the ss-graph, and at = 0 for the s;- -graph and the s.-graph. 

In Fig. 20(a), for any assumed a@ (or 8), the terminal points 
and Q are marked on the curve at equal horizontal one ta (or + 8) 


from 0. The e ordinates of - F, and Q are 8, 8», and s,. The midpoint M 


of the chord PQ is on the ° vertical through O, and MO is in middle ordinate 
of the are PQ. Bq. 24, 1MO and By = 5 vilding A; or By 


positive if the chord. PQ is ove oO. Eq. = = 8 MO—yielding 
the equivalent slope for the curving range PQ. 
- From the horizontal line through O, lay off the ordinates 4 OM on the 
respective Q-verticals. The curve thus obtained is the Se- graph, showing the 

variation of equivalent slope s, with amplitude. — Since 6 is pr oportional tO Se, 


| 

Vie 

— 

' 

19 and Eqs. 21 to 24 are approximate short cuts which may be suffi- 

bi 

™ 

a 

as 

| 
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this i is also the 6- amplitude graph, showing the variation of 6 with wn 7 


‘The limiting amplitude occurs where this 6-amplitude graph crosses the zero 
horizontal axis; this is the point of demarcation between stability and insta- 


bility. It is given directly by MO = 45,. 


Fig. 20(a) shows. the construction for the general case, case, with the e center of 


oscillation O off center, at an angle of incidence 
_ For the symmetrical case, with the center of of oscillation oO at zero and the | 
slope ct curve symmetrical, the construction ‘is simplified (Fig. 20(b)). 

ordinates of the 6-amplitude graph (from the horizontal line through 0) are 
simply ‘one fourth of the corresponding o ordinates of the : slope graph. vii oe 

‘Figs. 20(a) and 20(b) are drawn for an unstable range, yielding a maximum > 


limiting amplitude. If "Figs. 20(a) are inverted, and inst insta- 


Limiting 


47 


a 


Fie. 20. —Consrruction ‘OF 5-AMPL FROM SLOPE GRAF 


bility are reversed and the intercepted amplitudes a are minimum limiting 


as a parabola; accordingly, for convenient identification, it wil be referred to 
as the 5. parabola.” | The base is the limiting amplitude (given by the fore- 
going formulas or ‘graphic constructions) and the altitude is 5,, the value of 
“4 In Fig. 20, the 5-values (varying with amplitude) are represented by the 
ordinates between the 6- parabola and the horizontal zero axis of the slope graph. 
7 If the zero axis is on the convex side of the 5- -parabola (Fig. 20(c)) the 6-ordi- 


nates i increase with amplitude and there is no limiting (maximum or minimum) 
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amplitude. In such case, a ‘6-parabola below w the axis will: represent 

of unlimited amplitude, and a above the 


procedure for plotting response curve a given 6-graph 


Values of of | a, 


15 2 4 
‘Values ves of Values of 


PING 
- 


The points of the 6- graph a are iciaiiiaaai horizontally 
give the altitudes of the 6- 5-amplitude curves. These parabolic curves have 
a common base—namely, the limiting amplitude am, , determined as in Fig. 19 

g Fig. 20; and they may be drawn as simple pasabeles or, for greater accuracy 


form (with: ordinates) i is determined as in Fig. 20. JA line 


| 
| 
: 
| 
=) 
| 
,UCTION OF VERTICAL AMPLITUDE RESPON 
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‘representing the positive damping (5s) i is drawn across the curves, horizontally 
if the: positive damping is assumed constant , (independent of amplitude). The 
points where this transversal cuts the parabolas give the respective pea 


a. Using a vertical axis at V/(N 6) = 2 (or r any convenient multiple of 7), 
each ordinate a, is multiplied geometrically by 2 


the ordinates 2 7./b of the required lmeresponse graphe 
If the positive damping 4, increases with amplitude, the graphic construction 
illustrated in Fig. 21 Tequires Reg ‘modification. A A convenient 


the base of the 6- parabola. of this parabola the graph 
gives directly the amplitude response, Moe 


52+ = 


Limiting amplitude Ob 


Fic. 23.—Grapuic Metnop FOR TonstonaL J AMPLITUDE RESPONSE 
—_ a Comparison of Figs. 21 and 22 with known amplitude re response curves 
‘indicates that the positive damping 6, is not constant, but increases (at a 
diminishing rate) with amplitude. Increase with the square root of “the 


The corresponding graphic construction for angular omeiivedes' is illustrated 

in Fig. 23. In this case, by Eq. 16, 6 has two components—the contribution - 
G2 82, here denoted by 42, and the contribution F484, denoted by that 

+ corresponding limiting amplitudes a are the constant 

the linear varying (Cm) 4) the latter given (Fig. by Bn/k. With 


| 


the 8, determines the required amplitude de response, Qo. 


| 
0 
le 
rag 
( 
» * 
| 
«wih 4 
-e : 
\ — 
\ 
| 
| ordinate of il 
| 
line — 


BRIDGE OSCILLATIONS 


| 


given by Figs. 18 to 23, ‘inclusive, are for a section model, 


or the mean amplitudes over a segment at of span. Multiply | by 4/3 t to obtain 
‘the corresponding amplitude : at midsegment. 


CouPpLeD OSCILLATIONS 
 ‘Wetieal. and Torsional Oscillations, with Damping Included.—Let: 


1 


oe 


Then Ex Eqs. 1 10 ) reduce to 


substitute | C= 


in which H; = 


GG.) — F's — + + 


— Gi Gs) 


a2 = By definition of and d Xo, an 


— 


fc 
7 
| —C, — k qa | Ca 
F and = F ld 
Apply the determinant, erms. The 
‘in which Hy = Fitkd’*'~ 
= 
Eq. 26a, 
an 
Bs 
— 
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(2) 


2 a2/82 + 4H 


in which the last term pu He is relatively negligible. By Eq. 27d, 
frequency w will normally have an intermediate value between w; and w:. 
_ ‘With a few trial values of k and interpolating, ‘the values of X, and Xo are 
found to satisfy Eq. 27a, and the coupled frequency w is then given by any of 
the remaining —_ 27. The bead (flutter) velocity i is then given — 


For zero damping, t the foregoing solution is simplified by dropping ‘the 


dy-term and the dy-term from the equations. 

Coupled Oscillations with Damping: Numerical Evample: —As shown 
Figs. 14 and 16, the flat plate (d/b = 0. 02) is a stable section against both — 
vertical an and torsional oscillations. Such sections ; have ] potential instability 


following data are data are assumed 


867 a: = 0.625 = 0.0318 s, = 5.25 5 8 = — 0.33 


1.552 a.=2.00 b= s=110 0.55 q 


Let the st structural damping. g (in both types of oscillations) be 6, = 0.05. 
‘te, n= = 92 = = 0.016, d= = = 0. 383, and = = 0.571. The damping raises Vc, 
and therefore lowerske 
Estimating k by short-cut formulas (subsequently presented), the e F-values 

and G-values in Table 2(a) are read from Figs. 11 and 12, adding — -T ri/83_ 
80 to the F,-values to F’;. The computation of and Xo by 
Eqs. 26 is recorded in Table 2(b). To 27a, X,- - 4. 77 78.62 


Hence, by interpolation i in Table 2 — 


1= 0.504 
736 —0.159° 


any and 27d, w = 141. Hence, Ve 83 
les per ‘By either | of Eqs. 29, the amplitude ratio is s 2 b) 
— 0.546 — 0.1827 = — 0.58¢ ° -, The angular oscillation dominates. ‘The 


-vertical oscillation leads by ¥ = 0. 32 = 0. 051 cycle. The center of oscilla- — 


“mil 


“tion j is shifted 0. 58 “he toward the leading im The complete solution | is self 
“hooking when th the t two vo Bags. 29 yield identical values 0 of 2 


— 
* 

tio (including the phase difference W) is given by Eqs. 25, palit 
| 

ae 

id 

1b) 

X 
| 
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0. 347 
144 Y= 


Again by E . 29, the amplitude ratio is 2/(ab) = — 0.314 = 0.081 
y qs. 29, th amplitude rz ratio is b) a 


0.32602, 
TABLE 2.—Nomerican EXAMPLE; 


From Figs. and 1 


0970 | 1.015 | 2.690 0.935 | 0.170 | 0.025 | 0.475 
0.40 | 0.953 | 1.032 | 2.633 | 0.895 — = —-0.050 | 0.600 


Eq. 2 27a a may be written in the form: 


4% in which A and B are known constants of the section; thus, in 


— 
From Egg. 26a and 30, 


n with B, Eq. 31a may be reduced 


Without re requiring t ig the final value o of k, Eqs. 31 offer ¢ an in abridged so solution, n, yield- 


ing values of X 1 and X 2 that are fairly constant over a substantial range of 
values k. and xX 2 | known, the @ is given by Eqs. 


173 — with the dy-term 
he solution is similar, but wit aluesofkare 
nero damping, the s ions. Somewhat higher v 
For the from the equations. The solution yields to 
and the ds-term dropped 475, ay A 
k=0. 
— 
0.556 
364 | 0.551 
> 


—requiring ‘only for the determination of and yielding a fairly constant 
value of w over a substantial range of values of k. 


To test the | practical accuracy of the abridged solutions given by Eggs. 3 0 
to 32, they are applied to the data of the preceding numerical example a 
A = 4. 77 and B = 78.62) and the results are given in Table 3. _ tes 


Eq. 31a or Eq. 31b, the true value o of ki is giv ven, with sharp definition, by _ 

TH DAMPING 


0.514 | -0.031 | 21.71 | 33.80 10.27 105.77 
0.504 —0.162 | 16.56 77.94 
0.499 -0.269 | 13.12 | 20.22 88. 


plying th the relation (from E Eq. 266): 


preceding numerical example, from true value is found to be 
3.—CouPLep OSCILLATIONS WITH ‘Dampine, 


SimpLiFiep APPROXIMATE 


| Eqs. 32 Eq. 3la Ea. 30 | Eqs 7 


0.504 


14.91 —0.162 | —14.94 
(0.499 


14.92 | —0.269 | —14.97 


0.349 «(0.504 —14.94 (141 =| -0. 159° 14.94 | 7. 36 


| IE | 7.65 1. 41 14.91 


TABLE 


— 


(a) Vauue BY Ea. 33 = 
‘ 


Eq. 35a 


0.970 4 085 | 1.015 | 
0.963 | 1.097 | 1.023 | 7.50 
0.953 | 1.106 1.032 


Dropping the relatively small term H s/ x Xi, Eq. 33 yields, approximately, 
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| 
| 
| 
PLE; Co 
tio, 
a) 
ria 
04 | 7.21 | «1.42 
1.41 — 
.—COoMPUTATION OF THE FREQUENCY Ratio, k 
030 | 21.71 | 33.80 | -14.91] 4.54 
2a) 0.35 16.56 | 25.92 | -14.94 | -0.11 0.346 
0.40 13:12 | 20.22 | —1497] 0.346 | 
25) 


substituting the. expressions for and Jay and dropping | the minor 
terms, Eq. 34 yields 


which, for z zero damping, reduces to 


| 

- s. » 35 offer a short cut for preliminary approximate estimates of k. Fort the 
a of the preceding examples, the results are listed in Table 4(b). With 
-* 1.41, the approximate value of k= 0.346 (instead of the true value of 


k= =0. 349) yields VY, = 83. 5 miles 3 per. hr r (instead of the. true » value of 83 


niles per hr). 


Kgs. 35 may be written i in terms of the ratio of vertical | to torsional insta- 
b=: (or. stability), 6 1/52, a8 as given by Figs. 14 and 16. - By Eqs. 13 : and 16, 


for zero structural ural damping, 


Fit+ k dy 


1432 


~The added factor indicates the rec reduction of k- seeuniiai of V.) by structural 


damping, 


of the ratio ai/a2, in these formulas, 


= section would be unstable (for at all wind 
As the difference between the frequencies (w2 — 1) | is reduced, the critical 
velocity i is lowered. Actually t the lower limit of is given b by the discriminant, 
of Eq. 26b, H?; 4H, =0, or 2X1 = This occurs with We slightly 
mating If the differences between the frequencies were reversed, ‘with 
or exceeding we, then B would be e negative : and, by Eqs. 37, k would be i imagi- 
nary, indicating that the section would be stable, with no critical velocity at 
which s self-excited coupled oscillations could develop. 
Coupled Oscillations: Comparison with the Theodorsen Equations. — For the 
A applications to coupled oscillations, Mr. Theodorsen’s equations yield 
lower: values of k and corresponding higher values of the critical velocity V e 
Part of the difference is due to the application. to a a real section (d/b = 0. 02), 


0. 
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and = If, however, the 

C, - substituted for the writer’ 
4 


instead: of an idealized case with $1 


‘Theodor sen coefficients 


regarded as recondite and formidable. 


Instability Graph for Coupled Oscillations. —The critical velocity y yielded dl 


the equations for coupled oscillations is the wind velocity of zero net amplifica- 


TABLE 5.—Noumerica Exampie, Consrruction or Insrapiuity Graph 


0.462] 1.46 | 65 | 


_ 
0.05 | 0.349] 1.41 | 83 
0.10 | 0.304} 1.39 
0.15 | 0.274 1.38 
my 
0.10 
0.254) 1.36 | 100 | 12.30 
0.287 1.35 | 116 | 13.28 - 
- 0.223) 1.34 | 122 | 14.10 
= 0; tion or damping, represented by 6-6 = , Or ete = 6,. - Consequently, the 
7 - critical velocity computed for any given value of ‘the (positive) structural 


"damping 5, is also the wind velocity that ‘produces an | equal negative damping _ 
es potential amplification ratio 6. : By calculating the critical velocities for 
different: values of 6, the instability curve for coupled | oscillations i is obtained. 
shown in the preceding section, the approximate short cut, Eq. 35a, 
|| pias. a solution very close to the exact values. This provides a quick 1 method 
bs for. computing and plotting the instability gr graph (or negative damping curve) 
forcoupled oscillations, 
Coupled Oscillations: I nstability y Graph (N umerical Example) —F or or the same 
section (d/b = 0.02) and data as in the preceding ng examples, with trial k = 0. 30 7 


for the new values, Eq. (35a yields the values of A shown in Table 5, for ‘ae 


different v: values of 5, and Eqs. 32 yield the ¢ corresponding values of w ali Plotting 


) 
and s4a check is obtained. writer's equations for couple 

oscillations (both exact and simplified) provide a convenient and completely ee 
a generalized solution for the application of any of the theories. With the other  — = 
sets of coefficients substituted, their unifying relations and identities permit a 

1e 
7 
ae 
— 
htly 
sige 
yield 
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these values, graph for coupled oscillations will be as shor shown i in 


> practically pr over the span or ane, with common nodes, SO that. the 


amplitude ratio n/a is uniform throughout the length. the vertical and 


torsional modes differ materially, a a correction is ; required i in the analysis. = 
The two force equations for for a unit won such as s Eqs. 10 or 25, soins wy be 


‘represented in abbreviated form pit. (A 2) a = and J 

+ (B) a= 0—for ‘the vertical and to torsional forces, respectively. = 
2 n/b = = + (x); = (x), i in which u; and are time 


tions, and ¢,; and are the r respective (dissimilar) wave forms (along the 


‘span or spans). _ The virtual work done by the vertical and torsional forces, 


respectively (in an element of time), | over the spans will then be oo 


and 
U1 due dx + B, U2 2 | 


The vir- 


Write the numerical rm so as_ 


tor make gig 1 dx =1; | de =1; ; and dz = U 

tual work equations then reduce to (Ad) U1 > (Az) U uz = 0; and (B:) Uw 

* (Bz) uz = 0. Substituting the « expressions represented by be As, B,, and 


from Eqs. 


Xi + id 


Proceeding as with the solution « of Eqe. 25, a duplication of Eqs. 26 to 28 
is s obtained, with U-terms added in the expressions for He and Hy, (also in 


29) as the sole 


U2 (Fs Gs + F's G2) + +La — U2)(Fi F2 — Gi Gs) 


de Fy 


— 

Coupled Oscillations in Spans with Vertical and Torsional Modes Dissimilar. 
—The preceding analysis for coupled oscillations, derived for a unit section, 
— 
— 
— for 
— 
a 
— To 

dan 

equ 
sect 
— G2 Gs) 
— “and 
— 
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and —and Eq. 29 is replaced ~ 

Ci/k — (Cs — \ 

U2 


ae 


For in which and torsional modes are similar, U=1, 

Eqs. 39 and 40 reduce to the. ‘simpler | ones previously derived. wl 

Coupled Oscillations with Dissimilar Modes: N umerical Example.—Let the 
span lengths equal J = 2,800; 1 = 1,100; and (L = 1+ 2h); and let the wave . 


l 


Vertical ‘mode—main spar span: 
Vertical -mode—side spans: gi =e (0.740 
Torsional mode—main ‘span: (sin - — 0. 493 s sin in 3 x/l) 


Torsional mode—side spans: 


To normalize these expressions: 


4 = 1; hence, _ 0.830 


Hence, U dx = 0. 932. 


= the section data of the examples = = 0. 02), 1 zero 
‘damping, the solution including the U-correction, using Eqs. 39 and 40, yields 
k= 0.43, 1.46, ve =69 miles iles per hr, and .34 —0.089 i= —0. 35 


Gener al expressions have been derived for the aerodynamic | forces on 3 any 


0s scillating section and for the o opposing ; dynamic forces. The. resulting. general 
"equations yield solutions for all problems of ver ertical, torsional, and coupled 
including analytical determination and prediction of frequencies, 

critical velocities, rate of amplification, | and limiting amplitudes. Maximum 


‘simplicity and practical applicability have been. the governing considerations. 
Numerical applications, graphic methods, and simplifying short cuts have - 


The n methods of classical seieiiaaniin theory a are limited to idealized stream- 

; line sections and are not applicable to other sections, such as practical bridge 

3 sections. In order to overcome this limitation and to secure complete gener- 

ality of application to: any form of section whatsoever, the unknown parameters 
and coefficients i in the aerodynamic equations are supplied by static wind-tunnel 7 
tests on on section models, using tilted straight models to represent the effects of : 


1 
— 
_ 
a 
| 
335 
3b) 
in 
39a) 
9) 
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vertical velocity and angular. position, and horizontal curved models to repre- 
_ sent the effects of angular velocity. . With the use of curved models and the 
added consideration of the effects of varying angular position, 1, the results of - 
prior analysis of the problem of bridge oscillations have ) been extended and m | 
amplified toward greater finality and verified validity. 
- The coefficients representing the co correction for phase - have been derived 
ins simple, general form fre from 1 the pressure distribution graphs for the straight 
and curved models, by ; simply applying to the pressure ordinates the progressive 
phase difference across the width of the section. In classical rime el 
= the corresponding coefficients (for the united case of an ideal flat plate 
or streamline | airfoil) a are derived from the theoretical discontinuity or vortex 
tes extending from the trailing ec edge of the section to infinity, so that the 
: ‘coefficients theoretically represent a funetion of the complete past history ws 
the motion of the oscillating section. 
different concepts are not identical. "The classical derivation i is not 
to nonstreamline sections, and no method of adapting it so as to fit the known 
behavior of honstreamline sections is apparent. The co coefficients derived by 
_ the simpler concept do fit the known behavior of nonstreamline sections. The 
- equations of this paper, however, are written in completely general form so as" 
_ to retain their usefulness and validity for any substitution of coefficients derived : 
by other theories, existing o or anticipated. wil 
The analysis: here submitted is believed t to. be the first generalized aero- 
dynamic theory « of oscillating sections, ns, not limited to sections of any one type. 
It is also” believed to constitute the aerodynamic theory of bridge o oscillations — 
__ brought to greater completeness and generality than any prior presentation. 
It is submitted to the profession i in order that constructive review and discus- — 
; sion may further advance the solution of the problem toward ‘definitive com- ; 
pleteness and finality. 
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| URBAN EXPRESSWAYS TO 


\e 


had such learned of paper to find that the little 
adverse comment must indicate rather general assent to the basic p premise. at 
4 ‘Mr. Steigelman presents ar an interesting but | somewhat extreme projection of 
“the basic theory advanced in the paper when he recommends p prohibiting all 
private passenger cars from the congested areas, and proposes the erection of 
radial artery perimeter parking facilities as the innermost limits of of private 
8 vehicular travel. 7 Mr. Steigelman would place these travel termini around the 
edge of the downtown z zone and suggests either walking or surface mass trans- 
portation n equipment to accomplish circulation within the critical area. . At. 
‘inst, this would appear to be a radical step to insure freeing streets from con- 
“gestion However, in view of the steady, inexorable i increase in motor car 
- registration, who can deny that some day such a drastic move may be necessary? — 
any rate, although ‘searcely hecessary now or in the foreseeable future. 4 


(since such problems have a way of solving themselves), the “possibility 1 that. 


_ Some kind of prohibition ultimately may be required should be kept in mind by 
all city planners and devisers of traffic corrective measures. 7 Incidentally, it is 
a strong argument : against a too-liberal provision of urban | expressways for a 
local traffic when such facilities are not carefully balanced with 


Intercity would automatically be valuabie lor much cross 

traffic for intracity ; any ¢ case, the 


writer agrees: heartily with Mr. Lewis’ statement: * * the inner 


on this paper has appeared in Proceedings, as follows: April, 1949, by Elmer F. Steigelman, and Harold 
M. Lewis; May, 1949, by Joseph Barnett; and September, 1949, by Ralph Wesley Jones. — 


- * Public Works Officer, U. S. Coast Guard Yard, Curtis Bay, Md oa 


__ Nore.—This paper by Curtis G. Bradfield was published in December, 1948, Proceedings. Discussion sual 
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es should’ be so designed that they can take care 
_ 2. As stated i in the Paper (under the heading, “Relief for Local Traffic”), 


in the congested zone. 1 ‘ane controlled parking, as suggested by Mr. Lewis, 


iss as much a traffic flow deterrent as unlimited ont parking and there is no 
Justification for such retarding influences. 


_ 3. The writer feels it. is not putting “‘the cart before the horse” to. ‘state 
o- the 1e cities with the most ‘eteliodion have the » most extensive mass 
transportation systems. Cities become large for economic reasons, not 
q cause of local mass transportation facilities. In such cities, business people | 
t 


end to ‘group together i in the smallest possible . cluster of centers because, oH 


to now , No o communication device ever used is as satisfactory, in the final 
nied. as s personal face-to-face contact when | consummating, or arriving at, 

agreement. Somewhat belatedly, transportation usually follows this g group- 
ing to move the workers i in and out of these trade or profession centers. The ; 


originates first and then somehow blunders along 


large investments in properties would 
in an unreasonably short time. This process" is now occurring in 


oy The only method of maintaining the value (and taxable assessment) of the — 
- downtown area is to improve the mass transportation service to that area — 


ae keeping | with the demands of the greatest n number of people le desiring to to move in 


cand out of the district 


_ No discussion of wana expressways and urban mass transportation facilities 
~ could be e considered complete without the comments of Mr. Barnett to provide 
the reaction of those who lean toward a liberal provision of expressways. Mr. ; 
_— Barnett’s statement that there are about fifteen cities with combined city 
and satellite populations of 1,000,000 or or more i is, , of course, e, accurate. - How- ie 
ever, , when it is realized that the total population of ‘these metropolises is 


Ula 


a fourth of total population of of the entire > United 


wee 


tation of any consequence. Tf: these cities, s, which are of the utmost a 


of "these fifteen do not now have off-the-st -surface mass transpor-— 


tance to the national economy, fail to meet the demand for adequate transpor- 
tation or incur huge i indebtedness for a type of facility whi which d does: not halt the 


too- rapid tendency toward | blight i in their downtown areas, disastrous us municipal 


‘ financial conditions could easily arise. 


_ As Mr. Barnett so clearly indicates, surface ‘mass transit is isina vicious cy cle 


leading only to ultimate strangulation and, seemingly, is not the answer to 
the problem > It certainly follows that, if planners make it ‘still easier and more 
attractive for the motorist to bring his car into the zone of congestion, such a 


‘policy can only a accelerate 2 this eventual strangulation. The inescapable fact 
is that the greatest number of people desiring to enter and leave e the downtown 


| 
4q 

7 

d 
d 
7 
| 
: 

= 

. p 

e 


BRADFIELD ‘ON URBAN EXPRESSWAYS» 
areas daily iain use mass transportation facilities if the service were adequate. 
These are the people must be considered, whether should be 
encouraged to use private 
whether, instead, they ‘should be discouraged from ¢ the use 2 automobiles by 
furnishing: adequate 1 mass, rapid transit, is the crux of this entire discussion. > 
is not the “occasional traveling salesman, the service mechanic, or the 
4 insurance - collector who congests cities when he brings his very necessary car 
downtown. Instead, it is the intolerably large number of people who merely 
drive to and from their place of employment and leave their cars idle al all day. 
T hey are the m motorists that must be either accommodated or discouraged. 
The writer gathers from the comments of Mr. Jones that the paper was not — 
sufficiently lucid. Atleast, there was no intention that it should be interpreted : 
as meaning all downtown vehicular traffic is unnecessary a as Mr. Jones states. 
Instead, the paper seemed to the writer to be clear i in remarking t that much uch of © 
a the downtown vehicular traffic , arising from. motorists mer ely g going to places ces of | 
employment or keeping business appointments was not really mandatory ond 
these people could be served by adequate 1 mass ss transportation, 


aad Again, the p: paper would be quoted n more re accurately if Mr. Jones’ discussion ) 


ead ‘ “many” instead of “ “most” when commenting on the number of cities: 
imperfect grid patterns of street layout. 


‘Still again, the writer is unable to find in the paper any statement which 
could be commented on, thus: “Surely, the author’s statement that expressways q 


are ‘unnecessary inside cities, but. are » needed between cities is debatable.” 
7 Instead, the writer tried to state an opinion that ¢ expressways were , needed i in 
- in for intercity and cross-the- -city traffic (and such incidental local traffic as 

could advantageously use the facilities), but that the immediate need for 


expressways for purely local shuttle traffic could not be determined accurately 
or positively until all possible methods of handling this movement volume by 


improved mass transportation facilities had been thoroughly exploited. 
Mr. Jones further states: “* * * nor should the author confine his planning | 


to circumferential expressway belts as the only geometric pattern.’ ’ No part 


7 of this paper | recommends any expressway panne. _ Incidentally, the more 


important expressways Mr. Jones describes as | “radial” are sctually “tan~ 

Fr — ” as the writer recommends, and are basically for intercity or satellite’ 


The paper should not have ¢ given the impression that the writer _— 
prohibiting the taxpayer from using his downtown streets. 7 It did recommend 
discouraging: the use of these thoroughfares unless the need was — 

mandatory—a quite different approach. — In this connection, there seems to 
be no impropriety in a policy of. discouraging taxpayers from using various 

tax- “supported facilities unless they | have an “urgent need for them. Cases i in 

point are old people’ s s homes, welfare clinics, and mental hospitals. 
Unfortunately, Mr. Jones assumed the -writer’s statement 
elevated trackage along Market Street (Philadelphia, Pa. as 
the Market, Street Elevated | (rapid transit) west of the Schuylkill R River. 
Instead, the elevated trackage (railroad) along Market Street between the 
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BRADFIELD ON URBAN ‘EXPRESSWAYS 
~ Schuylkill River and Broad | Street Station was the intended —_— of an 
"element of interference with consistent city growth. The cited 
railroad fill) is known in Philadelphia as the “Chinese Wall” and the need a 
for its removal has long been a a familiar source of complaint o on the part. val 


Other: comments ‘the: writer feels may be justified when considering Mr. 


Jones’ discussion are: 


1. Citing frequency table to contradict the statement that. 
avers expressways does not appear germane. 


- Despite Mr. Jones’ skepticism, the writer believes it is 3 self-evident that 
- net cost to the taxpayer of off-the-surface mass transportation facilities 
confined to the downtown district, even if not i entirely self-liquidating, will be 
—- less than the tremendous, nonrecapturable outlay required for an 
extensive network of urban expressways c¢ covering every part of | a f a city except 
the downtown district. Expressways cost about $5,500,000 mile and 
approximately ‘$7, 500,000 mile. Obviously, the first cost of 
_ central zone s subways will be > immeasurably less } than a city-wide network of 
expressways, a nd even the most skeptical cannot help. but admit that some _ 
_ revenue will be collected to defray part of the rapid transit cost. . On the oe 
hand, the expressways will return nothing and will entail a constantly growing 
unliquidated maintenance cost. 
_ 4, The writer agrees W with Mr. Jones ones that “Cities must be planned today 
for the forseeable needs” ‘but cannot agree that they always. have 
: been * * *.” There is little evidence that past. city planning has been any 


sen tine ae ‘unrelated series of expedient measures arising from waves of 
indignant pressure on political groups. Furthermore, the oer: doubts that 
the future definitely indicates a return to a modern n version ¢ of the ae 

- horse and buggy—the personal automobile—and the | gradual atrophy of mass 
is little place in in modern areas for 


for moving large volumes of people place to place by 


= Mr. Jones’ analogy of a hand pump compared to a modern municipal — 

iisson supply supports the writer’ s contention. The hand pump was the 
personal facility—personal transportation, al albeit either horse and buggy or 

automobile—and the municipal is the ‘mass facility—mass 
«6b. "There seems to be ttle evidence to support Mr. Jones’ contention that 
subways cannot accomm odate all transportation equipment. ‘There 
are relatively fe' few surface transit vehicles operating on Manhattan (New w York, | 

) streets and the subways, although , successfully bene 


ey Despite Mr. Jones’ ‘statement, the writer er does not recommend that: 


ke all vehicles [should be] * * * taken from the streets lie ~ did 
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‘If this i is done, it is quite 
_ that Mr. Jones’ remark—‘* * * " expressways wi will be uneconomic ; expressways 


will also be unnecessary * * *_may prove tobetru. 


if the Paper w was not entirely clear in all sncendiny a or if erroneous impressions - 
were given, a brief restatement may be justified. _ ‘The paper was, in essence, a 
plea for | moderation in current city planning—an admonition that broad 
thinking on the entire ‘subject, ne not just the motorists’ woes, might be indicated 


before arbitrary steps are taken. 


The United States is largely a nation of extremists. Its citizens will endure 
and tolerate a condition for far too long and then, suddenly, wish to alter the. 
entire situation immediately. 

“Most cities are already burdened with a heavy tax load brought about , in 

_— instances, by past inadequate expediencies or hastily planned measures : 
which failed to cope with overevaluated or underevaluated problems. == 
This traffic problem and its effect on the e economic well-being of cities in. 
the United States is a matter of ‘great concern. _ A proper | solution of the — 
‘dilemma ‘must be found. . 3 The cost of either : an expressway network or a 
composite system with some expressways and some intensively 
_ mass transportation facilities will be huge and, if the wrong decision i is made, 


there could be serious economic repercussions. 0 
"Some guidance may be f found from past experience. . There w: was a a time w! hen 
personal transportation, in the | form of a horse and buggy, threatened to drive 
the lumbering omnibus from the city streets. However, instead | of this, the 
horse car and cable car (the rapid transit of their day) were provided and the = 
seeds of the present growth. of modern downtown districts sown. Again there 
has arisen a form of personal transportation that is strangling the the existence of | 
the Jumbering buses and street cars. The writer believes it is quite ‘possible 
¥ that the wiser course would be to concentrate planning efforts on rectifying the 
. shortcomings of mass transportation | and redeveloping it to the utmost before | 

~ saddling ; the taxpayers with still more debt lar, gely for the benefit of people who 
} desire a a luxury type of transportation to their places of employment. oe - 
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EFFECT | OF 1 DRIVING SOFT CLAY | 
_ By Stuart B. AVERY, JR., AND STANLEY D. WILSON, AND Bs 
GEORGE L. FREEMAN AND JAMES D. PARSONS 4 
SruaRT B. Avery, Jr.,? AND STanLey D. 30 Assoc. M. ASCE.— 
"Factual data, on n a subject that has long been are contained 
‘this paper; and the authors are to be commended for initiating a type of in- E 


vestigation that, if conducted | systematically on various types: of 
The investigation reported was made at the site of a structure omnis 
_by point-bearing piles driven through a clay | of low plasticity and slight sen- 
sitivity, to hardpan. Under» these conditions it is difficult to conceive | 


mental effects caused by pile driving. The writers’ interpretation of the dat 
differs from that presented by the authors. 

» Effect of Disturbance on Consolidation and Strength Characteristics o of Clays.— 

a ‘It is well established that. structural | disturbance of a mass of undisturbed, 

_ Sedimentary clay changes its consolidation characteristics. T his i As illustrated 
by” Fig. 28. ‘Typical ratio- -pressure curves and corresponding shear 
_ strength-pressure curves for a sensitive clay i in undisturbed, badly distu disturbed, 

laboratory remolded states are shown in Fig. 28(a). 
Driving piles into a , fully saturated clay, without pre-excavation, causes” 
7 the total foundation volume to increase by t the volume of the piles. — ‘In addi- 


: tion, the e following e effects will result in a naclay having a a sensitive str ucture 
(1) At least a a portion of the clay is badly disturbed at at unchanged water 
content. _ This means that the effective stress in the 1e clay moves from point A 
to point B, Fig. _— In this condition the clay is only partly consolidated 


(2) The shear strength i is correspondingly reduced from point C to point D. 


Shear Stranct 


___Norg.—This paper by A. E. Cummings, G. O. Kerkhoff, and R. B. Peck was published in Desember, 
1948, Proceedings. Discussion on this paper has appeared in Proceedings, as follows: May, 1949, by | 4 
Leonardo Zeevaert, Jacob Feld, Joseph C. Fegan, and Gregory P. Tschebotarioff; June, 1949, by DE a M 
W. Skempton, James R. Schuyler, and B. A. Kantey; and September, 1949, D. 

Cons. Engr., Contuiden, 

bel » Asst. Prof. of Soil Mechanics and Foundation ae ” , Harvard d Univ., Cambridge, Mass. 
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weight would reach full at point E. However, part 

of the w Ww eight of the overlying | clay is supported by negative skin friction, be- | 

| ive een ‘the piles and the clay, the consolidation would stop at an intermediate 
point, F, corresponding to full consolidation under the decreased load. — Te 
term ‘ ‘negative skin friction” was first introduced by Karl Tecengm,’ au Hon. M. 


_ ASCE, in 1935.) 


remolded ~ 
x) 
2 
> 


Shear Strength © 


= = 


Undisturbed 


MEDIUM-SENSITIVE CLAY NONSENSITIVE CLAY 
Fig. 28.—Laporatory Tests Errect or DistuRBANCE ON CONSOLIDATION 
SrreneTH CHARACTERISTICS OF SEDIMENTARY 


4) Thes shear strengths ¢ corresponding tc to points E and F are shown graphi- a 
| cally by points G and H, respectively. ~ It should be noted that the disturbed 7 
clay, after 100% consolidation under its own weight, has a higher shear strength 


the original undisturbed clay, 
The same relationships for a nonsensitive clay are shown i in Fig. 28(b), from 


4 


which it can be seen that even complete laboratory remolding would have little 
effect on either the strength or the consolidation characteristics of such a 

material 
Pas It was from an analysis of such laboratory test curves that A. Casagrande,‘ 
M M. ASCE, in 1932, warned that.increased settlements may result from driving 


“Die Setzung der Fundierungen und ihre Wirkung auf den Oberbau,” by Karl Ingenieur, 


; | oe  4“The Structure of Clay a and Its Importance i in Foundation Engineering,” by A . Casagrande, J J ournal, ; 
Poston Soc. of Civ. Engrs. April, 1932.00 
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- friction | Bias closely s} spaced, into a sensitive clay. His conclusion applied only — 

to the case of floating, or friction, | piles, and the writers question the _ 

-plicability of the following ‘statement of the authors (under the the heading, — 


4 “The preceding information indicates that the diving of piles at the 
“site where the investigation was made did not have the detrimental effects 
which were forecast on the basis of soil tests on undisturbed and remolded 
samples and the theory of consolidation as applied by Mr. Casagrande.” 


7 ‘This statement could easily be misunderstood bya reader to mean that Mr. 
C 


asagrande had forecast detrimental effects for point-bearing piles driven 
av 

project for which the authors have conducted their investigation. The writers” 

maintain that such a “forecast” could not be made on the basis of Mr. Casa- bs 

_ grande’s paper. In order to bring out differences between the two cases, the 

- Writers | propose to discuss the four distinct aaa of pile foundations illustrated 

Case I. Point-Bearing Piles Driven Through Nonsensitive Clay— —In this 

case, the e ground surface will heave’ simultaneously with driving, the volume of C] 

this heave being oer to the volume by the piles. The clay being. 


, negligible settlement of the ground surface will 0 occur even if if ‘the clay i is aly 
disturbed by the pile | driving operation. 

Case II. Point-Bearing Piles Driven Through Sensitive Clay.—As in case I, 
the ground surface will heave simultaneously with driving, but the disturbed | 
clay will immediately start to reconsolidate under its own weight. ‘ The con- 
solidation continues until sufficient friction has developed to carry the weight 
of the clay between the : piles by negative e skin friction. ‘The — elevation of 


— 


ao dinmater of the piles, and their spa spacing. ra ‘The more sensitive 
e clay, the. greater will be the volume of water lost as a result of complete 
under the existing overburden. Although with increased q 
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the greater the pile spacing the greater will be the effective load under which 


piles that can be carried by enantiee skin friction is also reduced. — Therefore, 7 
the disturbed clay must , reconsolidate. - ‘The writers’ experience indicates that 


4 even in extreme cas cases the settlements ¢ due to. consolidation are less than te 


— 


heave is to be expected i in any case. 


In cases I and II no detrimental effects can result, unless i in case IT the soi 
or between the piles is supporting a part of the structure, or the piles are designed 


: with so low a factor of safety that they settle or fail under the additional load 7 

of the weight of some of the soil between the piles. _ 
Ww 

III Friction } Piles Driven 1 I nto onsensitive Clay. —As i in ¢ case se I, the 
7 q ‘ground surface will heave simultaneously with driving and the volume of heave . 
be equal to the volume displaced by the piles. When the load is applied, 
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: October, 1949 1949 AVERY AND WILSON ON SOFT CLAY 
. the piles and the ground surface will settle together because of the consolidation 
of the lower part of the clay stratum under the additional load of the structure. 
Case IV. Priction Piles Driven Into Sensitive Clay. —As in case II, the 
urface will heave simultaneously with driving and the disturbed clay 7 


ground su 
Ground surface after 
“Original ground surface. 


Nonsensitive 
clay 


‘Volume of heave = volume displaced by piles __ 


~Volume of heave = volume 
minus of lost during consolidation 


displaced by piles 
 Point- bearing Piles Driven Point- bearing Piles Driven. 
 Throug gh 1 Nonsensitive Oy Through Sensitive Clay 


surface immediatly after driving 


Ground surface after tent has been placed. 


4 Ground surface after driving 


’ Original ground surface 


7 
mont. 


Nonsensitive 
clay 


Volume of heave = volume 7 on of heave = Position of piles 
displaced by piles volume displaced by after load has been 
of water lost Dy I | 
‘consolidation ful lly consolida idated 


Vv. 


“Friction Piles in Clay Friction Piles in Sensitive Clay 


Fie. 29. IN Four Canes or FounpatIons 


- will begin ‘to reconsolidate under its own a weight. — However, the the skin friction 7 

will soon build 1 up » to a value w which will enable it to ¢ carry the weight of the clay 
between the piles by negative skin. friction. In case IV this condition will 
"develop only along the upper part of the pile, whereas along its lower ‘Part 
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this day load will be tr: insferred into tl the ‘adjacent clay by 1 means of positive 
skin friction. Eventually a condition of equilibrium will be established such | 
that the upper part of the clay will carry only a small part of its own weight. 
In this condition, if the | piles are loaded, the load will be transmitted by skin 
- friction into the clay along the lower part of the piles, which then must start a 
new cycle of consolidation. | If the ‘skin friction along the lower part rt of t the 
piles is not sufficient to carry the r new w load, plus the clay load, the piles would 
actually fail by slipping. Otherwise, all between the piles w will develop 
additional consolidation which will produce settlements of the piles. . In addi- | 
tion, the clay stratum underlying the points of the piles will consolidate under 7 
‘In case one should expect that the settlements would not be and 
often smaller, than if the str ucture had been supported by spread footings. 


In case IV, how ever, the settlements i in a very sensitive clay may well _— con-— 


siderably larger than if no piles had been used. — 


Experiences in Driving Piles Into Boston Blue Clay.— —The foundation of the 
John Hancock Building” in Boston, Mass., is an example of ease II, 
‘Fig. 29. i he ‘soil conditions at the s site a1 are as follows: 20 ft of sand and gravel 

‘ill, 20 ft o of ‘organic | clay, 13 ft of medium stiff clay, 80 ft of soft | clay, 10 ft of 
4 hardpan, and then bedrock. T he organic clay is known locally as Charles: 
_ River silt and, with a liquid 1 limit of about 50, has physical properties compar- 
able to those of a a very soft clay. The soft clay i is typical Boston Blue clay 
with am a medium sensitivity. 7 The foundation consisted of 14-in., 117-lb H-piles, 
— with an n average spacing of about 5 ft, driven from slightly above the surface of 
the organic. clay to bedrock. Steel H- piles, rather than a full displacement 
_ type of pile, were used 0 on this project specifically to reduce to a minimum all 
undesirable effects. If full displacement piles had been used, an adjacent 
building would first have been | vin 

to settlements that would have certainly “caused “structural damage t to the 

adjacent se section of the building. _A follower was used to drive each pile the 
4 last 20 ft. The first 20 ft of excavation preceded pile driving, whereas the 
‘remaining 25 ft followed pile driving. 

A record of the vertical movement of the surface of the foundation clay at 
- middle of the : site is shown in Fig. 30, together r with a a cumulative ‘record 
of pile driving operations. ‘The ‘parabelic curve is the estimated heave due. 


to the excavation of about 20 ft of soil. J 


Be the time the second stage of excavation began approximately 80% of 


the piles had had b been driven. _ With full displacement, this 8 corresponds toa com- z 


puted heave of about | 6 in., whereas observed heave attributed to pile 
driving amounted to only about 32 i in. ‘The main reason for this difference is 
_ believed to be the consolidation of the clay caused by ' the disturbance of the 
clay structure during | pile s driving. | | This i is indicated by the settlement A: of 
- the observation points (see Fig. 30) from August 1 to August 13 when the 
_ driving operations were nied, by which time no piles had been driven 


__ #“*The Pile Foundation for the New John Hancock Building in Boston,” by A. Casagrande, Journal, 
Boston Soc. of Civ. October, 1947, p. 297. 
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than 5 50 ft from the observation points. Whereas one should have 


the clay ‘caused by actually ‘settled. From this 
unexpected observation it is concluded that even the small displacement and 
small strains caused | by driving steel H-piles into a sensitive clay cause addi- 
tional consolidation, 
If this conclusion is correct, it would follow that little or no online pressure 
oul develop against the bottom of the | foundation slab. - To date (August, 
1080), the observed soil pressures beneath the foundation slab have — 
between zero and about 7 lb per sq in., , confirming this conclusion, one 
a The New England Mutual (NEM) Building in Boston in the vicinity of 


ona floating foundation without piles. ‘The settlements of the adjacent streets 
during: construction of the NEM building were neglibible, whereas those 
adjacent to the JH building were large. Most of these settlements occurred | 
after pile driving an and before installation internal bracing. It is believed. 
that they are chiefly due to disturbance to the organic clay, ‘caused by pile 
driving, and, to some » degree, to disturbance of the soft clay. 
_ 1: In his discussion of the Casagrande paper,‘ H. A. Mohr, M. ASCE, described 
another point-| ~bearing pile foundation in the Boston Blue ‘clay w here the clay 

suffered large 2 disturbances although : all pile holes were partly -pre-excavated 


an auger 22 in. in ¢ diameter. The foundation consisted of | pre-cast, con-- 
_ crete: piles 20 in. square, driven « on 3-ft centers. > Because of the close spacing 
it was decided to auger every other hole, « drive these piles, and then repeat 

7 7 the procedure on t the intermediate holes. — No trouble was encountered i in the 
ss @ugering ; operation until after the first group of piles had been « driven; 1; but, 
oe 7 immediately after driving this group of piles, augering of holes, even at a 
_ considerable distance from the driven piles, and removal of material from the 
Disturbance of Clay During Construction of f Chicago (Ill.) Subway. —In- 

_ teresting observations on the settlements caused by disturbance to the clay 
structure in the construction of the Chicago ‘subway have been published by 
Mr." Terzaghi,** who stated: “The t tunnel lining * * + is surrounded by a lay er 

: of completely remolded clay, several feet thick. ” This badly « disturbed clay” 


has caused substantial street settlements. 


= Overloading of Piles Resulting from Disturbance of Sensitive Clays. —The 
d additional soil load which point- -bearing | piles are called upon to carry, as a 

_ 3 result of reconsolidation, may become critical for long ‘slender piles driven 
‘through : sensitive clay. An extreme example of this situation is in Mexico 
_ City (Mexico), 1 where the foundation soil is an extremely sensitive and com-— 
pressible bentonitic clay. Point-bearing piles are often driven 1 through this 
clay to a sand layer at a depth of about 100 ft. 7 _ As a result of general lowering 
- of the piezometer surface in the clay b because of pumping from m lower strata, the 


- clay i is consolidating causing continual subsidence of the ground surface. 


“Application of Soil Mechanics in Designing Building Foundations,’ by A. Casagrande and R. 
-Fadum, Transactions, ASCE, Vol. 109, 1944, pp. 383-490. 


4Shield Tunnels of the Chicago Subway,” Karl Terzaghi_ Journal, Boston Soe. of Civ. Engrs. 


uly, 1942, pp. 168-210. 


‘the John Hancock (JH) Building, has similar soil conditions but is supported | | 
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October, 1949 AVERY AND WILSON ON SOFT CLAY 
In connection with | a project. in Males City, pile pulling tests 
were performed to determine if negative skin friction would develop to a 
“magnitude : sufficient to 0 carry the load of the « clay, and also to determine the 
rate at which it would develop. ‘The clay at the site of the test is of volcanic 
‘origin, It had a liquid limit of about 400, a plastic limit of about: 100, and a 
natural water content at or r slightly below the liquid limit (void ratio about 
8. 8.0). The unconfined compressive strength of undisturbed samples of this clay — 
was about 0.3 ton per sq ‘sq ft. 
q The test piles, in groups of four, on 3-ft centers, drove very easily, showing 
only s a slight i increase in resistance with i increasing depth. The results of pulling 
tests showed ‘that, only 2 hours after driving, the shearing force required to 
break the piles loose was at least 0.15 ton per sq ft, which is of the same order 
of magnitude as the shear strength of the undisturbed clay. This shearing 
strength i is more than ample, at normal pile diameter and sp spacing, to carry the 
entire weight 0 of the clay between the piles. _ Furthermore, the shearing 
strength o of the reconsolidated clay, after the piles we were broken was found 
to be at least one half of its maximum value, in contrast to the negligible shear 


strength of the clay when freshly remolded. 
These tests explain v why the ground surface does not subside between piles, 
even though the ‘ground surface in the immediate vicinity | of a structure is 
subsiding. Howe ever, , overloading « of the piles’ may develop « even if the | clay 


between the piles is is not in the least disturbed except for a few inches surrounding 


the piles. It is conceivable that under these conditions long, slender — 


could actually fail i in the course of time. mp sre ree 
Disturbance of Clay at Detroit (Mi ich. ) Investigation. —The clay, described 
” the authors falls somewhere between cases I and II (Fig. 29). In view of the - 
preceding discussion it is apparent that no detrimental effects could possibly 
oceur, providing the piles are sufficiently strong to withstand nd the additional 


load of the clay between the piles as as well as of the structure. — - 
& Clay with a liquid lin limit of about 30, a plastic limit of about. 19, and a natural 
water content midway between these two values, , together | with the “modest. 
overburden, indicates two possibilities: (1) The clay i is well overconsolidated, 
probably by drying; or (2) The clay has relatively little sensitivity to remolding. 

In either case one would expect that pile driving would have little effect in 7 
increasing the settlements from disturbance to the clay. _ Some settlement 
‘must have occurred due to reconsolidation of disturbed clay as s evidenced ~~ 
the fact that computed heave corresponding to full displacement i is about 4 ft, 


whereas 1 measured heave was only al about 2 ft. a 


‘Thea authors rs state, under the heading, “Results of Tests”: 7 


“In | the homogeneous oes layer constituting the lower part. of the deposit 
_ the effect of pile driving appears to have been negligible. There is no 
essential difference between the average compressive strength of the soil — 


strengths obtained from P-1, taken 1 
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to F and in During t the 
pe ay 10 months, the clay further consolidated to some point between F 
and E, resulting in an even higher shear ‘strength than that of the original 


would 


‘Under the heading, “Interpretation of Results,” the state: 


“Since the water appears to have migrated exclusively in a horizontal 
direction, the process of consolidation initiated by driving friction piles 


“into the clay would not ordinarily result in settlement.” | 


a 
ae The direction. of the flow will ¢ depend largely on the relative permeabilities | . 
vertically and horizontally, but this has no effect on the ultimate settlement. ac 


Probably the excess water resulting from disturbance to the clay structure will | a 
first drain away from the vicinity of the pile in horizontal direction (unless the | ‘ di 
pile itself permits drainage). This lateral drainage may well cause slight swell- :% 
ing in the immediately surrounding ‘mass of clay. _ However, eventually the 
excess water must be drained off to enable the grain structure to carry the 
loads, and this he-ssiyn sg must result i in settlements, unless the piles are essen- : 
Mr. Casagrande, in discussing skin friction on 


“Unless te piles themselves facilitate the drainage, the consolidation will g ° 
proceed from the surface downward, so that the shearing resistance of the “ 
clay, and therefore the skin friction along the upper portion of the piles, will of 
increase before a similar increase can take place along the lower section of Jw; 
In contrast to this, the reported results indicate that, in in this case, the J de 
veconsulidation a the badly disturbed clay surrounding the pile proceeded : at po 
the same rate along the entire length of the pile. However, since the excess { ele 
_ water must eventually find its way to the surface, or to a more pervious stratum, 7 a 
and since it is the escape | of this water that causes the settlements, the Ww riters du 
disagree with the statement in the ‘ “Conclusions” that: — fol 
“The assumption that the settlement would proceed as a result of vertical pe 
migration of the water out of the clay, which formed the basis of the Casa- P 
grande settlement forecast, is not in agreement with the observations made the 
. In this connection it should be noted that in the early nineteen thirties it J °™ 
was still generally assumed that drainage during consolidation developed maiuly 
in a vertical direction. The original, one- -dimensional theory of consolidation § 
Mr. Terzaghi includes this as one of the basic a assumptions. Experience has AS 
singe shown that lateral drainage i in ‘Inany cases more important than vertical the 
drainage. ‘That applies particularly to settlements of “buildings having” a dat 
width of # the same order of magnitude as the thickness of the clay layers. im the 
water can 1 drain 1 horizontally more readily than vertically, i it will speed up the 
. consolidation process, but it will of course have no influence on the final sae ve 
_ % ‘The Structure of Clay and Its Importance in Foundation E Engineering,’ *by A. Casagrande, J ournal, ae 
Boston Soe. of Civ. Engrs., April, 1932, p. 187. 
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| 

uld | Comparison of settlements under a relatively small loaded area with those — 
a 

the | under a a large area, with identical foundation conditions and loads, was : afforded — 

rF 


in connection with settlement analyses of the Logan Airport*® in 
inal 7 These o observations proved conclusively that, under an area comparable in size 


to an “average building, consolidation proceeded much faster than when na 
similar load was distributed over a very large area.” - 


Conclusion. .—The “writers believe that including that 
_ reported by the authors, indicates that pile driving in a sensitive, sedimentary 
_ clay will badly disturb the clay immediately surrounding the pile and that even 
wey at a considerable distance from the piles, where the strain due to deformation 


is quite small, there can be sufficient disturbance to the clay structure to ¢ cause 

In 1932 Mr. Casagrande visualized that large distortions are ‘required to. 

disturb clay structures sufficiently to cause additional settlements. Evidence = 

P- hich has accumulated since shows that even minor strains may cause e sufficient 
disturbance to some sensitive, sedimentary clays so as to result in additional 7 


a settlements. _ Thus, the effects of of disturbance are in some cases more serious 
ssen- | than ¢ one would be justified to assume on the basis of the Casagrande paper.‘ . 
Certainly, a conclusion that the effects of disturbance are of small import and 
| - perhaps entirely nonexistent, as implied by the authors, i is not warranted either 7 


on the basis of their investigation, or in the light of other available evidence. 
al Additional investigations are required to determine the limits of | the zone — 
‘disturbance. writers do not believe that: the authors’ investigation 


) ‘warrants the conclusion — *'e the region of disturbance is limited to a 
gone within na few inches of the pile Until more data are available, 
the designers would do well to make conservative assumptions regarding 
ed at » possible effects of thi this disturbance on pile { foundations i in sensitive, sedimentary _ 
atum, , “With ith this investigation the authors have pioneered i in a new, but long over- 
miters ue branch of foundation investigation. Jt 1s hoped 


due branch of foundation investigation . It is is hoped that this work will be 
followed by similar tests in more sensitive clays ¢ of greater plasticity. 7 In such a 
ry _ investigations i it would be well also to include both consolidation tests and pore 

| "pressure m measurements in the c clay before and after driving of the piles (perhaps 
made i the piezometers could be built into the sides of the piles) to furnish additional _ 


i ‘information on the amount of disturbance and on the additional consolidation a 


that is taking place i in the clay between the piles. 
dation GEORGE L. FREEMAN, ASCE, AND JAMES _ Parsons, = Assoc. M. 
ce has | ASCE —Interesting, information has been made available in this’ paper, al- 


erties! though the writers are not in agreement with the conclusions derived from the 7 
sng a § data. They believe that the situation would have been entirely different had 
piles not penetrated the uniform soft clay bed into the underlying ‘sand 


up the gravel, 


i value 36 “Analysis of Pore Pressure and Settlement Observations at Logan International Airport,” by J. P. , 

Gould, thesis presented to Harvard Univ., at Cambridge, Mass Mass., in 1949, in fulfilment of the require- 
ments for the degree of Doctor of Science. 

37 Partner, Moran, Proctor, Freeman & Mueser, New York, N. 


38 Associate ‘Partner, Moran, Proctor, Freeman & Mueser, New Y ork, [a 
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Discussions 
ASCE, provided the ‘profession with a sound explana-_ 
tion of what had been previously observed relative to the detrimental effects of 
driving piles into soft clays. His paper has done much to curb the indiscrimi- 


_ use of friction piles as a structural support in such soile. 


the bed Pony whole for the several pati considered were as slp in Table 4. 
An examination of Fig. 5 5 indicates that the 2- -in. samples suffered some 
distortion during | sampling, and a similar distortion is noted in the s samples 7 
shown in Fig. 7 Iti is assumed the in. were entirely 


IN RELATION to TIME OF Driving 


» ERAGE WATER Contant Avenacr CoMPRESSIVE STRENGTH 
PERCENTAGES) (Tons PER SQUARE 


| Adjacent to piles Midway betw een | Adjacent to r Midway between | 
(P __(P borings) piles (C borings) piles piles 


24.70 
22.8 23.8 
22.0 


Boring A. This v: is the but not giv ven in the paper. 


assumed that a all the “undistarbed"” samples tested had : ‘suffered a a degree of 


ounding the e pile, as: shown in Fig. 4, remolding 
s the soil extending away from the pile surface for at least the diameter of © 
the pile itself i. Although samples s secured adjacent to the pile, shown in Fig. | 6, 
did retain their ‘ ‘laminated structure,” to some ne extent, they are visibly onl 

disturbed and their | strengths, as from tests, are the strengths of 

4 considerably remolded soils. Hs ring A been truly = | 
disturbed, their compressive strengths would n no doubt have been substantially 

‘more than shown by the tests, and there should have been a greater ‘difference 
between the strengths of undisturbed samples and samples recovered from the 

. The degree « of disturbance m midway between the piles n must have been nearly 

od the same as at the pile periphery, because the compressive strengths 1 month 

after pile driving are the same in both cases. WwW ith Tegard to the | 


 femnineted. clay, the a authors state (see ‘ ‘General Soil Conditions”) that: 

“In their undisturbed state, the samples had an unconfined compressive — 

strength approximately eight times that of the same samples when remolded 


“in | the laboratory at their natural water content.” : 
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_ This difference is not apparent in Fig. 8, except in the upper few feet of the 
= clay. The same reduction in soil strength after 1 “month, ir in the 
‘samples adjacent to, and midway between, the piles, may be purely a coind- 
dence; and the tes test results certainly do not warrant the assumption that at the | 
disturbed soil zone is limited to within a few inches ofthe pile oe 
_ The limited scope of the investigation, and the scattered seeatin of the few 
tests | made, do ; not warrant the conclusions presented i in the paper. — Certainly 
“more than « one boring | of each type should have been made before and after 
‘driving piles. ‘It would also be interesting to know the results of unconfined 
compression tests on the samples from boring C-11. The addition of consolida- 7 
tion tests to the test program might have yielded valuable information re- 
garding t the relative strengths and degree of disturbance of the soil in different. 


The authors state that, as a result of the pile driving operations, the ond 
surface between the piles heaved about 2 ft. Piles driven i in similar soil for-— wl _ 
mations frequently cause adjacent piles to heave, : sometimes enough tor require 
-redriving. _ When the piles are point bearing and closely spaced t the resulting 
settlement. of the soil between the piles may be negligible, « even when significant 
- heaving has been noted, because of the transfer of soil weight to the pile surfaces 
through friction. The resulting arching effect between piles is sufficient to 
restrain the natural tendency for the remolded or r disturbed soil to reconsolidate. _ 
foregoing observed phenomena, however, cannot ‘properly be construed to. 
“indicate that, if the - piles were not point bearing and were in effect friction piles, 
they, or the soil between the piles, would not settle as a result of the remolding 
‘ effect and added load. the contrary, settlements should be expected 
in all probability | would occur. The ‘authors’ conclusions would be far more 
convincing had the authors determined from actual level o observations whether 
the soil did or did not settle after pile driving. 


Although the degree of restraint against settlement afforded by the piles 
depends (among other things) on the physical Properties of the clay penetrated, " 
such restraint has been observed to occur in some of the softest of soils. . At 
LaGuardia Field, in New York, , N. Y., for example, where the soft organic silt 
is ‘approximately 60 ft thick, sattlomente of the order of 5 ft have been observed 
since 1937 when about 20 ft of ash fill was placed over much of the site. . On the 
other hand, in areas s wheye structures hi had been supported on end bearing piles” 
penetrating the fill and mud and supported in the underlying s: sand, observed 
- Settlements between the piles have been. of the order of 1.5 ft for the same 
~ period. | Obviously, it ; cannot be concluded that, if the structures had been 
built: on friction piles supported in the organic silt, settlements of the lower 
order would have resulted, or a friction ‘piles would be beneficial to the 


* * the 


‘region 0 of disturbance is limited to a zone Poser Misael of the p pile and the 
main Volume | of the clay mass is apparently undisturbed.” As early as 1915, 


Karl on the i influence of pile size and shape on pon 


4q 
| 
= 
a 
{ 
— 
of FF 
‘ 
a 
ng 
of 
mn- 
nce 
nth 
ded — 


FREEMAN AND PARSONS. ‘ON SOFT CLAY 


to surrounding soil, and pointed « out that both factors are in- 
_ fluenced by the type of soil involved. In varved soils where the distortion 
of the soil may be e readily seen it is not uncommon to observe severely « distur rbed 
zones extending at least a diameter beyond the face of the piles. — 
__ In most instances the o the observed effects of driving piles i into soft. clay forma- 
tions are more pronounced than ‘stated by the authors. The disturbance 


exhibited i in Fig. 31 was a of considerable concern in one such instance 


‘Fie. 31 —Excavation ADJACENT TO§PiLtes SHowina VisvaL Evipe 
DISTURBANCE Dus to Pitz Drivina — 


and is ¢ offered as ‘as direct ‘evidence aga against the conclusions submitted by the 
authors. ~The piles shown are -point- -bearing, 15-in. closed- ended, | pipe piles, 
driven approximately 100 ft of relatively soft varved clay, of glacial 
origin, into an underlying compact sand and gravel stratum. ‘The structure 
‘supported by the | piles is shown in 2 plan i in Fig. 32 and may be described as a 
heavy. orname ital type of building. Its entire weight, with the exception ae 
the ground ficor partitions and slab, is carried by the the piles. I It was considere 
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to excavate more than their weight in soil to ewe the grade established for the 
floors. The general grades of the ground surface and floor, together with the 
- - soil f formations, are shown in Fig. 33. The piles are on a considerably wider 
_ spacing than those described by the authors; and the footings are are spaced 10 
on centers along the foundation walls. 
Piles were | driven’ from June to August, 1940, and the ground floor was: 
laid in ‘November and December, 1940. ‘Shortly before the building: 


completed and ready for 0% ccupancy, cracks. appeared in the ground floor parti- 


Floor between | footings 


ground columns and walls of 


Structure supported on piles 


Boiler 


Yj 
Yj Yy and inorganic silt 7 


Fig. 32.—Puan Suowina Location 33.—ComposiTe VERTICAL SECTION OF Som 
Bortnas AND SETTLEMENT OB- | 
SERVATION PoInts IN FLooR 


tions and slab and it was found that the ground between the pile clusters wai Ww as 
settling. ‘The settlements could be traced to no other cause than the dis- 
turbance of the soft varved clay produced by ‘driving pipe piles through this this 


As a temporary solution, underfloor grouting or “mud jacking” was done 
and holes were drilled through the floor in August, 1941, for the purpose of 

establishing : settlement observation points in the surface of the underlying: 
varved clay bed. The locations of of these observation points are shown i in Fig. 


33, and the settlements that 1 were observed : are plotted ‘against the squa square =— 


“of elapsed time in Fig. 34. Because of the straight- line relationship s so estab- 
lished it was possible to estimate, for comparison with observed settlements 
in the floor, the amount of settlement previous to the installation of the settle- 
‘ment observation points. _ ‘These settlements occurred in spite of an actual 


net decrease in on the underlying formation was for 


’ conclusions regarding 
the effect of ae any piles into soft clays are consistent neither with the data 7 


presented, nor with current aud for the following reasons: im salts 
> 


 * The conditions investigated by the authors are not analogous to the 


situation i in which friction piles have been ‘supported i in clay: — MEHL 


| 2. Had the piles in this case been supported by friction in the clay bed, 
- seitlement, of both the piles and the surrounding ground w ould have occurred po 
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ry to support either the ground floors or the partitions since they im-— 
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_ 3. Laboratory. test data submitted by the authors, and their own own state- 
ments, indicate a loss of strength in n the: clay both adjacent to and between the 
driven piless 
os Illustrations: presented by the authors indicate sc some disturbance of the | 
soil even in > “undisturbed” samples, and great | disturbance in the samples 


secured from the region a adjacent to the piles; a 


Estimated settlement before 


readings were taken 


Estimated Total Settlement in Inches 


| 


et 


— 


© - Average of points 
x - Average of points 5,8,9,11,13 
26,37,39,.41& 46 
- Average of points B,E,G,H&J- 
4 - Extensometer readings 


< 


Square Root of of Time in Months 
Inpucz By REMOLDING RESULTING FROM Pie 


_ 5. The data presented in no way | contradict the original contentions of ff Mr. 
Casagrande! regarding the 1 remolding | effect of piles driven into soft clay; and 
6. To substantiate the authors’ conclusions more borings and more labora- 
tory tests are required, including tests of the soil midway between piles. The _ 
‘simple level observations required to prove the settlement or nonsettlement — 
of the ground surface should be included in any further similar investigations. 
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DISCUSSIONS” 


MODERN CONCEPTS APPLIED TO 
CONCRETE AGGREGATE 


Hurserr K. Cook, Tuomas B. KEnnepy,! J anp Bryant Martuen. 5—_The 
paper | by Mr. Blanks, describing the « extremely interesting studies of concrete. 
aggregate w hich have been made in recent years by! the Bureau of Reclamation, 
United States Department of the Interior, is a valuable contribution to con- 


4 crete research and will be studied by workers i in this field for many years. 


The emphasis that he has placed o on the importance of. studying aggregate 
- materials as ‘components of concrete rather than as isolated unconfined rock 


a and mineral particles is perhaps the most important concept « developed from 
recent research on aggregates. The importance of this: approach has been 
recognized by the Concrete Research Division of the W aterways Experiment 
‘Station (Clinton, ‘Miss. ), and its work likewise has been directed toward de- 
veloping suitable procedures for evaluating materials i in accordance therewith. 

‘The automatic accelerated freezing- -and- thawing ‘apparatus for concrete which ’ 7 
was described in 1945 (39) (40)°* was recommended as a method of testing aggre- 


gate, among other reasons, “because | it permits materials for concrete ti to be 
tested in concrete” (39a). It is s suggested, however, that in n emphasizing this 
‘point and stating that | experience shows that good rocks do not always 


produce s superior concrete, and weak rocks do not always cause deterioration,” a 

-_ Mr. Blanks may have implied, more strongly than the data permit, that the _ 
-Teverse is frequently true. _ The reviews of the significance of the soa 
tests: of aggregates (41) (42) suggest that, in the absence of information such ; 
as tevealed by more recently developed procedures discussed by Mr. 
Blanks, the standard “tests will generally, provide a satisfactory method of 

aggregates: from bad ag aggregates. Also, as stated by Harold 


Notse.—This paper by R. F. was published i in April, 1949, Proceedings. ‘tee 
* Chf., Concrete Research Div., Waterways Experiment Station, Clinton, Miss. = | 
>. 4 Chf., Physical Test Branch, Concrete Research Div., Waterways Experiment Station, Clinton, Miss. 
aad 'Chf., Special Investigations Branch, Concrete Research Div., Waterways Experiment Station, < 


__$¢ Numerals in parentheses, thus: (39), refer to corresponding items sama in the Bibliography (see Appendix 
of the paper and at the end of the discussion in this issue). 
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niques be recognized in ‘that they ; may serve as guides for further aggregate — 


In his discussion of ‘the effects of of thermal expansivity of aggregates on con- 
‘crete, Mr. Blanks has ‘suggested that great differences in the thermal coefficient 
_ between | aggregate particles and cement paste may cause failure in concrete. 
Studies conducted at the Concrete Research Division have sug ggested that 
‘such differences, when considered as occurring between coarse aggregate 
particles and the mortar matrix of the concrete, may provide a more useful | 


index than is provided by considerations of ageregate- -cement paste ‘relation- 
ships. A large number of studies, involving g determinations of the linear coeffi- 
cient of thermal expansion of rocks used as coarse aggregates and of cement 
mortars containing different natural and manufactured sands as fine aggregate, 
and accelerated freezing and thawing tests of ‘concretes composed of such 

coarse aggregate-mortar combinations, have provided the basis for statistical 
analysis which has yielded a strong inverse correlation betw ween t thermal coefii- 
cient differences and resistance of concrete to accelerated freezing an and thawing. 
The results of this study have not been published in detail but are mentioned | 
briefly in the papers by J. ho one (43) and Bryant Mather (44). Mr. 


Gilbert has stated (43a). oe 


TE me must be ‘concluded, therefore, , that a large difference between these 
coefficients [of expansion of the coarse aggregates and mortars] should be 
avoided since the stresses set up at the interfaces of coarse aggregate par- 
ticles and mortar can be sufficient to cause loss of bond and cracking from 
the interface which will eventually — watertightness in hydraulic 


structures.” 


Mr. Blanks has stated that speci! epee? heat and heat conductivity of aggregates 
strongly affect the rate of temperature re change i in | concrete and the rate. at wi which — 
heat is transmitted through concrete. He then states that these factors are 
of limited significance for small structural members. In n an investigation 
conducted at the Concrete Research Division (45) it was found that differences 
in n thermal diffusivity ‘between coarse aggregate and mortar appeared to be 
responsible for differences in resistance of concretes to accelerated freezing 
andthawing 
ail Mr. Blanks’ discussion of alkali-a aggregate ‘reaction is ‘extremely interesting 
and valuable particularly i in view of the fact that his laboratory has, perhaps, 


done work in t this field than any His 


Reaction”): an 


“Acid to Siiealitacadl volcanic rocks are present in some of these areas and, 
= where significant quantities of these rocks have been used as aggregate in 


combination with high-alkali_ cements, deterioration caused by alkali- 
aggregate reaction i is certain to occur. 
—appears, however, to be somewhat more positive and inchuive than the data — 


“warrant. e Iti As believed that 1 more than an editorial modification would bei <7 
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“Reid to intermediate volcanic rocks, containing significant: quantities of 
- siliceous volcanic glass, are present in some of these areas and, where signifi- 
cant quantities of these rocks have been used as aggregate in combination 
with high-alkali cements, deterioration caused by alkali- a reaction 
may be expected to occur.” 
The sentence should also be qualified by a further statement, such as the fol- _ 
lowing which is taken from a previous paper (46) (46a): 


“The possibility that the expected reaction will not occur a that it 
may be forestalled—as a result of proper attention to the chemical compo- 
oe sition of the cement, the use of an inhibitor, the fact that the deleterious 
constituent is present in an amount much greater than or much less than > 
the ‘pessimum’ amount, or the favorable and exposure of the 
a structure, should always be borne in mind * * * However, in the present 
a state of our knowledge these are mere possibilities—none can be counted 


on to give a safe guarantee of freedom from trouble if reactive — 


What is sorely needed by the tiene 1g profession is a series of 
- accelerated performance tests to evaluate the behavior of proposed aggregates 
in concrete. | One such test is the accelerated fr eezing-and- thawing test pre- 
viously referred to in this discussion. Another is a nonfreezing. 
test, by which the influences of temperature and moisture-content changes - 
could be studied. ‘A number of combinations of influences such as alternate 
wetting and drying, alternate heating and cooling, and more complex conditions : 
(such as heating-wetting wi with cooling- drying and heating-drying with  cooling- 
wetting - appear to play an an important part in the deterioration of concrete. 
_ Aggregates that will produce concrete of satisfactory durability \ when exposed 
to a nonfreezing natural environment, such as is found i in the southern United 
"States, may may not produce satisfactory concrete if exposed to freezing and 


CE. Scholer, M. ASCE (47), has been | pioneering in the aeiiiimaaiid of 
nonfreezing testing techniques for concrete. The Concrete Research Division 
7 and the division concrete laboratories of the Corps of Engineers have developed 
equipment which shows great pro promise in testing the aggregates, or cement, in 
concrete, under the influences of heating-wetting with cooling-drying and of 
i The device described by Mr. Blanks for measuring the: rate of evaporation | 
from saturated specimens and the fo force associated with capillary r movement of 
water into the material has promise of becoming : a valuable tool in helping to 
- find | the reasons for the lack of durability exhibited by some materials. In its 
present state of development the device © appears to have the limitation of being 
used with only y a relatively sm: small specimen (1.5 in. in diameter and 1 in. thick). 
rs - Although the results obtained from such tests may be valid when the specimen ~ 


is representative of a material of constant pore characteristics, the interpreta- _ | 


tion of the data obtained by testing specimens from a heterogeneous gravel 
deposit: or quarry, from concrete specimens containing such materials, 
rapidly becomes more complex; and the usefulness of the data i in 1 terms: of 


- establishing a an index of durability for such materials becomes a matter of the 
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- adequacy of the sampling and the number of tests made. The question is is also 

_ raised as to the validity of application of the data obtained on a specimen of 

the size used in the test to mass concrete containing aggregate 6 in. in size. 
_ The writers agree with the statements | made by Mr. _ Blanks i in the “‘Con- 
clusions” that, at the present : state of knowledge of the evaluation of concrete. 
oe _ aggregates, t the petrographic examination may serve to establish the testing 7 

‘program necessary for each aggregate. e. The petrographic examination, ex- 
ample, would be valuable. in determining the number of pressure- -deficiency 
tests which might | be required in the the investigations of the heterogeneous 
materials referred to in the preceding pa paragraph. ee 
ait. The general acceptance of air entrainment in concrete may ev entually lead 
to a general overhauling of current concepts of aggregate suitability and 

~ method of aggregate testing; : also, the i increasing us use of manufactured sand with | 

crushed stone as aggregate for concrete necessitate modifications of 


present- day of what and test limitations, 
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AN ENGINEERING CONCEPT OF FLOW IN PIPES © 


By RALPH W. POWELL, J. C. STEVENS, GEORGE AL WHETSTONE, 
H. HIcKox, AND JAMES F. HALSEY AND WILLIAM M. 


he ad ASCE. —To present the ‘more recent knowledge 


regarding the loss of head in pipes, in a form that the practicing engineer can 
_ understand and use, is most commendable. _ The writer finds much in the paper : 
with which he is in entire agreement including : several worthwhile original ideas. _ 


The point that old experiments ; contained many ‘more values at low velocities 


than at high velocities, and that this led to the notion that f always | decreased 
with increasing velocity, is very good. Fig. 5 certainly shows the short-comings 


of the “exponential” formulas that resulted. Pe 
The explanation a as to why the “dip” at the transition from smooth to rough 
flow, which J. ‘Nikuradse found with uniform Toughness, does not occur 
in “natural” pipes, is pra praiseworthy; but, the author’ ‘idea that most. pipes 
- _ with which en; engineers deal are like the p pipes in ‘Tables 1 and 3, is questionable. 
Apparently these pipes reported in the paper had roughness particles of varying 
» but with th the variation very uniformly distributed along the length. As 
the author shows, this leads to a sudden change from smooth pipe flow to rough 
pipe flow ata a certain critical Reynolds number, which is simultaneous through- 
out the pipe. Mr. ‘Harris points | out that, if if certain reaches of the pipe are’ 
rougher than others, ‘the total drop i in head for the entire length will indicate a a 
friction, f, intermediate between the f-values of the separate reaches, and there- 


fore will | show a transition such as found by! Cc. F. Colebrook.” 36 > Assuming that 


= of John R. Freeman,”’ Past- President and Hon. M. “ASCE, and “—— 


Norr.—This paper by Charles W. Harris was published in May, Proceedings. Discussion on 
i this paper has appeared in Proceedings, as follows: September, 1949, , by E. L. Harrington. ow 


% Prof. of Mechanics, Ohio State Univ., Columbus, Ohio. 
2 


Ath. 


*‘*Turbulent Flow in Pipes, with Particular Reference to the Transition Region Between the Smooth 
and Rough Pipe Laws,” by C. F. Colebrook, Journal, Inst. C. E., Vol. 11, 1938-1939, pp. 133-156. 


_ %“Experiments Upon the Flow of Water in cee and Pipe Fittings Made at Nashua, New Hampshire, . 
June 28 to 22, 1892, by John R. Freeman, A.S.M.E., New York, N. Y., 
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others show such a seem seems, therefore, that engineers should be 


conservative and use graphs**? 129° based on M1 Mr. Colebrook’s transition— 


_—for naturally roughened | pipe, until ‘this transition region has been investi- 

"gated further. (In Eq. 7 represents the equivalent sand size roughness, 
originally called k by Mr. Nikuradse, and must be in the same units as D.) 

> —ifei in Eq. 7 is made e equal al to zero, Eq. . 7 becomes: the Nikuradse formula 

for smooth pipes. Eq. 1 is in remarkable agreement. with this formula for all 

values of R from 36, 000 to. 3,000,000, the maximum difference i in the value of 

S by the two formulas being 0.0002; but they begin to deviate 1 rapidly for higher 

—_ of R. | _ The author states that the present knowledge extends only to a 

a Reynolds number equivalent to a a velocity of 5 ft per § sec in an 8-ft pipe, which 


; 4 58° F gives a a Reynolds nu: number of 3,200,000; but Fred C. Scobey, M. ASCE, 


has reported tests on the Ontario Power Tunnel at Niagara Falls in Canada at 
20 ft per sec and in 18-ft- diameter pipe. Assuming a temperature of 55° F, this 
gives a Reynolds number of 27,500,000 and f was found to be 0.00694. _ The 
Nikuradse formula gives 0. 00703 and the Harris formula, 0. 00792. 


should not be for Reynolds : numbers. greater than 3,000,000 and, 


A “the Nikuradse formula i is s the best. ‘Fig. 15 gives a p plotting of Mr. Nikuradse’s S 
formula in the middle e part of the r range. Tf an algebraic formula is desired, 
the writer ‘suggests 


with 


e author seems troubled by the fact that the ‘Nikuradse for mula gives a 
nero value of f for an infinite value of R. ‘There i is something logically « dis- 
- turbing i in this fact, but when one e computes that it is necessary to go up to | 
a Reynolds number of about five hundred million to bring f down to 0.005, and — 
" to nearly a million million to bring it down to hal half that amount, the practical 
engineer need not be greatly worried. 


__ The derivation of f Eq. 3 is based on the e assumption that the ‘ “index” velocity, 
i which divides smooth pipe flow from rough pipe flow, is the same e for all 
diameters of pipe as long as the absolute roughness remains the same. This 
- fact, according te to the author, was shown ‘ “substantially” by the experiments 
_ listed in Table 3—that i is, the “breakaway point” came at ‘ ‘somewhat less than 

ft per sec.”” Ing Kq. 7, let R=o and divide by 2, then: 


* ‘Friction Loss in Pipes,” by Ralph W. Pow jell, , Experiment Station News, Ohio State Univ., Columbus, | 

Fluid by Hunter ter Rouse, , John Ww iley & & Sons, New ‘York, N.Y. (1946, 
yo 18 “The Flow of Water i in — Pipe,” by Fred Cc. Scobey, Bulletia No. 852, U. 8. D. re Ww ashington, . 
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«which is ‘Mr. Nikuradse’ s equation for rough pipes. es. On the o other hand, 
making « e = Oand simplifying, 
oR 0. 40 


which is the Nikuradse equation for smooth pipes. Solving Eqs: 9 simul- 
| taneously to determine the critical point where smooth th pipe : flow conge to 


rough p pipe if there i is no transition, 


Hence, with constant absolute roughness and viscosity, V; varies inversely as 7 
. square root of f. Thus, if V; is 3.50 ft per sec for the 1-in. a in n Table ? 
8, it would be 3.50 V 0.0247/0.0360 = 2. 90 ft per sec for the 4 -in. pipe and 
— 8.50 ¥ 0. 0247/0. 0183 = 4.06 ft per sec for the 4-in. pipe. © This is too much 


variation to be ignored. From this point o of 1 view Eq. 4 changes to 


fi — 0.0061 / 


J. C. Srevens,*® Past- PRESIDENT, ASCE.— “ The writer is glad to find that 
the old and tried 1 Chezy_ formula still has a champion. Int this phase of hy-. 

fe much confusion has resulted from the use of r many exponential for- 

7 ‘mulas, some of. which have become so entrenched in current thinking that it 
— will be e difficult. to eradicate them, but eradicate > them w: we must if progress is 7 
be made. Empiricism must ive way to > rationality and the complex con- 
cept to the simpler concept. 

io Of special interest to the writer were the head measurements in Table 1 
wherein differential water levels were determined to 0.00037 ft. Suc refine 
ment is not at all impossible; remote registering instruments in which water 
surfaces are recorded within on ane error of +0.006 in. (0.0005 ft) by electronic 
probes have been specified for the W aterways at V icksburg, 
 Itis ‘the twilight zo zone ne between nonturbulent flow ond fully turbulent nt flow 
that has given rise to these many exponential formulas. _ For nonturbulent 
Pasies wherein lost head varies w with the first power of the velocity and the 
velocity distribution is parabolic, the Poiseuille formula is mathematically 

. correct and needs no empirical coefficients. For fully turbulent flow, | the 
Darcy-Weisbach variation of the fundamentally ‘correct Chezy formula has 


extreme simplicity in its favor. The use of a specific friction factor r which 


Cons. Engr. (Stevens & Koon), Portland, Ore. | 
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depends solely on roughness and is independent o of pipe diameter is @ novel 

concept. If it was originated by the author he surely deserves signal credit 


for r it. . The next step follows logically; the friction factor for the particular 


pipe ii in question is: readily found from the specific friction factor : and the pipe 
diameter. It may be obtained directly from Fig. 


& many niiditad exponential formulas that are re applied quite ect 
far beyond t the limits of the experimental data on which they are based. This 

ean readily become a most inaccurate expedient, as illustrated in Fig. 5, for 
values beyond the velocity that. indexes fully turbulent flow. — This example, 
however, appears to be an extreme case. ~~ 


TABLE 4.—Losr Heap hy OnE THousaND FEET oF 


Authority head, hy, in feet,* equals: 


. Velocity, in feet persecond | 30 _30 100 


18 2g 


v2 


L Vis 
1.000343 
0.00036 


V2 


Manning 2. 87 X 0. 012? Da 
= 


Ww illiams-Hazen.. 


Percentage range of div: ersion fror from the Harris 0 to +5 0 to —16 Oto —21 


_ The writer was mildly surprised 1 to find that for { ‘the velocity re ranges | normally 
encountered i in engineering practice the many formulas give . results not too 


far apart for a given pipe. . To satisfy his curiosity in this regard he | tabulated 
the flow for eight formulas for | new cast-iron pipe, two diameters, ‘and three 


velocities. The results are give ve in Table 4. _ The percentage range « of diver- 


gence item the Harris formula (line 9, Table 4) is compared in line 10, Table 4. 
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WHETSTONE ON ON FLOW PIPES 


A. W HETSTONE, Assoc. M. ASCE. —Under the heading, “Artifi- 


October, 1 1949 


Natural Surface ces: Uniform Sand Grains,” ’ the author 
hile | “Since the tendency to dip occurs at higher velocities for the smaller 
follows that only many roughness elements are of one 
ely, With this statement and its elucidation tl the has Tesolved t the long- 
Lhis standing « conflict between the painstaking x and undeniably accurate work of 
, for | Mr. Nikuradse, and its equally obvious failure to accord with | experience | on 
uple, Bcommercial pipe, 
addition to this fundamental contribution to the basic understanding of 


flow, the author has rendered a service to engineering by his consolidation 


of known results, by his establishment of a criterion (his value of R,) for dl 

transition point to purely ‘turbulent flow, and by his | discovery | of the ‘ “specific: 

Without commenting ot on the relative merits of Eq. 1 and the more usual 

Blasius expression for this range (a close examination of the experimental data __ 

being 1 necessary to dc do this } justly), it would s¢ seem to the writer that the author's ee 

pattern of reorganization 1 and simplification of pipe friction data into three Bees 

on a (log f — log R)-graph, two of them universal and the other a horizontal - 

line with its ordinate depending on the specific roughness coefficient, deserves 
dit may be that the factor 

0.0061 or some of the values of fi in the tabulation under the heading, “Dia- — 


" grams for Smooth ar and Rough Pipes: Specific Roughness Coefficients for Pipes,” ° 
sid may be subject to change under closer scrutiny; that is a small matter. te 


- important i thing As that designers now ha have leans rational | basis for the 
engineering treatment of | pipe friction. 
3.5 
rr Hickox, M. ASCE.—In promoting the use of the Weisbach formula 
. a for head loss in fluid flow, as opposed to empirical formulas, Mr. Harris has 
ad done the engineering profession a service. The empirical formulas are based on 
oad experiments covering a limited range of flow conditions; ; and, although they » — 
la undoubtedly correct within this range, they are not correct when extrapolated 
far beyond it. For this Teason t the use of the W eisbach formula with friction 
a coefficients. depending o on the type of surface and the Reynolds number should 
0-21 be encouraged. = However, there are a number of items in the paper with which 
? 
the writer disagrees. Although some of these undoubtedly result from personal 
___ J viewpoints, others seem to be of such fundamental importance as to warrant: 
mally gf The author’ 's s concept of viscosity does not appear to be in harmony ny with the 
+ too kinetic theory. - Since viscosity is that property of a fluid which enables it to — 
lated ‘tansmit a a shearing stress, the supposition of a fluid Tesistance existing in a 
three uid of zero viscosity, or of zero shear, is simply a contradiction of terms. 
diver- The “mechanical interlocking” at the wall is merely momentum exchange 
ble 4. Associate Prof., Dept. of Civ. Eng., Texas Technological College, Lubbock, Tex. 
- Associate Director, Eng. ‘Experiment Station, Univ. of Tennessee, Knoxville, le, Tenn. 


— 
— 
— 


sequently, the absence of both i in the of 
the fluid and at the boundary. Accordingly, the he concept \ of a a zero friction 
factor at an infinite Reynolds number caused | by 2 zero viscosity is not at all 


inconsistent, 


Eq. 1 is not t supported by actual test data. Reference i is 


» which f was about at a "a Reynolds number of 1 10, 000 000. Eq yields 
f= = 0.0087. _ Reference is also made to the data presented by J. N. Bradley, 
“Assoc. M. ASCE, and S. I P. Wing,™ M. ASCE, which show values of f for the 
tunnel of the Ontario Power Company at t Niagara Falls i in Canada, a tunnel of 
ap a power plant at Livenza, Italy, and the conduit of Englewood Dam in Ohio. 
The Reynolds number for all these | data ranges between 3,000,000 and 


+40, 000,000. - The values of f are lower than those given by Eq. 1 but are in 


substantial agreement with the von K4érm4n-Prandtl smooth | pipe curve. 
tbe : Iti is difficult to understand the author’s conclusion that complete turbulence 
7 exists as soon as the condition for smooth | pipe flow i is no longer { fulfilled. To 


- arrive at this conclusion, it i is necessary to ignore the v work of J.  Nikuradse,* r 
of C.F. Colebrook*® and C. M. W hite, 36 and of the’ “writer. Test data 
by each of indicate | a ‘zone between the ‘Tegion na smooth pipe 


— flow and that of complete tt turbulence where | the friction i is not constant, but 
varies according to some law not yet fully ‘understood. a 


crt The author’s statement (under the heading, “General Principles”) ‘that 


lliaal viscosity plays ‘no part in the loss of head i in the > fully turbulent: flow 
ost frequently encountered in moderately rough pipes and channels” is 
somewhat misleading. The flow most frequently encountered is not fully 


7 developed turbulence with a constant | friction factor. Most of the test data 
presented by y the many experimenters in the field of fluid friction show a a if riction 


factor Z varying with Reynolds number . This result holds true for all the 
empirical formulas in which the exponent of of the v velocity is other than 2 Most 
of the empirical forn mulas that have been | presented are correct for the data on 


_ which they were based, and to dismiss them as inadequate because they | do not 
represent a constant factor i is at least, arbitrary. | 
these , equations fc for what they : are, explaining their divergences, and assigning 


them their proper place i in ‘hoon field of | pipe . friction | knowledge. Tr 
_ The author also states (under the heading, ‘Rough Pipes of Same Material | 

and Various 1 Diameters: Open Channels in Relation to Closed Pipes’’) that: 
nd Various radius is not used for pipes because it is an unsatisfactory corre- 


lating factor.” Throughout the paper he has used diameter. for “correlating 


ne, Va —_ Tunnel,” by G. H. Hickox, A. J. Peterka, and R. A. Elder, Trans-_ 

-——s- 88 "Laws of Fluid Flow in Rough Pipes,”’ by J. Nikuradse, English translation in The Petroleum Engi- 
neer, March, 1940, pp. 164-166; May, 1940, pp. 75, 78, 80, and 82; June, 1940, pp. 124, 127-128, and aad 
July, 1940, pp. 38, 40, and 42; ‘and August, 1940, pp. 83-84 and 87. 
$6 **Reduction of Carrying Capacity of Pipes with Age,’’ by C. F. Colebrook and C. M. White, Journal 

Inst. C. E., No. 9, 1937-1938, pp. 381-399. 
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Pipe flows; and, since the hydraulic radius of a pipe is 7 one fourth of it 


diameter, it is equally as as good a correlating factor, = © 
The author would be helpful if he a how to differentiate between a 
- series of localized | losses and a rough pipe. Ifthe writer’ s understanding of the 


‘i paper is correct, ‘a smooth | pipe may be 


come a rough pipe if there are enough _ 
“localized disturbances. The dividing line between these two concepts is 


James Hatsry Jun. ASCE, anp M. Owen.**—A simplifying 


‘unification, in Mr. Harris’ opinion, of existing methods of analyzin ng flow in 
closed conduits has been presented i in this paper. The author’s concept of this” 
problem, in which three distinct types of flow are characterized by straight — 
tine functions of friction factor, f, versus the Reynolds number, R, on a log- : 


-arithmice plot, as shown in Fig. 18, can be examined further in the light of 
_ additional data not included in the paper. OS 


10 10° 


Wee. 18. —Fricrion Factors For THREE Types or Fiow A 2-In. STEEL Pier 


~The three modes of flow in Fig. 18—those of nonturbulent: or 
1 


laminar flow, partly | turbulent fi flow, and fully turbulent flow—were first noted 
by. J. Nikuradse*® and are generally accepted by the profession. (The 
lin 


imits of the three types in Fig. 18 are up to R= 2 ,000, from R = 2 000 to 


Re: 16,200, and greater than R = 16,200, sempoctively. ) | stated by ih 


author, nonturbulent or laminar flow is subject to rational analysis, partly 

lished experimentally, full is characterized constant 
frietion factor, 

The author’ s experimental evidence Indicates that these den types of 
fi 

— from partial turbulence (that is, smooth pipe curve) to full turbulence | 

ag This: manner of transition does not : agree with the published work of C. F. 
Colebrook, 26 Hunter Rouse,‘® M. ASCE, and Lewis F. Moody," all of whom 


turbulent flow is governed by the smooth pipe law which has been well estab- 
' ow were obtained i in 2-in. and 4- -in. pipes. _ However, his data. show an abrupt 
present a gradual transition curve between partial and full turbulence hon 


Research Asst. Prof. of Fluid Mechanics, Univ. of Illinois, Urbana, Ill. 


_ Research Asst. Prof. of Theoretical and Applied Mechanics, Univ. of Illinois, Urbana, Ml. 
Forschungsheft, Verein Deutscher ingenioure, No. 361, 1933. 


Evaluation of Boundary Roughness,’’ by Hunter Rouse, 2d Hydraulics Conference, 
No. 87, Univ. of Iowa Studies in Eng., City, Iowa, 1948. 


“Friction Factors for Pipe Flow,’ "by Lewis F. s F. Moody, ‘Transactions, AS Ss. | Ve 66, 1944, Pp. 671. 
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_HALSEY AND OWEN ON FLOW IN PIPES 


tending over a of the number be 


I; ¥ 


(in which & is the wall particle | size and the other terms are as previously 
defined)—was plotted in Fig. 18 to show the generally accepted gradual 

_ transition from partial phe turbulence. . The value of k/D (the relative 

- roughness factor, in Eq. 12) was determined by substitution of the averages | 
value of f = 0278 as obtained from greater than R= = 200 in 


developed flow i in rough pipes. it will be noted that the 

Colebrook 1 function as applied to the author’ S data does not agree with the 

: straight line curve 1 used by the author until a Reynolds ake of more than 
300,000 is approached. This finding is consistent with those of other inate 


‘ gators—namely, ‘Mr. Rouse, ‘© who shows full turbulence to be developed when 


‘When a applied to the author’s data for a a 2-i in. pipe, Eq. i gives R= . 321,000 | 
as the ‘point. at which full turbulence r may be. expected to be developed. arene 
‘Since the author’ s data and curves did not agree with 0 other ‘investiga 
tions,’ 26,40,41 the writers undertook to ‘clarify the 2-i ‘in. pipe curves by : making 
tests on a 30- ft section of moderately rusty 2-in. steel pipe. This 
‘ae ty of Illinois i in n Urbana | since e 1928. The test section ws was s preceded by a 
- straight section 28 ft long and followed | by another. straight ‘section with .- 
a minimum ley length of 60 ft. _ The hi head loss was measured with a a 5-ft differential 7 
gage using m mercury, carbon tetrachloride, or a mixture of carbon tetrachloride 
and gasoline (specific gravity, 1.067) « as the gage fluid. 7 Test runs were over- 
_ lapped so that any error in reading small gage differences with a . heavy gage 
fluid was compensated for by large differences with a lighter gage fluid. The” 
pee renin measured | by weight, with a a filling time of 3 min or pene for 
“nearly all runs. _ The results of these tests are plotted as Fig. 19. oc 
An abrupt. change from partial | to full turbulence, such as the : author 
obtained, is not present i in ‘Fig. 19. A few minor deviations appear, but it is 
believed that these n may be attributed to experimental error. From values of 
R = 334,500 to R = = 465, 000, the friction factor, f, has a range of only sane 


es the average value of f = 0.0256 for full turbulence (see Table ai — 


_ 4l¢ For a copy ¢ of these tests data address the writers 9 t 121 | Talbot 
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HALSEY AND OWEN ON FLOW - PIP ES 


Using the value of f = 0. 0256, E Eqs. 12 and 13 were solved and the Cole- 
‘brook function for the writers’ data was plotted in Fig. 19. The Colebrook - 


curve falls below the writers’ curve, merging with it at a value of R of approxi- > 
mately 350,000. | The application of Eq. 14 | to the writers’ data gave a value — 
of R= - 450,000 as the point : at which full turbulence was developed. _ 


TABLE 5. —Tesr Data FROM 2- In. STANDARD ‘Steer Pree; 


‘Run M Velocity | Measured | Rey nolds Fric tion 
factor (ft per head 


| ® | | 


40.40 
39.60 02: 
38.70 0.0252 
37.42 0254 


36.00 | 430, | 0.0253 
34.60 | 418,000 | 0.0257 
33.30 412,000 | 
31:50 401,000 

29.60 | 390,000 


Ona 


368,000 
361,000 
355,000 
359,000 
357,000 
| 
348,000 


ao 


1.142, . 
1.078 | 0.1128 
0.978 0.0936 
0.854 | 0.0730 
0.744 | 0.0560 


0.660, 
0.569 
0. 495 


& 


0. 0270 
0. 
60.0: 


0.2760 | 


0.2285 


_ The fact that the Colebrook curve falls above the curve presented by Mr. 
i and below that in Fig. 19 must not be considered a defect as there is a 
considerable variation in the roughness and transition curves in each class of 

Pipe. The Colebrook function is meant to be the average of many pipes in 
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: msec R 
23.95 
3.80 
23.30 
the (22.25 50 | 4.79 | 1.800 | 91,000 | 0.0290 a We 
an 51 | 4.41 | 1.525 83,900 | 0.0288 
| 52 | 4.38 | 1.514 0.0291, 
sti- O38 | 3:88 | 1.191 | 73,500 | 0.0292 
0255 || 54 | 3.88 | 1.207 | 73/500 | 0.0296 
en 4 - 0258 || 55 | 3.77 ‘1.159 71,600 0.0300 
— 0254 || 56 | 3.37 | 0.922 64,000 299 
15] 19.59 26.75 363,000 0258 || 57 | 3.36 0.906 63,900 295 = 
14) 58 | 2.995 0.735 | 56,600 302 
(19.45 | 26.35 0256 |] 59 | 2.62 | 0.576 49,800 308 
17 19.40 | 26.25 0258 60 | 2.61 | 0.575 | 49,600 311 
19 19.00 25.10 0.389 | 40,300 318 
18.90 | 24.95 | 0-384. 40,100 316 
21 0258 || 64 | 175m | 0.272 | 33.100 00327 
334,500 0258 65 | 1.638 0.238 31,100 | 0.0328 
ga- | 337;000 0257 || | 
ing 320/000 0259 || 66 | 1.631 | 0.237 | 31,000 | 0.0328 
25 | | 313,000 0260 || 67 | 1.63 | 0.2395 | 31,000 | 0.0332 
| 291,000 | 0.0261 69 0.1724 | 25,900 0.0342 
ni- 280,000 | 0.0261 || 70 | 1219 | 0.139 | 23,100 
y a 29 244,000 | 0.0263 || 71 21,700 | — 
220,000 | 0.0265 || 72 20/500 
‘al 200,000 | 0.0268 || 74 | 16,200 
190,500 | 0.0268 || 75 14,100 
ide is.000 | 9.0200 || 
wl 173,000 | 0.0269 76 0.0457— 12,550 
rere 35 155,500 | 0.0272 || 77 10,760 
age 149,000 | 0.0273 || 79 | 0.417 7,890 
87 139,500 | 0.0273 || 80 | 0.3685 | 7,000 
The 38 130,200 | 0.0274 || 
124,700 | 0.0278 |} 81 0.01392 | 6,360 
for 40 129,000 | 0.0284 82 0.01152 | 5,840 | 
| 83 0.00951 | 5,220 | 
41 123,000 | 0.0281 || 0.00710 | 4,350 0.0500 
hor 116,500 | 0.0280 | 85 0.00368 | 2,960 | 0.0553 
1% 
1% 
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the same class. The variation of roughness and transition curve in any one 
class of pipe precludes the use of the curve in the paper since it is obvious that ; 
4 any standard of design having wide : application would necessitate the use of. 
= average curve. - _ The writers believe that the author’s method of analyzing — 


flow i in conduits represents an oversimplification | that can result in large 


— 


40 


Fie. 19.—. 


errors when predicting the flow in’pipes in n the transition re; region on between partial 
and full turbulence. _ This oversimplification is not desirable since most flow 
__ problems in engineering practice occur in the transition region bounded by 


the smooth pipe curve for partial turbulence , and aciauaea of Eq. 14 for full 
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DISCUSSIONS | 


CONTINUOUS ARCHES AND BENTS “ANALYZED > 


BY COLUMN ANALOGY 
‘Dit 


iscussion 


By LESLig J. MESZAROS, s, A. A. EREMIN, A AND > J. J. POLIVKA 


Leste J. Meszaros, Assoc. M. ASCE.—Mr. Yu has developed a m e 
direct and easily understood approach | to a problem that has been solved in = 


several ways. Among such treatments, for the record, parallel development 
by Hardy Cross," Hon. M. ASCE, the work of Mr. Cross with N. D. Mor- 
gan,” M. ASCE, and the work of Alexander Hrennikoff, “ Assoc. M. ASCE, — 


should be given. 
A. A. EREMIN,® Assoc. M. ASCE.—In this paper - the author has as developed 7 


~ the column analogy method of stress analysis in continuous rigid frames and 
arches. The advantage : of th the method is that it eliminates a solution of si simul- 


taneous equations. —Itis regrettable that, in illustrating the method, the author 
has introduced a rigid frame with straight members which may be solved v with 
less effort by various graphical or analytical methods. To illustrate this point, 7 
i: and also to check the method, the writer has solved the frame in Fig. 4 (“Ex- a 
mple 1”) by the graphical distribution of moments. ape 7 


The distance a to the y point of inflection (in Fig. 9) is determined by | 


: in which 8s is a constant that varies with the shape and elastic properties of 


‘mem 


members and can be ¢ computed analytically'* or graphically. 7 The moments 


‘Nore. —This paper by T. M. ‘Yu was published i in February, 1949, Proceedings. “Discussion o on this 
a has appeared in Proceedings, as follows: June, 1949, by Jack R. Benjamin. _ 


ASCE, Vol. 88, 1925, pp. 1196-1206. 


Frames of Reinforced Concrete,”’ by Hardy Cross and N. D. Morgan, John Wiley & 
Sons, Inc., New York, N. Y., 1932, p. 227 and p. 316 et seq., particularly p. 338. 


“Analysis of Multiple Arches,’’ by Alexander Hrennikoff, Transactions, ASCE, ‘Vol. 101, 1936, a 

Associate Bridge Engr., California State Highways, Bridge Dept., Sacramento, Calif 
\6 Transactions, ASCE, Vol. 111, 1946, p. 887. 


a. “Simplified Graphical Distribution of Moments in Rigid Frames,” by A. A. Eremin, Sacramento, 
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in the frame with the joints restrained from translation were constructed 


cross lines in Fig. 9. 
4 The basic moments i in a a member, produced by the transverse rse displacement 


byt the heaeeiaied displacement of joint D after distribution are shown in Fig. 


10. : The final tToments are obtained by superposing the moments in Figs. 9 


10, w vhich. are close to the moments computed by the author. 


“in which = ;r = =;and E = 


a 
M=13. 55 


VW F 
30 ft— —30 30 ft. 
Fig. 9.—D1acram oF Benpinac Moments 
FEET) 


1c. 10. OF Sresway 
aT Potnts D anp 


_ The advantage of the e column a analogy method applied to a single span is sin 


=. 
the simplicity « of the rule of signs. In Mr. Yu’s method, applied to continuous — 


a frames, the rule of signs is quite involved a aad therefore his method is more 


complex, than other methods i in common use. 


——— POLIVKA, 18° M. ASCE. —As an effort to simplify ‘the analysis 0 of f struc- 


: tures that ar are statically indeterminate to a high degree, this paper is a valuable” 
contribution. It constitutes an extension the column analogy method to. 
7 structures of more than three redundant elements by ‘replacing all parts adja- 

cent to a selected : span (or to a structural part) with three or less redundant 
elements—substituting members of definite elastic properties. 4 Asa ——, 
column analogy _method, this procedure is new; but the principles 
enumerated by the author in the “Synopsis” ’ (steps: (a) to (d)) have been 
E. available since 1920. ” - Furthermore, the column analogy method i As basically 
an n adaptation of the ellipse of elasticity method introduced and developed by | 
C. Culmann” in 1875 and W. Ritter* in 1888. The author limits his applica- 
tion to a general solution of continuous arches: and bents and the tracing of 


lines, whereas a more basic. approach!? will permit the solution 
4 cases in which the effect of shear an and direct ‘stress on on elastic deformations i is to 


1920, Vols. 4, 5,and6. 
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Ia the paragraph preceding Example 1, Mr. Yu states that: sodiedlad it is 
not ‘necessary (nor i is i it possible) to know even the location of ** *” the 
axis of the substitute member. That it is possible to know the location of 
this axis is demonstrated by Fig. 11, which represents the central frame — 
in | Fig. 4, 4 


37.5 


100961 


clastic area ag—that i is, the elastic weight of the substitute 


- equal to 21/208 = 0. 100961. . (The e author computed 0.1005.) It is ; concen- 
the centroid Op, along the vertical axis of the column CD 


can from Table 3A. For example, moment a y? of the 
member DD’ is 19. ‘Since ap y’p = 19. yp = V19.1/0. 1005 
= 13.785, which checks the author’ ’s computation. — Similarly, from Fig. — 
Ll, rp = = (80/21)? X 13 = 188.66214, checking with the author’s value of | 
19. 110. 1005 = 190.04975. 
In introducing Fig. 1, the author states: “The extension Lof the column — 
analogy] is not applicable to frames with closed panels 12 
‘demonstrates. the analysis a two- “story frame with rectangular panels 
ABEGDC, one of which is s closed. J A horizontal force’ H=1 acts at | point E, 
as shown. The characteristics of the substitute members BB’ and DD’ 
ean be computed from data in Table 2C. In this | case also, a more accu- 
tate result is obtained by ratios. For. example, referring to Fig. 10, 
tz = ap = (17 — 9)/104 = (30 — 22)/104 = 1/13 = 0.0769231. _ The author 
reports: 0.163 — 0.087 = 0.076. Since 17 — 9 = 30 — 22, the elastic cen- 
toids of the suhetiiuhe members are at the wsiiipeinte of columns AB and CD. 
tz of the elastic area is equal to az rp. The 


displacement along BD is equal to 1,200 1,600 _ = 41.025641, 


by 
ent 
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15) 
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TABLE 6.—Eastic ‘Loap AND ForMULA FOR mi, 


Corrected to contreld.. 


1222 POLIVKA ON ARCHES AND 

i= 

a 1 co" 12— i, 

pe |1.00| -15| 0 | -15 o | 2 | = 
EG | 0.50 +10] 0 | +5 | 37.5 | a 
ts of internal stresses (forces): mi 


a8 compared with 2 X 20.5 = 41 in Table 2C. WwW ith the knowledge ol 
ixw x ap x and (d 10) 10 = r’p, owe can compute: d = 160/3; 
= 1,300/3; a and i, = 100/3. The complete analysis of the two- -story frame 
given in in Tables 6 and 7 . Table. 6 is similar to Table3A. 
TABLE 7.—Moments M = m, EXAMPLE 


iat 


~ X13 AB= ~ 181 X13 181 


-POLIVKA ON ARCHES AND BENTS 1283 
7 


of a continuous frame with closed panels is the Vierendeel 
open web truss (see Fig. 13). _ The characteristics of the substitute members 


are 2p = 64/3; ts, the second moment, is 4/3; and, | a, the elastic area concen-— 


on Uke 


wien 


8400 _, 
= 442.10 


47.37 7400 _ 389.47 
32 ft = 


‘Fie. 


ory 


trated at the — of the chords, is 1/16. The analysis’ is similar to that 
in Fig. 12, and the results in Fig. 13 check the solution of this problem pre- 
sented by A. ‘A. Anirikian, M. . ASCE, in 1942. = 
7 Corrections for Transactions: In February, 1949, Proceedings, on pages ges 173 
and 174: In Table 2B, center diagram, reverse the direction of the vertical and 
horizontal l force arrows and change the corresponding minus sign to plus. _ For 
the vertical force Y =1 in Table 2B, under the columns for P, My, nd Mz, 
Teverse the signs; and i in line 6, under Y = 1, reverse the signs of each quantity 
except 0. In Table 2C(a), under tz = 1, change the minus signs to plus sin 


ae ‘Analysis of Rigid Frames (An Application of Slope Deflection),’’ by A. Amirikian, U. 8S. Govt. 
Printing Office, Washington, | D. C., 1942, 
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7 APPLIED COLUMN THE 


Discussion 


By JAcop KAROL 


Karo,” M. ASCE. —This comprehensive paper by Mr. Shanley is 
an admirable ‘presentation of ‘the theory underlying the present methods of 
stress analysis of compression members i in aircraft design. The mathematical 

-4 development of the nondimensional stress-strain diagrams and column — 

: may | seem somewhat ‘ ‘impractical” to the civil engineer who uses allowable 
— stresses in the design of columns of steel, concrete, or timber that : are fixed by 
the design specifications at values well within the elastic range. In. aircraft 
- 4 structures, however, the design stresses are in the inelastic range. Since air- 
planes. are generally static tested for full design loads, | it was “necessary for the 
sad stress analyst to develop a , design theory w which, in the ideal case, would result — 
a structure that failed at slightly beyond the design load. , the 

- detailed study, by the author, of the several phases of applied column theory 


vas 


is undeniably ‘practical.’ 


The writer will confine his discussion to Sections 6 and 7. — 
The method of correcting test data t to o “standard,” ” as indicated i in. ‘Fig. 12, 
is correct only for columns that fail by primary instability, since the » el 
; related column curves are derived from the tangent modulus. The horizontal _ 
scale marked “not used” designated as “effective i in the original Lock- 
heed report (see Fig. 12). This value is actually the important variable in 
general case, as will be evident subsequently. 
It is supposed that series of compression s of sheet stiffener combina- 
tions have been made. _ Since the sheet and stiffener may each be subject to 
— local buckling, the column curve plotted from the test data will not be affinely 
related to the tangent modulus column curve of the specimen material. — lt 
may be noted in passing that determination of ;* and L/p for a sheet a 
a combination i is quite a problem in itself, but it is assumed that they have been 
determined somehow. —Itis desirable to plot the test data on a graph, which 


also has the Euler curve plotted thereon, to determine e the point ¢ of intersection. 


Nore.—This by F. R. was published in June, 1949, Discussion on 


_ ‘Paper has appeared in Proceedings, as follows: September, 1949, by William R. Osgood, and J. V. du Plessis. 
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KAROL ON COLUMNS 
The procedure for correcting the data for enn stress f, to standard is as 
{ollo 
me : fe be the : yield strength f for standard ma material; K f. be the yield 
strength for specimen material; f, be the stress for (L/p):, as reduced to stand-— 
ard; fe be the stress for (L/p):, read” from the curve; and fs be the stress for 
(Lip)s = (L/p), INK, read from the curve. Then, for the ratio 


The foregoing is a verbal ‘statement of the. mathematios 
7 expressing the affine” relationship between the test column cu curve and the 
derived” column curve for the standard material. interested reader 
readily prove this for himself by assuming ¢ that the test column curve is mathe- 
matically expressible in terms of f, and K f., as indicated elsewhere,'® and by 
calculating the stresses fi, fe, and f;. Hence, ‘in the general case, correction “4 


me 2) 


§ standard is a function of the shape of the test column curve, the yield strength 


yis "tio, amd Lip, 
? ie The author’s extended treatment of the effects of eccentricity in n pin-ended_ a 
columns is undoubtedly an indication of the importance of the subject. After” 
rejecting the secant formula inapplicable to materials not having pro- 
able nounced yield point, he ‘suggests t using interaction curves to solve the problem, 7 
1 by and presents in Fig. 15 such a series of curves for axial loading an id bending - 
raft “due to third- -point loading. He assumes correctly that, for an initial moment. - 

air- due to constant eccentricity, the : sag would be greater. He further states that - 
the | ‘the shape of these curves may be primarily : a function of the kind of loading — 7 
oe “involved, which is also. true. Hence, a number of sets of interaction curves» f 
the would be necessary in design, unless there is some other approach | to the 
eory Problem, 
Be After following the author’s lead and verifying the statements referred to 
= in the preceding paragraph, it occurred to the writer that the problem could — 7 
12, be generalized so that one interaction curve would apply for various loadings, — . 
inely t perhaps for any shape of cross section and any type of stress-strain diagram. - 
yntal The reasoning was extremely simple. For any column subject to combined 
‘ock- @ bending and compression, the maximum : stress on the section is the sum of the 
nthe ™ ®Xial stress and the total bending stress. The total bending stress is itself the 
«sum of the primary bending stress and the secondary bending stress due to 
tii, | latera al deflection. _ By us using the total bending stress as a parameter, the type - 
ct to of loading is automatically cared for. If the total stress on the section is” 
indy known n, then, by deducting the axial total bending stress is obtained. 
. It ~The primary bending stress can then be determined by using | the available — 
been The problem resolved itself into finding an expression which would ‘equal 
which the critical column stress for pure compression | and the modulus of rupture for : 


pure bending . The expression was originally determined to fit the test data 
in | Fig. ‘14 for the entire range of R’, and L/ p=3 30. It was later found that 7 
the. same expression checked quite well any point taken at random from 
Fig. 14. _ This might be considered as an indication of its ‘possible. e generality. 


p* % “Strength of Metal Aircraft Elements,’”’ ANC-6, Army-Navy- Civil Committee | on Aircraft Design 


nteria, U. S. Govt. Printing Office, Washington, D. C., 1942. 
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expression: 
foe = for + (m 


‘The total stress in combined compression. and is. given | by the 


in which 


‘Dividing both sides of Eq. 49 by fn 


+ (m 1) (0.839 tan 8)... 


mR if m — 1) (0.839 tan @).. 


Eq. 
for prohlen— 


—can be determined for any value of 


The remaining problem is to determine the primary bending stress f, from 


“the total bending stress for. For any distributed loading (including: the test 


-third-p oint loading), as well as for initial curvature, the stress is given accu- 


(56a) 


and, for con: constant —— vat axial load, 


f= 140.2724 


(68) 

56 ‘are in ‘Fig. 20(a) Thus, the problem of combined 
sa load and bending due to eccentricity or transverse loading i is solved from 


Moo 


anc -< 
+ 
Tor various values Of m. since m 1S Known 
= = 
“aig 
(56) 
all 
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1228 -KAROL oN 
_ Two examples of the use of Fig. 20(a) applied to the data in Fig. 14 follow: 


(1) Given: Bea 30; and, in kips per square inch, foo = 120, * hed 60, 
£ = 40, and E= 29, 200. _ The problem i is to compute fy. From the assumed 
data, m - = 120/60 = 25 7s 320 kips per sq in; @ 40/320 = 0. 125; and 
wr = 40/60 = 0.66 667. _ From F Fig. 20(a), Ru = 0. 38 ; hence, fue = = 0. 38 x 120, 
7 = 45.6 kips | per si ‘sq in. From Fig. 20(a), folfor = 0.875; hence, s =0. 875 
x: 45.6 = 39.9 kips per sq in. ‘From ‘Fig. 14, test fo = 40 kips per sq ing 

(2) Given = 120; and, i in kips per: square inch—fro= 120, Sor= fe = 20, 
a and f= = 50. Itis required to compute fy. - From the given data, m = 120/20 
‘= 6; folfeo = 50/120 = 0. Guess = 0.498 =a. From Fig. 


folfor = 0.502; hence, Ry = 0.83. From Fig. 20(a), R’. = 0.498, ‘which checks | 


the guess. Hence, | ety 20 x 0.498 = = 9.96 kips per sq in. - From Fig. 14, 


“test fe= 10 Kips per 
«dt is asa ‘that es addition of stresses in the inelastic ran range is a 
convenient fiction; ; but, in the foregoing treatment, Eq. 49 h has the merit of 

being a close approximation to the experimental truth and of indicating | the 

- important va variables in a problem not solved hitherto in the general case. Since 


from pure compression to pure be bending, it seems logical enough to be generally 
valid. Furthermore, if is considered represent not only the ultimate 
modulus of rupture, but also a limiting bending stress i in general, the usefulness 
of Fig. 20(a) is extended . Thus, in in aircraft design, it can be used for limit 
loads as well as for ultimate loads. _ The w writer | was able to check the test. data 
by J. A. Van . den Broek,*1 M. ASCE, on 1-in. mild steel round bars for both 
initial yielding and Miinae capacity by using the yield point as the limiting 
bending stress in the first case ¢ and 1. 7 times the yield point as as | the modulus of of 


rupture in the second case. 


a. The column theory, ‘from wi Ww hich the secant formula is derived, is se ie 7 


_— for comparison. "The total stress is given by the expression: 


Dividing through both sides of Eq. 59 by m fer, 


— 


as stress ratios, Eq. 60 becomes 


= 


Eq. 61 is a a family of straight lines, as represented in Fig. 20(0) for various” 
values of m » | t is obviously correct only for a value of m =1 This is true in ; 


_ the short column range for structural steel, in which the > limiting stress for 


. both bending and compression is the yield point s stress. For values of. m 


greater | than 1, it gives less conservative results | than Eq. 52. The relation of 


4 
primary to total bending stress is given by Eqs. 56, and use of the chart is 


similar to the method presented for Fig. 2a). 
_ -& “Rational Column Analysis,” by J. A. Van den Broek, The Engineering Journal, December, 1941. 
al 


Eq. assumes gradual change in the total stress as the loading changes 
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By COURTLANDT EATON 
CourtLanpT Eaton,’ 2 —The long strides that have been made 


‘in aerial photography i in recent years are brought out i in this informative paper. 

Of particular interest to the writer are the applications of aerial photography 

to the production of contour maps. About 1930, the Fairchild Aerial Surveys 
under contract w ith the Los Angeles County (California) Flood Control Dis- 
trict made two aerial ; surveys in particularly difficult terrain—one, a ‘a contour 


and west forks; and the other, a contour map of a quarry site i in the lower 
LS canyon about 600 ft high with a practically vertical face. = 


# F rom the San Gabriel Canyon survey some fourteen possible dam sites" 
were located, from whieh, my te field reconnaissance e and office studies, the sites 


for field surveys. Cost data are not available, but. the cost was 
a was possible by ordinary field ‘methods. The c contract called for a high 
en of accuracy and field survey spot check areas indicated that aerial con-— 

tours were more accurate than field surveys. This accuracy was due to the: 


facilities, under a aerial methods, for following 5 up small creeks and draws, which 
in the aggregate add up to substantial areas. — OS 
quarry survey was equally satisfactory a: as” developed by spot checks. 
q 
Because the vertical face made a field survey extremely difficult, | hazardous, 
and costly, it would have involved delays that were inadvisable. 
Aside from costs and time s savings, there are many advantages i in aerial sur- 


Ay 
rious -veys | well known to engineers who have used ape 7 ou advantage i is s that 
ue in aerial surveys can questions of accuracy 


_sraphie re ene such claims either can be reconciled « or Ww will net be made. aa 


—This paper by Leon T. Eliel was published in May, 1949, Proceedings. — I es ae 
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EATON ON PHOTOGRAMMETRY Discussions 
Another advantage is that aerial photographs bring o out it details 3 not visible 
:" ground inspections—for instance, rights of way where. photographs give 
- evidences of old watercourses through improved properties. The writer re- 
members an aerial survey of lands 1 now w occupied by the Santa Anita (Calif. ) 
‘race track which showed evidences of an 1 early race track said to be one built 
by Lucky Baldwin and later changed i in direction for better sun conditions. — 
_--- Sinee there. are also many other advantages i in aerial surveys, topographic 


¢Z veys by present archaic and laborious methods probably will become obsolete 


for major areas, 
_ The writer’s only “criticism” of ‘the pap paper is that it is too short. In his 
closure, it is hoped that Mr. Eliel will elaborate on this most interesting g subject. 
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HARLES | JAEGER, | AND ) Louis: 


CHARLES research work reported in this paper is to be 
commended. because of its general interest. is 6The equation of momentum is 
one of the most powerful devices for computing uniform and nonuniform flow, 
and any effort to check the theoretical results by laboratory experiment should - 
be supported and encouraged. ‘The principal results of the study are sum- 
marized in Figs. 3 and 7. _ Fig. 3 and Eq. 15 are in the form presented by the | 
writer in 1936.33 “Assuming that Eq. .15 can be written directly for for profiles 
(1) and (3), i in Fig. 1(¢) (assuming 


vous 


Substituting a’ 3 = 2 an 0 in Eq. result gives a momentum 
equation similar to Eq. 15. The factor a’; & (a’; > 1) takes account of the 


nonuniform velocity distribution and tz h accounts for a pressure ‘distribu- 

tion that may be different from the hydrostatic pressure on the sill or weir— 


being a coefficient. Since m=> =; an nd Eq. 27 


2” 
becomes 


_ Nore.—This paper by John W. Forster and Raymond A. Skrinde was published in April, 1949, 7 = 
Proceedings. - Discussion on this paper has appeared in Proceedings, as follows: September, 1949, by En- 
§pecial Lecturer on Technical Hydraulics and Hydr. Structures, Imperial College, London, England 
ons. Engr., English Electric Co., London, England. 


18‘‘Der Mischungsvorgang bei ploetzlichem Querschnittsuebergang,” by Charles Jaeger, 


4 

| 

see eee eee ee (28a) 

| 


JAEGER ON HYDRAULIC Discussions 

curve 

The curves ¢; = - = in function n2 are represented in Fig. 11 which is com- J gtatic- 

Dey to Fig. 3. A straight line in Fig. 3 becomes a parabola in Fig. 11; the - 

urve > = In ig. curve =— in Fi an urve —- = 0 in cas 

-eurve = is the «Fig. 1 = ea,” | incas 


hydra. jump) i in . Fig. is the = =01 in Fig. ll. Fig. 12, | (0) 

of a negative sill, should be interpreted in a similar manner, as can also the i 

—- of river bed enlargement. 4 Figs. 11 and 12 represent all possible 
ons of Eq. 27, for positive as w wis as for negative sills. Fig. 3 represents . 


solutions o or nega 


Values of 


| { 


— 
= 
only part « of this range e of solutions. It i is a satisfaction to the writer thus to 
4 find his own studies confirmed by more recent laboratory experiments; and he 


‘is in agreement physical interpretations presented | tall 


e (a) In Fig. 7 ; the experimental | curves are offset somewhat to the right of 


the theoretical curves. _ The authors attribute this offset effect to the non- ‘ 

uniform distribution of velocities. On the other hand if a’, > 1 the 
quantities pi and F actually decrease and the theoretical curve . will then, in me 
fact, move to the left. The writer has found that we ‘1 and a’ 3 | are always very 
nearly equal to 1 (One does not need to know the value of as in n profile ( (2), y 


= to meet the experimental curve. To the writer, this translation of the 


44 “*Technische Hydraulik,’”’ by Charles Jaeger, Birkhiuser, Basel, 1949, pp. 150, 152, 158, and 159, 
Fi igs. 88, 89, 94, and 95 95, Tespectively. 


| 
“Fig. 1(c). ) If & > 0 (as it probably is), the theoretical curve is moved to the 
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curve means that the pressure distribution on the vertical sill is not hydro-| 
static—it cannot be hydrostatic because of the necessity of curving the stream 
‘tubes and because of the turbulence ne near the sill. In some experiments by 
students at the Imperial College, London, England, values of ~ > 0 were found - 
in cases where the occurrence of bridge piers narrowed the section, = 
-(b) There are other ¢ parts o of Figs. 11 and 12 that are not represented in 
Fig. 3. These characteristics have been investigated, thus indicating that real 
hydraulic solutions are available: for these problems as well. A. K Koch and- 
M. Carstanjen'* have done some research on the problem. of streaming flow a 


“over a positive sill and . L. de’® spent a time Anvestigating conditions 


Values of 
_— 

W hen analyzing the energy of the flow that is represented by some of ‘the 
other solutions of these » diagrams s, the investigator finds that ‘some parts of 
the curves are not adaptable to real hydraulic solutions. 
© There is need for further research work to explore the possibility 
further solutions for the different curves in ‘Figs. 11 and 12. : To o date, ole 
parts of the entire have been investigated. example, m more accurate 
measurements of a’ , and will probably show an even closer agreement 
with theory. On the whole, the results found by the momentum equation 


“Von der Bewegung des Wassers und den dabei auftretenden Kriften,”’ A. Koch and M. 
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have agreed remarkably well with the results of tests, an n of these Le 
tests: is still indicated. unusu 


The momentum equation was introduced into hydraulics by J. J. B. Bélanger" “mend 
in 1849 and J. A. Bresse!® in 1860. They established a a comprehensive e e theory, sentec 


introducing the function— | 


—or “total momentum.” Experimental research f followed by H. Bazin, A. H. rise 
Gibson, Mr. Koch, J. -Einwachter, J. Smetana, Mr. Escande, Nebbia, data 
Boris A. Bakhmeteff, Hon. M. ASCE, _and K. W Voycicki investi- § 

‘gated the case of the drowned flow below a gate and H. Favre studied the sag 
lateral spillway and the lateral flow of water into a primary stream. The “no gr 
equation of has also been used for the computation of nonuniform 


flow (translation w. waves). Model tests have demonstrated how the height of 
; the wave crest depends not only on the variation of the discharge (as in the 


ease of the wave e body) but also on the variation of impulse. No wen, i 


"research work has been done on this interesting problem. 
7 One important case that was solved o on . similar Ii lines by Eq. 8i is the compu- 


agree! 

~ tation of partly working surge . tanks, as proposed in 1930 for the ] power ‘station tion. 

ine W ettingen, Switzerland. 21 In that particular case, the computations for - 

_ translation waves had to be combined with surge tank computations. These 

are some of the problems w whose solution by means of the momentum equation — 
represents a considerable amount of successful research work, _ * 

Figs. 11 and 12 also indicate the particular « domain where the « equation of | 7 
momentum offers no solution. Physically, this means that the hydraulic j jump 


becomes unstable and moves upstream or downstream. Similarly, there are 


certain areas in which the Bernoulli equation does not apply. Are these two. ae 

_ domains, w ith their imaginary solutions, the same? What is the connection s 


7 led to some further r research as reported by the writer in 1931, 2 1943.23 1947,4 
1948, and 1949.2 
‘When the investigator is s dealing ig experimentally with particular problems 
¥ outlined in this discussion under topic (c), the 1 more general aspect of the 
equation of momentum described in the paragraphs following topic (c) should 
+¢ kept i in mind because it is there that the explanation yn of many of the diffi- 


encountered isto befound. 


17 ‘‘Notes sur le cours hydraulique,”’ Bélanger, Mémoires, Ecole Nationale des Ponts et 


de Mécanique Pt. 2, “Hydraulique,” by J. A. Bresse, Mallet-Bachelier, Paris, 


19 ‘Technische Hydraulik,” ” by Chesies Jaeger, Birkhauser, Basel, 1949, pp. to 


- between ‘ “total momentum” and energy? These theoretical problems have 


2 ‘Uber die Eigenschaften von Schwiillen und die Berechnung von Unterwasserstollen,”” by Meyer- 4 he 
and H. Favre, Schweizerische Bauzeitung, Nos. 4 and 5, 1932. show] 
“‘Analogien zwischen Stiitzkraftminimum und Energieminimum in der Hydraulik,” Y by Charles witho 
‘Contribution a l'étude des courants liquides & surface libre,’ by Charles Jaeger, Revue Contre de that t 
**Te limpulsion totale et de ses rapports avec totale d’un courant liquide a libre,” 
by Charles Jaeger, ibid., Nos. 37-41, 1947. beab 
Steady Flow in "Open Channels: The of Boussinesa,’ by ‘Charles Journal, Inst. ea 


_C. E., February, 1948, p. 338. 
‘Technische Hydraulik,” by Charles Jaeger, Birkhauser, Basel, 1949. 
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M. Jun. ASCE. results of this paper present 
unusually. good data on the hydraulic - jump, and the authors are to be com- 
mended for the manner in which they conducted the tests. The curves pre- 
sented lead the closer to rational of The 


‘more complete analysis than n that given in n the paper. 


4 - — Abrupt Rise.—The authors state that thee experimental curves for the abruy upt 


rise (Fig. 7) are parallel, straight lines. 7 An equation to fit. the | - experimental | 
data was easily found by plotting d3/d, against h/di, with the parameter being 


straight- line curves of values of F all having approximately the same oe 


‘The resulting equation fits the experimental curves within + 5%, and had 


‘no greater er ror than the spread of the rege! points: ds = dz — 1.40 40h; or 
ds/dy = do/dy — 1.40 = (F) — (h/d:). To complete Eq. 14—- 


] 1 ay 
—which includes the “necessary variables F, dy and “Fig. 13 shows the 
“ @ agreement between the authors’ experimental « curves and the empirical equa-_ 
or 
on Author’s Experimental Curves, Fig. 
Points from Equation 3 30, 
on 
> 
ve 
ms 
13. —V ERIFICATION OF AL Equation with EXPERIMENT AL au Dara oF Fic. 7 
ris, 
(dy. +h) = 0.40 h = Ad.. . -(31) 
showing t the resulting drop i in \ the water surface below that required for a jump 
Without an abrupt rise. se. Within the limits of the authors’ studies, ‘it follows 7 
that the ratio h/Ad = 2. 5. It may prove to be a useful rule of thumb that the 7 


drop in the v water surface below that level required theoretically. for dz would 


be about 40% of the height of the abrupt rise. 


™ Asst. Prof. of Civ. Eng., Carnegie Inst. of Technology, Pittsburgh, Pa. 
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_ Fig. 14 is another plotting of the authors’ data’ (based on the ec curves in 
7) to show this drop. The care with w which the experimental work must 

have been done i is evident. Denoting K as the ratio (ds + h)/ ds, in 1 which dy 

is the depth computed from Eq. 4, these observations follow —_ ate 

() af ratio K is decreased by we h/dy (F constant) or by decreasing 


@) TI ‘The j jump and the drop i in the water are little affected by the 


Fishigh; 
(3) whe ratio K K approaches 1. 1. 00 for all values of F and h r/ dy when F becomes 


(4) If h were increased (F and d, constant), ds; must be lowered if the j jump 

not to move i in the basin; and 

5) An empirical equation for K follows: 
+>=1- 


2b) 


designer might consider the ratio 0 of importance when the tailwater 
is low and a depressed basin is considered. | ~The ratio of the depth of tail- 
= which may be used compared | to the depth dz is computed for the plain 


jump. Iti is not to be consider ed a a new - variable, but may be useful in design 


as shown in the form of Fig. 1 15. Calling on the empirical data, 


_ Fig. 15 shows the agreement between the experimental data and re curves: 
computed from Eq. 33a. This was facilitated by choosing d= =1, values: 
being taken directly from Fig. 7 or being | computed as a function of FRO 


~The difference between the theoretical and experimental curves 
plained by the authors to be due to a nonuniform. velocity distribution w hich 

at the ratio X/(h + = 5. They a are correct in stating that, w hen 
X/aAt+ ds) approaches 10, the mami data will check closer with the 


_ theoretical curves. However, nonunifor m Vv velocity distribution will not ex- 


plain all the discrepancy. 
_ The omission of a dynamic force on the face of the abrupt rise is not pel- 


missible. Although the > force i is small i in some cases, it it becomes exceedingly 
important when the j jump is close to the rise. 
As long 4 as X i is greater than the e length « of ‘the j jump, the assumption | can 


be made that the jump i is not ¢ changed by the 1 rise follows. Thus, 1s, the 


rey 
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K = 1—7 — = a, F)..... 
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0 
_d3 = do — 1.40h or — = 1 — 1.40> = K’..........(33a) 
or tl 
eque 
prov 
a 
— 


October, 1949 -LAUSHEY ON HYDR 


a lower level. on the | rise. This problem can also be analyzed easily by the 


energy following basic. equation 1 applies a: as long as X>Lys: 
29 4 a ds 
_ kinetic energy; and K is the coefficient for contraction loss at the rise. In 7 
fact, the exact theoretical curves can be duplicated by Kg. 34 ifa =1 and 


= 0. The writer computed a half dozen points which checked exactly. 


ol 


4 


d3 


> 


m Fig. 
4 Values Computed 
Empirical Equation 33(a) 
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Values of F * of F 

Fig. 14.—PERCENTAGE Drop IN THE WATER SuR- _ Fic. 15.—VeRIFICATION OF EXPERIMENTAL Data © 

Face THE ABRUPT RISE AS Al Func- AND Empirical EquaTIon FoR DEPTHS ON 


TION OF Fanp H/di AND BEFORE THE ABRUPT RISE 


‘Iti is a basic fact that the hydraulic jump can be analyzed by the energy 
or the momentum principles provided the head | loss is known when the energy 
equation is applied. Thinking m may be less confusing if both a are. applied to 
‘provide a a better understanding « of the problem. However, in both the energy | 
equation | and the momentum principles, it is imperative to know the velocity 
distribution to compute, correctly, the true momentum and energy at point 3, ‘ 
thus ¢ deriving an | equation that checks experimental data. 
‘The abrupt r rise can easily be analyzed by the are wing equation as long 
4 aS long: 


4 @tk) at + K) 
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 : he writer computed fifteen | values of a-+ K from the experimental data. 


a This coefficient averaged 5/3 and varied very little from this value, which is 


a reasonable one. In Fig. 16 points are shown computed from Eq. 35, w hich 


"show the close agreement, | w hen a K is to 5/% 


ate 
Values of the 
corr 


prin 
ve for 
men 
of F 

1G. 16. BETWEEN Curves, 7, AND Mopirie 

_ Eq. 15, based on the momentum principle, must also be modified to include 
I ple, 
7 _ the « correct coefficient 1 for momentum at point 3 and, therefore, E Eq. 1 
be e modified as follows: 


2 \ P 

C 
_ 


ad (37) 


_ Actually Cp varies, but the —— is small within the experimental range, > 
the variation would produce a negligible difference. average value of 1.25 

= computed from the data for a ’. However, it varied considerably more> 
- than a + K, and it seemed to be a function of the Froude number. The 


variation of average values of a’ with F is as follows: 
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‘The use of the proper value of a’ ‘and of the ens force on the rise 
corrects Eq. ‘15. By using values of a = 1.25 and_ of the drag force with 
Cp= 2.0 0, points | were computed : and they are shown for comparison i in mm 16 
Peel 
course, the correct coefficients would nee exact agreement. 


Bewever, the 
be than) ‘at section 3—either the introduced 
could provide the lump correction, or further investigation of velocity dis- 


tribution could provide more exact values. 

_ The writer has emphasized the fact that, when the abrupt rise is placed 
after the j jump, the problem i is simple and need not be confused directly with — 

the hydraulic jump. If F and. are ‘known, then and V2 follow, and the 
correct drop dow! n can be computed easily by either momentum or or energy 


principles provided the velocity distribution and the contraction at the 


weir should hangeable. is important to si superimpose he experi-— 
mental curves of Fig. 7 at The agreement 
between two is in Fig. 7 where it is ‘noted that the values for the 
‘ratio 


dy 


Values of 


N 


Val f 


17. OF THE RISE AT WITH THE Fie. 6 
The two theoretical curves (at X/d2 equals about 10), Fig. 2, also agree’ 
reasonably well. W ‘hy: should not curves for the weir and abrupt rise at critical 
depth be the same for other ratios, as long as x > La The ‘same barrier 
should: cause approximately the same velocity distribution. and force on 
- the rise, irrespective of what happens immediately after: the barrier . These 
beg geal on the close correlation between the limiting | case of the abrupt 
rise, ds = d,, and the weir suggest that one analysis (the abrupt rise) should 


‘Suffice for both. , The: main criterion is s whether the flow is above « or at or — 
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critical depth at the rise. the analysis should be the same, 
a although slight variations occur due to the difference between 
the length of the rise in the two cases. «dt would be most acceptable if if the 


limiting case of the abrupt rise would parry the terrifying Eq. 18a, based no 
a. estionable coefficients, for the weir. . Nothing should be that complicated. 
nut lth, The weir analysis is unnecessary, therefore, except | to show the variation in 
hai in terms of x/ds fc for values of F. The he authors state that, when the dow n 
stream slope i is steep al and the depth i is less than critical, a weir or abrupt L rise 
_ of the same height a are re interchangeable. — ‘Fig. 17 supports this statement. —— 
An attempt to predict the change in the experimental curves of the abrupt 
‘tise (Fig. 7) when (ds + h) is not equal to 5, is possible by combining the 
data f for ‘the weir with the abrupt rise. A first approximation for similar” 
curves « of other ratios of a/(ds- + h) in Fig. 7 can be obtained as follows: 
- Project points from Fig. 6 to the line ds =d0 of Fig. 7... ; Then rough curves 
can be sketched in parallel to the existing ones of Fig. 7. 
Additional Variables. —The fact that an -backwater curve ‘exists 
stream from the j jump n may be the cause of some error if the depth d; is s measured 
at a constant a" upstream from the rise, which may be upstream from the 
toe of the j jump. an As a matter of interest, if there were no backwater curves 
upstream and downstream from t the j jump, it could probably ta take | place at any- 
; i here : along the channel and would likely be unstable as to position. — To: move 


jump upstream, the dow stream depth must be of this 


a= 


discharge being If is measured ‘the depth at 
7 the toe of the j jump will be higher than expected and the Froude number should — 
really be lower. _ This condition, combined with a slightly lower h/d; than oe 

puted, would move the experimental points toward closer agreement with the | 
curves. Hence, for precise work the depth 1 should be measured 4 

_ the toe of the j jump, especially when the ratio X/d2 varies considerably. es - Atany 


= the backw ater curve is a a ‘separate emg w oa solution can be handled 


‘The ratio. do/d, must be as a function of Ordi- 
or V4/(g dy) will be known by ihe engineer. A rapid method for 
‘approximating dh/dy * = f (F?) mentally i is as follows: Set down values of ds/d, 
| from 0 to n. Values of F* are then tabulated quic 


F*-value— 


Pe 


: ™ value of elds wi will be the difference bet aia the preceding a and the existing 
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«dt ‘would be of interest if the authors would state whether they wea 
‘ds, which was assumed to exist before the rise, to see if it agreed with the — 
computed value. In other words, is the j jump preceding a | barrier (as long as as 


X > Ly) any different in dimensions or energy loss from a similar jump in a 


levelchannel? 


is The methods presented for the abrupt rise and weir cannot be used vi when 
the X-distance is less than the length of the j jump. . The theoretical depth d2 ; 
would certainly not exist upstream from the rise, and possibly not at all. it 


would be a difficult problem. The momentum equation would be most ap- 


plicable, but a study of the velocity distribution, dynamic force on the face al 7 
the rise, and assumption of hydrostatic pressure would be absolutely necessary. 


_ Research ¢ of this kind may have value in the design « of even shorter satisfactory 


stilling basins in certain instances. 


os A complete s study « of velocity distribution i is necessary in the problem pre- — 
sented by the authors. As long as > Ly, it should be possible to express 


the coefficients a, a’ , K, and Cp as functions of X/d, which would complete 


the problem by oe considering the ) placement of the jump with respect to 7 
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ISSIPPI RIVER BRIDGE AT DUBUQUE, : IOWA 


BYR.N. BERGENDOFF AND JOSEF SORKIN 


RN. BERGENDOFF® a AND JosEF SorKIN,®*° MEMBERS, ASCE.— The p: primary 
= of presenting | this paper, giving a detailed description of the design 


and construction of one of the major structures across the Mississippi River, was 

: to » establish a . permanent record regarding this bridge in the belief that it may 

prove helpful to colleagues in the profession, in the planning of structures of 

comparable magnitude. Ina sense the urge to present this paper was a reflec- 

tion of the writers’ own experiences and difficulties | encountered in the usual 
7 exploratory work, in planning major bridge crossings and in pursuit of deter-— 

- mining best suitable types and features for a given site. There is an abundance 

of theory on design analysis of of every | known type of structure anc and, although it 


is valuable i in the final design | of a structure, theory alone « offers little help in| 
studies of comparative types, particularly when the time element 


is limited as has been: the usual case in the past several years. 
~ The writers were particularly pleased, therefore, i in the concurrence with this 
logic as indicated in the opening paragraphs of Mr. Jones’ discussion. It is 
also noted that Mr. Jones concurs in the pleasing appearance of the bridge. 
In this instance the structure as located is in full view of the inhabitants of the 
| cities served on on both banks of the river as well as to railway traffic a across the 
river and along the east and west banks. The designers felt particularly 
- gratified when the structure received the 1943 Bridge Award of the American ene 
7 Institute of Steel Construction in Class I for the most beautiful bridge ge opened 7 
to traffic in that year, 
4 Mr. Jones: questions the inherent economic advantage in the of 
this | type of bridge | as compared to. the conventional cantilever or or ‘continuous 


= construction. He ‘concedes, however, that “Unfor tunately, no o data can ri 
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Even bids received on alternative types, although | vn 
nee serve to indicate the economic preference in one particular | case, would 
not permit drawing general conclusions, as s unfortunately too many elements 
‘significant | at the moment may be totally irrelevant at some other time, « or ata 
different geographical location. It is well recognized that, aside from such 
fundamental costs as those of material and labor, the contractors’ bids also 
reflect in no small measure the lack or abundance of prospective work, the 
chance location of the major part of their equipment necessary to the job, and 


various other factors of similar nature. ey 


. In the design of the bridge a 50% increase in stresses was allowed - 
erection including wind.  Tothe best knowledge of the writers, such an increase 
is conventional. 7 Mr. Jones’ statement that “This increase was the occasion 
for considerable worry + + Ww) is fully appreciated; but it is believed. that the 

‘structure through various construction stages was not devoid of ample ss safety 
factors. . The erection stresses, although higher than normal, were still well 
below the elastic limit of the metal, and the structure was sectioned for a 30-lb 
wind against a maximum recorded wind v elocity equivalent to an 11-lb wind. 

Le The writers cannot agree with Mr. Jones’ deduction that in conventional 
‘cantilever and continuous truss construction by cantilever method ‘mem- 
bers carry the moment from cantilevering only the half length of the — 
span; instead of h half | the length of the main n span, as at Dubuque.” » In canti- 
lever erection of a cantilever truss the temporary chord members, introduced - 
at the hinge of of the suspended span, , necessarily produce a moment at this point, 
which carries s along the anchor arm to the support; in essence then, the final 
“result i is the same as for that of the Dubuque bridge. For purposes of record 
it may be of interest to note that the total amount of metal added for erection 
“of the Dubuque bridge i is 150 tons (not 345 tons, 8, a8 erroneously : stated in Mr. 

Balog’s discussion) which in knowledge of the writers compares quite 
favorably with conventional cantilever or continuous bridges of similar ‘Mmagni- 
tude, ' The Mississippi River Bridge : at Greenville, Miss., gree to traffic in 
1940 and designed a also by the he writers’ firm has approximately the same length - 
of channel span, but i is of conventional cantilever type. The final weight of 
_ metal added for erection purposes is comparable with that in the Dubuque span. 7 
Py Mr. Floris would apparently prefer to use deck truss construction for the | 


* by making th the top chord of the end spans straight, while mie 


“thei bottom chord. This arrangement permits the lowering of the mid- 
piers, if the suspended floor of the center span is also raised to the level of 


the top chord of the end spans.” 
Mr. ] Floris states that such an arrangement w would nome in ‘Additional naviga- : 


tion increase in cost is not ¢ juite clear 


“clearance. “Additional horizontal clearance | would necessarily r mean an increase 
‘in the bengthh of the center s span. - Additional vertical clearance would r require 
either an appreciable steepening of the approach grades or an increase in the 
“height 0 of approach, span piers, or both; in any event there i is reason to believe 


construction would result in ‘additional cost. Furthermore, vertical 
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Se clearance requirements under the west end span eliminated the possibility of 
such an alternative in this scase. The | writers appreciate Mr. Floris’ reference 

_ to three publications in German as regards other applicable methods iy! stress 
analysis. Since, however, Mr. Floris explains “These methods * * have 

* * * little favor among authors and stress analysts” in the United 
States,” ” the writers are pleased to to learn that they alone did not miss the 
= rtunity of using these methods. 7 Mr. Floris could make a real contribution 

to the profession if he would translate these publications into English. The 

_ customary procedure used i in the design | was found 1 neither “cumbersome” nor 
“laborious.” All basic data. preliminary to the development of Eqs. 17a, 17, 
_ and 17c an analyzed by slide rule without any significant loss of accuracy. 
_ The solution of the three simultaneous equations by any accurate method 
certainly cannot be classed a as laborious. / It is the writers’ considered opinion 

3 that in most instances t the difference between one method of analysis and num- 
-erous others which are generally available is not. material. ‘The best method 

_ for any individual analyst is the one that best appeals to his logic. The time 

i element, if any, is usually « of secondary importance and i is generally compensated 


by elimination of e1 errors, , because of the analyst’s ¢ s clear understanding of the 


theory, 


_ Mr. Rust ‘emphasizes the particular importance of stable foundations for 


= greater t than that required, averaging well over r twice the Sorat 

capacity. ‘The writers fully agree that statically indeterminate structures’ 

ee may well be avoided in instances where considerable settlement may be antici-_ 

7 pated. Its should be recognized, however, that nominal settlements are of no 

material significance. As was s pointed out in the p: paper, a settlement’ of 1 in. in 

one of the end piers would result in a change of the reaction at that pier on the 

~ order of 10 kips, which i is small, “‘percentagewise’ ’ and its effect upon the str esses 

in the truss members would be equally insignificant. - Similarly, the same 

settlement of the intermediate piers would be of no great importance. Actual 

observations, since completion of the structure, ‘hee disclosed no evidence ot 

: Mr. Dohn Taises the question of secondary stresses in the truss members, 

particularly near the intermediate piers. Contrary to Mr. Dohn’s presump- 

tion there are few if any unusual features in the framing | which would cause 

secondary stresses of "extraordinary magnitude. —«dt is also difficult to follow 

Mr. Dohn’s deduction that the use of silicon steel v will necessarily result i in 

@ greater secondary stresses. . Geometric cambering to neutralize the secondary 

y 7 ‘stresses was used for the Dubuque bridge. — It is to be appreciated that in long 


—- highway bridge ae stresses due to live load. are re relatively > small 


| _ supports; and, therefore, ‘the effect. of of live loading is relatively ‘small. As 
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of primary loads a very few men mbers are subject to agente re to wind. : 
Mr. Balog’ s protracted discussion ‘Tepresents | a great deal of f effort and 
time. | ‘Such effort could be a distinct contribution to the original ps paper. Un 
fortunately, Mr. Balog’ s comments are so thoroughly misapplied as to merely _ 
confuse the factual data presented i in the original paper. The discussion can 
be criticized on the following elemental points of recognized sound saiaeeenatinal 


|. Construetion costs are not, and never hav re constant. 


Minimum “amount t of material does not necessarily represent - the most 7 
economical design. . If so, wide flange rc rolled beams would have no place in 


3. Foundations on rock available at. nominal depth can be e designed « cheaper | 


than foundations necessarily supported on long friction piles. i oe 
4 . The « cost per linear foot of a span 845 f ft long is more than the cost per — 


linear foot of a span 740 ft long. 


. 5. The ability to compute stresses and strains is not the sole qualification 

| ‘of a good designer. Knowledge backed by practical experience with problems 
dealing» with properties of materials, fabrication of parts, function of details, - 
limitations ns and adaptability of of erection equipment, costs of material and labor, 
restrictions of War Department requirements, etc., is necessary to permit the 
designer to solve his problems properly. 


«6. In any design, due consideration must be taken of relative costs of labor 


and materials. . European practice favors more labor and less material because 
labor i is cheaper in Europe; the reverse is true in American practice. ae oe 
Mr. Balog makes an extremely questionable co comparison son of the cost of the , 
Dubuque bridge with that of the Bettendorf (Iowa) bridge. After arriving at 
a misleading difference in the cost per square foot of roadway by ignoring 7 
elemental facts, Mr. Balog accepts it as prima- -facie evidence for his thesis. — 
The Bettendorf bridge was constructed in 1935 when the Engineering pil 
Record Construction Cost Index was 196.44. ‘The Dubuque bridge was 8 con- 
structed. at 1941 costs when the Engineering News-Record Construction Cost 
Index was 257. 84, a difference of 31.5%. The main span of the Dubuque 
bridge is 845.17 ft, whereas that of the Bettendorf bridge i is 740 ft. An in- 
crease of of the Bettendorf span to the length of the Dubuque span would result 
in approximately a 15% increase in the pull on the cables, height and load on 
| the towers, and in the pull on the : anchorages. — Furthermore, the entire Betten- 
dorf bridge, including approaches, is | founded on shallow ‘Tock foundations. 
At Dubuque subsurface material for a depth of 400 ft is sand, silt, and clay. 
It is obvious that, transposed to the Dubuque site, the Bettendorf bridge 
would d entail materially greater foundation costs. On the other hand, founda- 
tion costs for the Dubuque bridge would be if that 
structure were transposed to the Bettendorf site. 
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liate @@ paper, that the bridge, as finally constructed, compared favorably in cost to a _ a a 
suspension span. Mr. Balog overlooks these factors altogether and seemingly 
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bases his conc conclusions | as regards s similarity of the two bridges on the mere fact | mitte 
that the ratio of cost of substructure to superstructure is 0.52 in the case off econc 
_ either b bridge. The identity of these ratios is purely accidental, and has § truss 
nothing to do with comparative economy of the structures. Misuse of per- as W: 
_fectly good data is ‘no Teflection on on the Bettendorf bridge. , This bridge was § conce 
a "designed by e eminently | able engineers ; and built by competent contractors. to m 
The type and other features of this bridge were undoubtedly dictated by sound engin 


Mr. Balog takes Seaton issue w with comparative ec economy studies of the divul 

‘Dubuque bridge, stating: thus:s engin 

torsic 

: _ “The cantilever, being a statically determinate system, cannot utilize the | the oF 
= tenacity and ductility of the steel ; ** * and * * * must necessarily § 

cost more than the continuous truss.” | of su 


Although this argument is applicable to short s span n structures, i in which the live § struc 
load f forms a significant part of the total load, it is totally ‘irrelevant i in spans struc 
of the magnitude under consideration. "(For the Dubuque bridge, the T 
2 load is 15% of the total load.) In long spans, introduction of articulated § ough! 
joints at inflection points of of the continuous truss span, | thus converting it into spans 
— a cantilever system, will not in the least detract from the economic aspects bal 85-ft 
ductility” which in Mr. Balog’ opinion are inherent in con- § unbr: 
tinuous systems only, 
Mr. Balog is concerned with the position of the mass centroid i in the Du- § of tk 
--_ buque bridge. _ However, since later in his discussion he takes i issue with the mem 
«Spacing of the trusses and maintains that the trusses could just as well be & cant 
spaced closer, he apparently i is quite satisfied with the stability | of the bridge. | = 
Thus o one part of his discussion answers the other. His argument relative to tedly 


the mass centroid of the structure, if carried to its logical conclusion, would § were 
—- through braced arch, including the Hellgate Bridge (New § to th 
= . Y.), Sydney (Australia) Harbor Bridge, : and | many other equally § of ¢p 
"prominent structures. — It also leads to the implication that, since the Empire differ 
Building i in New York could not possibly have the stability of a single: que 


7 story frame, skyscrapers are follies. To emphasize his point Mr. Balog evel Bi frst. 


_ cites the dangers of earthquakes. Such logic has no substance. © ba question 575 | 
‘is not which type of structure is infinitely more stable, but rather, “Is the by ¢ 


‘structure stable?” 


By the way of further comparison, Mr. Balog p. proposes two other types 0 


structures in addition to the suspension bridge as shown in Fig. 10. 5 As regards fern 
oe the bridge indicated i in Fig. 10(a), a structure quite ; similar was actually studied ] th , 
for the Dubuque bridge. — Owing to the span length involved, 1 no adva antage | 
was found in this type of construction and it was dismissed primarily due to 

- anticipated erection difficulties. — Because of the similarity of this system | tothe @ — 
one finally adopted, the writers did not deem it pertinent to mention this study —_ 

in their original paper. Subsequent to the Dubuque project, consideration — 
7 ‘was again given to the use of a girder type tied arch at another site where the amt 
span was materially shorter. Designs of a trussed arch similar to jo that at at with 
Dubuque and of a girder arch similar to that shown in Fig. 10(a) were sub- ord 
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fact | mitted to competent fabricators and erectors for their views as regards relative _ 


se of J economy of the two structures. Their conclusions, too, were in favor of the — 
has | trussed arch, , principally because of necessary ‘maintenance of channel traffic, 
per- a8 Was also the case at Dubuque. — Erection matters apparently are of no 
was ff concern to Mr. Balog. _ To state “Erection considerations, being subject also. 7 
tors. § to many unforeseen circumstances, should not govern the design,” is for the 
ound | engineer to adopt the proverbial habit of the ostrich, hiding his head in the - 
_ ‘sand hoping that the contractor will complete the job ey cost which will not om 
f the divulge the ignorance of the engineer. Evidently Mr. Balog envisages bridge 
- “engineering as simply a matter of wrangling with “tenacity and ductility,” 
__ Ei torsion, stresses, and strains and throwing the mathematical brain child into 
oti the omniscient hands of the contractor, and hoping for the | best. The outcome 
of such procedure i is obvious. contractor may either build the structure 
7 at a cost fully reflecting all the follies of the design or he may redesign the 
elive structure where any similarity between the original | design and the final 
spans structure is purely coincidental. 
live The. continuous { truss unit propounded by Mr. Ba Balog i in in Fig. 10 ‘10(c) is — 
lated oughly contrary to all recognized concepts of economic proportioning of truss 
; into. spans. _ The writers cannot agree with Mr. Balog that a truss with 85-ft panels, 
cts of 
type In this case again the mere weight | of metal is only a part 
» Du- § of the problem. The cost of handling and — transporting the long, — 
h the § members and the field costs of splicing and erecting them are clearly, of signifi- 
all be Bcantimportance, 
ridge. § As regards the design of the stringers, Mr Balog misinterpreted the admit-— 
ve to tedly ambiguous statement in the paper. Actually, although the crossbeams 
vould were designed — as continuous beams" over the stringers, the load distribution 
(New to the e stringers ¥ was substantially i in accordance with the 1935 specifications 
ally of the American Association of State Highway Officials (A.A. S.H.O. ). The 
mpire differences in weight arrived at by Mr. Balog as indicated in Table | 8 are 
ingle- due to serious oversights of certain factors and effects by Mr. Balog. In oe 
even & first place, the weight of the stringers s and crossbeams is 530 lb per ft instead of 
— 575 Ib. _ Furthermore, half of the sidewalk including its live load is supported 
wd by th the stringer adjacent thereto. The floor systems indicated by Mr. Balog 
_ as type (a), , Table 8, ignore this load completely. Similarly, the type (6) floor, 
ok ‘Table 8 (Mr. Balog’ 8 design n), would obviously require a thicker slab. — “The 
udied additional cost of concrete, reinforcing steel, and the cost of additional metal in 
ntege the floor beams and trusses to carry the « extra dead load would more than offset 
ue to the gain in 1 stringer weight. ‘The same applies i in varying gg degree to type (c) 
+ the and type (d) floors, from Table 8, where i in addition to these factors, Mr. Balog | 
study | fails to account for diaphragms s required by ‘specifications and for the additional 
ation weight of connection angles as as well as for the cost of additional field riveting. 
“ the i [herefore, Mr. Balog’s comparisons are irrelevant since they lack common basis 
at at "a the 1 writers’ design, and i ignore de details considered essential by engineers, — 
sub- design as well as for construction. Mr. Balog’ 's conclusions stated ‘thus, 
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‘comparative ‘investigation of the type of floor system 
= by the authors and a five stri inger arrangement. * and even then 
type (a), Table 8, cost more,’ oo 


do not correspond to the numerous studies made by the writers. In fact, “the 


writers have found that designed on a common basis, with all p pertinen nt elements 
properly appraised, al all | conventional floor systems result in the ‘same cost for 
all practical purposes. — As regards the use of silicon steel for stringers, Mr. 
Balog fails to account for the difference in cost as compared to carbon steel. 


—_ from the difference i in ‘direct costs, , use e of silicon: steel for floor systems 


‘fore subjected to numerous stress cycles, is s seriously questioned by many engi- 
neers. In this connection it is pertinent to note the specifications of the Ameri- 


can Railway Engineering Association thus, 
“On account of the possibility of fatigue failure, the allowable unit stresses 


in Article 1801 shall not apply to members that will be subjected to a 
large number of stress cycles during the life of the bridge and having either 
reversal of stress or low dead load stress. For such members the allowable 
- unit stresses shall be reduced, and in floor systems shall not exceed those 


Specified i in Article 301 (for cai carbon ‘bon steel).” 


_ In ¢ discussing ¢ the lateral and sway ay bracing Mr. Balog « cites several European 
Spano where sway bracing an and lateral bracing : at the lower arch chord w ere 
d. Undoubtedly, however, this is only a part of the story. Cross. 

frames ean can be eliminated if strut type portals are provided i in lieu thereof. It 

is entirely possible that such portals require less metal. _ However, the thing of 

utmost importance is not the type but the ultimate cost thereof. — Tti is well 9 

recognized that European practice conserves metal by taking peer of 

4 the cheaper labor. _ This, in fact, demonstrates the universality of the engineer's 


. ‘ingenuity. However, it is is not ot by mere accident that these practices neither 


were originated, nor have found application, in the United States where the a 
‘relative costs of labor and materials differ radically from European’ costs. J Lig 
Indiscriminate « compar ison of these practices i is s absolutely i improper. — 
As a simple observation as to the advantage of unobstructed view, the 


writers | s stated: “The three- -span continuous truss structure * * * affords to 

traffic a a view unobstructed by diagonal truss members across the $45-ft main 

channel span * course, , the writers presumed the pertinent line of 

; sight of from 4.5 ft to. 5 ft | above roadway level. _ This observation is simple, 
factual, and clear. Mr. Balog complains of the obstruction caused by the 
bracing of the arch, which is from 25 ft to 100 ft above the floor. 7 
Mr. Balog claims his alternate designs show 20% less metal. “This is 

contrary to results obtained by the writers in preliminary studies ‘made f for the 
: “crossing and the writers believe such extreme claims are as erroneous as the 

~ basis of the other criticisms of Mr. Balog’s discussion. elie 

a he writers wish to express their appreciation to the constructive contri- 


butions of the discussers, which helped in crystallizing some of the phases 0 of 


described in the paper. 
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